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1.  Statement  of  the  Problem. 

An  investigation  of  the  microscopic  behaviors  of  low  molecular 
weight  sorbates  in  polymers  is  proposed  using  nuclear  magnetic  resonance 
(NMR)  spectroscopy.  The  study  will  focus  on  both  gases  and  condensed 
penetrants.  The  diffusion  of  small  molecules  into  polymers  has  been  the 
subject  of  many  investigations,  because  of  the  importance  in  industrial, 
medical  and  biological  applications.  Some  examples  are  the  membrane 
separation  of  gases,  controlled  drug  delivery,  control  of  polymerization 
and  the  production  of  barrier  materials  for  packaging  and  chemical 
defense.  The  diffusion  of  penetrants  in  polyisobutylene  is  of  special 
interest  since  this  material  displays  low  permeability  for  a  polymeric 
rubber  above  the  glass  transition.  The  choices  of  penetrants  are  for  the 
gas  system,  carbon  dioxide  and  for  the  condensed  sysem,  toluene. 

The  purpose  of  the  investigation  is  to  explicitly  characterize  the 
molecular  details  of  penetrant  transport  in  terms  of  the  structural  and 
dynamic  nature  of  the  sorbed  species.  From  the  polymer  perpective  we 
also  wish  to  characterize  the  effect  of  the  sorbed  species  on  the  polymer 
properties,  both  structural  and  dynamic.  In  this  regard  the  effect  on  the 
glass  transition,  sub-glass  relaxations  and  polymer  integrity  will  be 
important.  Also  the  potential  cooperativity  between  penetrant  motion  and 
polymer  chain  segmental  motion  will  be  addressed.  The  spectroscopic 
approach  here  is  intended  to  complement  the  macroscopic  studies  by  other 
research  groups  including  those  at  the  Materials  Technology  Laboratory 
(Watertown,  Ma.)  and  the  Natick  Army  Laboratory  (Natick,  Ma.).  In  a  general 
sense  the  molecular  level  approach  here  will  help  to  develope 
interpretational  models  to  accommodate  both  macoscopic  and  molecular 
information. 

The  specific  experimental  approaches  used  to  characterize  both 
structure  and  dynamics  in  these  rubber/penetrant  systems  involve 
experiments  measuring  NMR  relaxation  phenomena  and  relating  those 
measurements  to  the  molecular  motion  of  the  species  probed:  either 
polymer  or  penetrant.  This  is  accomplished  using  well  established 
motional  models  for  both  small  molecule  reorientation  and  polymer 
segmental  dynamics.  In  addition  the  NMR  Pulse  Field  Gradient  technique  is 
used  to  directly  measure  the  penetrant  diffusion  constant.  This 
experiment  has  the  advantage  of  not  perturbing  the  system  and  the 
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determination  is  not  model  dependent.  Also  any  complexity  in  the 
diffusion  process  is  revealed  (eg  heterogenious,  anisotropic  or  restricted 
diffusion).  The  timescale  during  which  diffusion  is  measured,  and  thus  the 
lengthscale  probed,  is  at  the  experimenter’s  control. 

Both  the  temperature  and  concentration  dependence  will  be 
measured  and  interpreted.  Current  models  for  molecular  dynamics  will  be 
combined  with  the  concepts  of  Free  Volume  theories  to  characterize  the 
diffusional  behaviours.  The  comprehensive  data  base  that  the  relaxation 
and  field  gradient  measurements  yield  on  all  the  components 
simultaneously  should  provide  insight  into  both  the  applicability  of  the 
theories  used  and  any  cooperativity  between  polymer  and  penetrant  at  the 
local  structure  level. 
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3.  Summary  of  Results. 

The  results  and  conclusions  of  the  project  are  contained  in  the 
following  compilation  of  manuscripts  both  published  and  submitted  for 
publication  as  indicated: 

(i)  “Nuclear  Spin  Dynamics  of  13CO2  Sorbed  in  Polyisobutylene  Rubber” 
Z.P.  Dong,  B.J.  Cauley,  A.  Bandis,  C.W.  Mou,  C.E.  Inglefield,  A.A.  Jones,  P.T. 
Inglefield  and  W.-Y.  Wen.  J.  Polymer  Science.  Poly.  Phys.  Ed.  2li1213 
(1993). 

(ii)  “An  NMR  Study  of  Segmental  Motion  in  Polyisobutylene  and  the 
Relationship  to  Translational  Diffusion  of  Sorbed  CO2”  A.  Bandis,  W.-Y. 
Wen,  E.B.  Jones,  P.  Kaskan,  A.A.  Jones,  P.T.  Inglefield,  and  J.T.  Bendler.  J. 
Polymer  Science.  Poly.  Phys.  Ed.  2Z,  1707  (1994). 

(iii)  “  NMR  Study  of  Segmental  Motion  in  Polyisobutylene  and  the 
Relationship  to  Translational  Diffusion  of  Sorbed  CO2”  A.  Bandis,  W.-Y. 
Wen,  E.B.  Jones,  P.  Kaskan,  A.A.  Jones  and  P.T.  Inglefield.  Polymer 
Preprints,  ACS.  1,  427  (1994) 

(iv)  “A  Nuclear  Magnetic  Resonance  Study  of  Dynamics  in  Toluene- 
Polyisobutylene  Solutions:  (I)  Penetrant  Diffusion  and  Fujita  Theory”  A. 
Bandis,  P.T.  Inglefield,  A.A.  Jones  and  W.-Y.  Wen.  J.  Polymer  Science.  Poly. 
Phys.  Ed.  (submitted). 

(v)  “A  Nuclear  Magnetic  Resonance  Study  of  Dynamics  in  Toluene- 
Polyisobutylene  Solutions:  (II)  Ventras-Duda  Theory”  A.  Bandis,  P.T. 
Inglefield,  A.A.  Jones  and  W.-Y.  Wen.  J.  Polymer  Science.  Poly.  Phys.  Ed. 
(submitted). 

(vi)  “A  Nuclear  Magnetic  Resonance  Study  of  Dynamics  in  Toluene- 
Polyisobutylene  Solutions:  (III)  Segmental  Motion”  A.  Bandis,  P.T. 
Inglefield,  A.A.  Jones  and  W.-Y.  Wen.  J.  Polymer  Science.  Poly.  Phys.  Ed. 
(submitted). 

(vii)  “A  NMR  Study  of  Penetrant  Diffusion  and  Polymer  Segmental 
Dynamics  in  Rubber”  A.  Bandis,  P.T.  Inglefield,  A.A.  Jones  and  W.-Y.  Wen. 
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Polymer  Materials  Sci.  and  Eng.  ACS.  71,  217  (1994) 

(vili)“A  NMR  Study  of  Penetrant  Diffusion  and  Polymer  Segmental  Motion 
in  Toluen-Polyisobutylene  Solutions”  A.  Bandis,  P.T.  Inglefield,  A.A.  Jones 
and  W.-Y.  Wen.  ACS  Symposium  Series,  A.  Marek,  Editor  (submitted). 


Manuscripts  to  follow: 


Nuclear  Spin  Relaxation  Dynamics  of  ^^C02  Sorbed  in 
Polyisobutene  Rubber 


Z.  P.  DONG,  B,  |.  CAULEY,  A.  BANDIS,  C.  W.  MOU,  C  E.  INGLEFIELD, 
A.  A.  JONES,  P.  T.  INGLEFIELD,  and  W.-Y.  WEN 

Department  of  Chemistry,  Clark  University,  Worcester,  Massachusetts  01610 


SYNOPSIS 

Recently  we  presented  the  d3mamics  of  ^®C02  molecules  sorbed  in  silicone  rubber  (PDMS) 
ascertained  from  spin  relaxation  experiments.  Results  of  a  similar  investigation  for  ^^C02 
sorbed  in  polyisobutene  (PIB)  are  presented  in  this  report.  The  spin-lattice  and  spin-spin 
relaxation  times  as  well  as  nuclear  Overhauser  enhancements  (NOE)  were  determined  as 
a  function  of  temperature  and  Larmor  frequency.  The  relaxation  mechanisms  found  to  be 
important  for  ^^C02/PIB  system  are  intermolecular  dipole-dipole  relaxation  and  chemical 
shift  anisotropy  with  a  minor  contribution  from  spin  rotation  relaxation.  We  have  deter¬ 
mined  the  parameters  which  characterize  correlation  times  for  ^^C02  collisional  motion, 
rotational  motion,  and  translational  motions  in  the  PIB.  The  self-diffusion  coefficient  of 
5.15  X  10"®  cm^/s  obtained  from  the  nuclear  magnetic  resonance  (NMR)  data  is  close  to 
the  literature  value  of  the  mutual  diffusion  coefficient  of  CO2  in  PIB  at  300  K  obtained 
from  permeability  measurements.  In  contrast  to  the  case  of  CO2/PDMS  in  which  a  broad 
distribution  (characterized  by  a  fractional  exponential  correlation  function  of  the  Williams- 
Watts  type  with  a  =  0.58)  is  observed,  a  sharp  distribution  with  a  fractional  exponent,  a, 
of  0.99  is  found  for  the  CO2/PIB  system.  Instead  of  assuming  an  Arrhenius  type  temperature 
dependence,  we  used  a  Williams -Landel-Ferry  type  temperature  dependence  and  found  it 
to  be  better  suited  to  describe  the  behavior  of  this  system.  PIB  is  a  densely  packed  “strong” 
chain  polymer  which  responds  gradually  to  the  temperature  variation  and  gas  sorption.  In 
contrast  PDMS  is  a  relatively  loosely  packed  “fragile”  polymer  with  a  propensity  to  exhibit 
rapid  dynamic  responses  to  the  temperature  change  and  gas  sorption.  ©  1993  John  Wiley  & 
Sons,  Inc. 

Keywords:  carbon  dioxide  •  spin  relaxation  •  polyisobutylene 


INTRODUCTION 

Recently  we  presented  a  study  of  the  nuclear  spin 
relaxation  of  sorbed  ^®C02  in  silicone  rubber  (poly¬ 
dimethyl  siloxane ) }  Our  analysis  showed  that  the 
single-mode  hypothesis  can  be  used  to  interpret  the 
experimental  spin-lattice  relaxation  times.  However, 
this  hypothesis  with  a  single  correlation  time  was 
found  to  be  incapable  of  predicting  the  nuclear 
Overhauser  enhancements  and  spin-spin  relaxation 
times.  On  the  other  hand,  if  a  distribution  of  cor¬ 
relation  times  was  assumed,  we  were  able  to  interpret 
all  of  the  experimental  data  in  a  reasonable  manner. 


Journal  of  Polymer  Science:  PartB:  Polymer  Physics,  Vol.  31, 1213-1221  (1993) 
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We  have  quantitatively  analyzed  the  spin  relaxation 
mechanisms  and  determined  the  parameters  which 
characterize  the  correlation  times  for  ^^C02  colli¬ 
sional  motions,  rotational  motions,  and  translational 
motions  in  the  silicone  rubber.  A  fractional  expo¬ 
nential  correlation  function  of  the  Williams- Watts 
type  with  a  fractional  exponent,  a,  of  0.58  was  found 
in  the  description  of  the  ^^C02  dynamics.  This  ex¬ 
ponent  is  close  to  the  distribution  of  time  scales  as¬ 
sociated  with  segmental  motion  in  the  silicone  rub¬ 
ber  with  a  value  of  a  in  the  range  of  0.47-0.51.^ 

We  have  extended  our  nuclear  magnetic  reso¬ 
nance  (NMR)  investigations  to  another  rubbery 
polymer,  polyisobutene  ( PIB )  and  determined  the 
spin-lattice  relaxation  time,  Ti,  spin-spin  relaxation 
time,  T2,  and  nuclear  Overhauser  enhancement 
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(NOE)  of  ^^C02  sorbed  in  the  polymer  as  a  function 
of  temperature  and  Larmor  frequency. 

Polyisobutene,  an  elastomer  with  properties  and 
behavior  considerably  different  from  those  of  sili¬ 
cone  rubber,  is  chosen  for  our  investigations.  PIB*s 
wide  working  range  makes  it  a  convenient  elastomer 
for  application.  Its  glass  transition  temperature  ( Tg 
=  205  K)  is  suitably  below  the  temperature  at  which 
the  elasticity  is  desired.  It  has  good  chain  mobility 
but  not  so  good  that  the  chains  can  readily  form 
crystals.^  PIB  is  classified  as  a  “strong’^  chain  poly¬ 
mer,  while  PDMS  is  considered  as  one  of  the  “frag¬ 
ile”  polymers  according  to  Angell  et  al.®  In  addition 
PIB  is  an  elastomer  notable  for  low  rates  of  diffusion 
of  small  penetrants."^  The  low  diffusion  rates  have 
been  attributed  to  efficient  intermolecular  packing 
and  low  total  fractional  free  volume  in  PIB.®  Some 
interesting  results  we  have  obtained  are  reported 
below.  For  example,  the  rotational  correlation  time 
at  300  K  for  CO2  in  PIB  is  about  800  times  longer 
than  that  in  PDMS.  The  fractional  exponent  a  is 
0.99  for  the  ^®C02/PIB  system  in  contrast  to  the 
corresponding  value  of  0.58  for  the  ^^C02/PDMS 
system. 


EXPERIMENTAL 

The  polyisobutene  sample  used  for  the  study  was 
obtained  from  Cellomer  Associates,  Inc.  of  Webster, 
NY  in  the  form  of  thick  sheets  of  high-molecular 
weight  (Cat.  #40E,  Lot  #02).  The  sheets  were  cut 
into  cylinders  7  mm  diameter  and  11  mm  in  length. 
Four  of  these  cylinders  (44  mm  in  height)  were 
placed  in  a  heavy-wall  NMR  tube  ( 10  mm  o.d.  and 
1.45  mm  wall  thickness) ,  then  degassed  under  vac¬ 
uum  for  several  days.  Measured  amounts  of  ^^C02 
(99.8%,  Merck  &  Co.)  were  introduced  into  the 
NMR  tube  and  the  CO2  pressure  was  monitored  with 
a  transducer-electrometer  system.  The  NMR  tube 
was  cooled  to  the  liquid-nitrogen  temperature  and 
then  sealed.  The  volume  of  the  PIB  sample  was  1.69 
cm^  and  the  volume  of  the  tube  above  the  sample 
was  3.96  cm®.  The  total  amount  of  carbon  dioxide 
in  the  tube  was  2.16  X  10“®  mol  and  it  exerted  a 
pressure  of  about  10  atm  at  room  temperature. 
( Concentration  of  the  dissolved  ^®C02  in  PIB  was 
about  1.5%  by  weight  at  room  temperature.) 

Carbon-13  spin-lattice  and  spin-spin  relaxation 
times  (Ti  and  T2)  were  measured  at  Larmor  fre¬ 
quencies  of  22.6,  50.3,  and  75.4  MHz  on  Bruker  SXP 
20-100,  AC-200,  and  MSL-300  spectrometers,  re¬ 
spectively.  A  simple  180°  r  90°  pulse  sequence  was 
used  to  measure  Ti  while  90°  r  180°  sequence  was 


used  to  measure  T^,  Gated  NOE  experiments  were 
performed  in  the  AC-200  spectrometer  while  con¬ 
tinuous  decoupling  was  employed  on  the  SXP  spec¬ 
trometer  with  a  decoupler  power  of  6  or  12  W. 


RESULTS 

Values  of  the  spin-lattice  relaxation  time,  Ti,  for 
^®C02  sorbed  in  PIB  at  various  temperatures  and  a 
Larmor  frequency  of  50.3  MHz  are  presented  in  Fig¬ 
ure  1  as  an  example.  The  plot  of  Ti  versus  temper¬ 
ature  shows  a  broad  minimum  with  a  Ti  value  of 
2.9  s  at  around  290  K.  It  turns  out  that  the  values 
of  Ti  alters  only  slightly  with  the  change  of  Larmor 
frequency  in  the  range  22-75  MHz. 

The  values  of  spin-spin  relaxation  time,  T2,  for 
^®C02  sorbed  in  PIB  at  two  Larmor  frequencies  are 
shown  in  Figure  2.  In  contrast  to  Ti,  there  is  no 
minimum  and  the  values  of  T2  increase  rapidly  with 
the  increase  of  temperature.  The  values  of  nuclear 
Overhauser  enhancement  (NOE)  for  ^®C02  in  PIB 
are  plotted  as  a  function  of  temperature  at  Larmor 
frequencies  of  22.6  and  50.3  MHz  in  Figure  3.  The 
NOE  values  are  somewhat  scattered  but  they  seem 
to  increase  slightly  with  the  temperature  and  lie  in 
the  range  1.06-1.43.  Experimental  Ti,  72,  and  NOE 
data  are  summarized  and  listed  in  Tables  I,  II, 
and  III. 


Figure  1.  Spin-lattice  relaxation  time  Ti  at  50.3  MHz 
for  ^^C02  (10  atm)  sorbed  in  polyisobutene  as  a  function 
of  temperature.  The  line  is  a  fit  based  on  the  Brownian 
diffusion  model  for  temperatures  above  Tg  ( ^  285  K ) . 


SPIN  RELAXATION  OF  SORBED  '^COj  1215 


Figure  2 .  Temperature  dependence  of  T2  at  two  Larmor 
frequencies,  50.3  MHz  (♦)  and  75.4  MHz  (O)  for  ^^C02 
( 10  atm)  sorbed  in  polyisobutene.  The  lines  are  fits  based 
on  the  Brownian  diffusion  model  for  temperatures  above 
T^(285  -  290  K). 

Interpretation  and  Analysis 

Glass  Transition  Temperature 

The  glass  transition  temperature  of  PIB  as  measured 
by  differential  scanning  calorimetry  (DSC)  or  dy- 


Figure  3.  Temperature  dependence  of  NOE  at  two 
Larmor  frequencies,  22.6  MHz  (  + )  and  50.3  MHz  ( ♦ )  for 
^^C02  ( 10  atm)  sorbed  in  polyisobutene.  The  lines  are  fits 
based  on  the  Brownian  diffusion  model  for  temperatures 
above  Tg  (280  285  K). 


Table  I.  Experimental  Ti  and  NOE  Data  for 
(10  atm)  Sorbed  in  Polyisobutene  at  22.6  MHz 

‘'CO2 

Temp.  (K) 

r,  (8) 

NOE 

233 

_ 

i.06r 

253 

3.39“ 

1.17“ 

273 

3.60“ 

1.21 

293 

3.04 

1.20 

313 

2.93 

1.43 

333 

4.02 

— 

353 

5.16 

— 

373 

6.09 

— 

•  These  data  are  those  at  temperatures  close  to  or  below  the 
glass  transition  temperature  (280  K)  of  the  CO2/PIB  system  and, 
therefore,  not  used  in  our  model  fit  to  the  theory  which  is  appli’ 
cable  to  rubbery  polymers. 

Note:  r,  =  205  K  +  85  K  -  10  K  =  280  K. 

namic  mechanical  techniques  is  205  K.  These  meth¬ 
ods  are  sensitive  to  low  frequency  response  (ca.  1 
Hz)  while  the  frequencies  which  are  characterized 
by  Ti  are  the  Larmor  frequencies  of  carbon-13  nuclei 
and  protons,  and  both  of  these  lie  in  the  mega-Hertz 
range.  We  can  use  the  WLF  equation  to  estimate 
the  temperature  shift  (AT^)  associated  with  the 
higher  frequency  of  the  NMR  measurements. 
Mandelkern®  has  used  this  approach  to  characterize 
the  relationship  between  the  glass  transition  and 
carbon-13  relaxation  for  a  number  of  polymers.  The 
actual  form  of  the  WLF  equation  to  be  employed  is 

log(w2/wi)  =  Ci(ATg)/{c2  +  (ATg)) 
where  Ci  =  16.8  and  C2  =  108.6  for  PIB.  For  coi 


Table  II.  Experimental  Ti,  T2,  and  NOE  Data  for 
^®C02  (10  atm)  Sorbed  in  Polyisobutene  at  50.3  MHz 


Temp.  (K)  Ti  (s)  T2  (s)  NOE 


253 

4.2“ 

— 

273 

3.49“ 

0.070“ 

1.27' 

293 

2.89“ 

0.440“ 

1.14' 

303 

3.13 

— 

1.12 

313 

3.10 

1.27 

1.16 

323 

3.38 

1.55 

— 

333 

3.70 

2.21 

1.29 

343 

4.80 

3,34 

— 

353 

5.20 

3.80 

1.33 

373 

8.04 

4.82 

— 

•  These  data  are  shown  at  temperatures  close  to  or  below  the 
glass  transition  temperature  (285  K)  of  the  CO2/PIB  system  and, 
therefore,  not  used  in  our  model  fit  to  the  theory  which  is  appli¬ 
cable  to  rubbery  polymers. 

Note:  7i  -  205  K  +  90  K  -  10  K  =  285  K. 
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Table  DI.  Experimental  Ti  and  Tg  Data  for  ^^C02 
(10  atm)  Sorbed  in  Polyisobutene  at  75.4  MHz 


Temp.  (K) 

T,(s) 

r,(s) 

260 

4.62* 

6.92  X  10'^ 

280 

— 

0.0667* 

300 

1.94* 

0.635* 

311 

2.93 

0.769 

322 

3.14 

1.14 

331 

3.37 

1.30 

340 

3.96 

2.13 

350 

4.87 

3.06 

360 

5.64 

4.14 

•  These  data  are  those  at  temperatures  close  to  or  below  the 
glass  transition  temperature  (290  K)  of  the  CO2/PIB  system  and, 
therefore,  not  used  in  our  model  fit  to  the  theory  applicable  to 
rubbery  polymers. 

Note:  7;  =  205  K  +  95  K  -  10  K  =  290  K. 


of  1  Hz  and  CO2  in  the  range  22-75  MHz,  the  values 
of  LTg  fall  in  the  range  85-95  K.  This  will  bring  the 
spin-lattice  relaxation  estimate  of  Tg  up  into  the 
range  of  290-300  K. 

On  the  other  hand,  the  glass  transition  temper¬ 
ature  is  expected  to  be  lowered  by  the  sorption  of 
CO2  ( 10  atm)  and  the  subsequent  plasticization  of 
PIB.  The  depression  of  Tg  by  this  mechanism  is  es¬ 
timated  to  be  about  10  K.  A  combination  of  the  ef¬ 
fects  of  CO2  plasticization  and  the  high  frequency 
measurements  would  bring  the  Tg  into  the  range 
280-290  K  for  the  ^^C02/PIB  system.  This  range  of 
the  glass  transition  temperature  seems  to  be  con¬ 
sistent  with  our  observations  that  (a)  the  T\  data 
show  a  change  in  slope  in  the  range  290-310  K,  (b) 
the  T2  value  decreases  suddenly  from  about  0.5  s  to 
less  than  0.1  s  when  the  temperature  is  lowered  from 
290  K  to  280  K,  and  (c)  as  the  temperature  is  raised 
from  273  to  293  K,  the  dominant  relaxation  mech¬ 
anism  of  CO2  changes  from  CSA  to  dipole-dipole  at 
a  Larmor  frequency  of  22.6  MHz  (see  Discussion) . 

Relaxation  Mechanisms 

Let  us  focus  our  attention  on  the  sources  of  relax¬ 
ation  for  ^^C02  above  the  WLF  shifted  glass  tran¬ 
sition  in  PIB  and  analyze  the  experimental  data  ob¬ 
tained  at  or  above  293  K.  Three  mechanisms  are 
contributing  to  the  NMR  relaxation  of  ^^C02 1  spin 
rotation  (SR)  relaxation,  chemical  shift  anisotropy 
(CSA)  relaxation,  and  intermolecular  dipole-dipole 
(DD)  relaxation.  The  same  equations^  as  used  pre¬ 
viously  to  analyze  the  relaxation  of  ^^C02  sorbed  in 
silicone  rubber  are  applied  to  study  the  relaxation 
for  ^^C02  in  PIB.  In  addition  to  assuming  the  Ar¬ 


rhenius  type  temperature  dependence,  we  also  em¬ 
ploy  the  WLF  type  temperature  dependence  for  cor¬ 
relation  times. 

log[T(T)/T(To)] 

=  c?(T~To)/(c?  +  T-To)  (1)^ 

where  =  298  K,  c?  =  8.61,  and  c§  =  200.4  K.® 
Instead  of  fixing  c®  as  a  constant,  we  treat  c®  as  an 
adjustable  parameter  and  let  it  play  a  role  similar 
to  AH,  the  energy  of  activation  in  the  Arrhenius 
type  equation. 

In  the  previous  study  of  ^^C02  in  silicone  rubber, 
it  was  necessary  to  assume  the  distribution  of  dy¬ 
namical  CO2  motions.^  We  shall  also  use  the  Wil- 
liams-Watts  type  correlation  function  to  depict  mo¬ 
tion  of  CO2  in  PIB: 

0(t)  =  exp[-(i/Tp)“]  (2) 

In  eq.  ( 2 ) ,  Tp  is  the  central  correlation  time  and  sets 
the  position  of  the  distribution  on  the  time  axis. 
The  parameter  a,  which  falls  in  the  range  of  0  <  a 
1,  determines  the  width  of  the  distribution. 

There  are  a  total  of  37  data  points  and  our  model 
contains  at  least  six  potentially  adjustable  pa¬ 
rameters.  They  are  correlation  times  tsr(300  K), 
tcsa  (300  K),  and  tdd  (300  K)  for  spin  rotation, 
chemical  shift  anisotropy,  and  dipole-dipole  inter¬ 
action,  respectively;  width  parameter  for  distribution 
(a),  distance  of  closest  approach  for  CO2  to  PIB 
protons  (b),  and  either  the  WLF  parameter  (c®)  or 
the  Arrhenius  activation  energy  (AH), 

The  fitting  procedure  now  comes  into  play  and 
the  values  of  Ti,  T2,  and  NOE  computed  by  the 
equations^  are  compared  with  the  experimental  val¬ 
ues.  Then  the  values  of  the  six  adjustable  parameters 
are  changed  in  an  effort  to  minimize  the  sum  of 
squares  (the  least  squares). 

The  Best  Fit  Results 

The  resulting  fit  has  root-mean-square  residuals  of 
11%,  which  is  also  the  approximate  value  of  our  ex¬ 
perimental  error.  The  results  of  the  fits  for  the  T 1 , 
T2,  and  NOE  data  are  shown  in  Figures  1,  2,  and  3, 
respectively.  Note  that  in  all  cases  the  fits  are  for 
the  rubbery  regime  above  Tg  only.  The  parameters 
of  the  best  fit  are: 

rsR(300  K)  =  2.9  X  10’'^^  s 
tcsa(300K)  =  1.3  X  lO'^^^s 
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rDD(300K)  =  2.6X10-®s 
a  =  0.99 

cl  =  8.6;  AH  =  28.2  kJ/mol 
b  =  3.66  X  10  cm 

Since  the  distance  of  the  closest  approach  for  car- 
bon-13  of  ^^C02  to  PIB  protons  can  be  estimated 
from  the  van  der  Waals  radii  and  the  result  is  close 
to  the  above  value,  b  can  be  eliminated  and  no  longer 
considered  as  an  adjustable  parameter.  The  WLF 
parameter  cl  determined  by  the  fitting  process  turns 
out  to  be  8.6,  which  is  practically  identical  to  the 
literature  value®  and  therefore  it  can  be  eliminated 
also.  The  remaining  parameters  are  just  four:  three 
r’s  and  a. 


DISCUSSION 

The  overall  relaxation  rate  at  each  temperature  and 
Larmor  frequency  is  a  sum  of  the  rates  from  three 
mechanisms: 

1/Ti  =  l/ri(CSA)  +  l/Ti(DD) 

+  l/Ti(SR)  (3) 

In  order  to  see  the  contributions  of  each  mechanism 
to  the  overall  relaxation,  calculated  values  of  the 


Figure  4,  The  chemical  shift  anisotropy  contribution 
to  the  spin-lattice  relaxation  rate  at  three  Larmor  fre¬ 
quencies,  22.6  MHz  ( - ),  50.3  MHz  (-•-),  and  75.4 

MHz  ( - )  as  a  function  of  temperature. 


Figure  6.  The  dipole-dipole  interaction  contribution  to 
the  spin-lattice  relaxation  rate  at  three  Larmor  frequen¬ 
cies,  22.6  MHz  ( - ),  50.3  MHz  (-•-),  and  75.4  MHz 

( - )  as  a  function  of  temperature. 

rates,  l/Ti(CSA),  l/Ti(DD),  and  l/Ti(SR)  are 
plotted  in  Figures  4,  5,  and  6,  respectively,  as  func¬ 
tions  of  temperature  at  three  Larmor  frequencies. 
Though  our  calculations  are  based  on  the  experi¬ 
mental  data  for  ^®C02  in  rubbery  PIB  in  the  tem¬ 
perature  range  303-373  K,  we  have  expanded  the 
range  of  our  plots  to  cover  a  wider  temperature  range 


Figure  6.  The  spin  rotation  contribution  to  the  spin- 
lattice  relaxation  rate  as  a  function  of  temperature. 
T*i(SR)  is  independent  of  Larmor  frequency. 
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Table  IV.  Contributions  of  the  Chemical  Shift  Anisotropy  Mechanism  to  the  Overall  Spin-Lattice  Relaxation 
of  ^^C02  Sorbed  in  Rubbery  Polyisobutene 


Temp  (K) 

22.6  MHz 

50.3  MHz 

75.4  MHz  •  ‘ 

Ti  (CSA)  (s) 

Percent 
of  Total 

Ti  (CSA)  (s) 

.  Percent 
of  Total 

Ti  (CSA)  (s) 

Percent 
of  Total 

233 

1.31 

98 

0.99 

99.7 

0.93 

99.9 

253 

1.52 

91 

0.50 

99 

0.33 

99.8 

273 

4.01 

61 

1.00 

96 

0.51 

99 

293 

10.7 

27 

2.32 

83 

1.10 

95 

313 

25.9 

12 

5.32 

61 

2.40 

85 

333 

54.9 

6.8 

11.2 

41 

4.96 

70 

353 

104 

4.4 

21.1 

27 

9.36 

54 

373 

181 

3.1 

36.7 

19 

16.2 

40 

393 

293 

2.2 

59.4 

13 

26.3 

29 

233-393  K  for  the  purpose  of  seeing  some  features 
of  interest  which  may  be  taking  place  at  or  slightly 
below  the  spin-lattice  relaxation  estimate  of  Tg  ( 280- 
290  K).  The  estimated  percent  contributions  are 
listed  in  Tables  IV,  V,  and  VI,  and  also  plotted  in 
Figure  7  as  a  function  of  temperature  at  22.6  MHz. 

As  can  be  seen  from  Figure  7,  the  CSA  and  DD 
contributions  to  the  spin-lattice  relaxation  are  sub¬ 
stantial  and  change  greatly  as  a  function  of  tem¬ 
perature  in  the  range  250-330  K.  As  the  temperature 
increases,  the  dominant  mechanism  alters  from  CSA 
to  DD,  and  the  two  curves  (%CSA  vs  T  and  %DD 
vs  T)  cross  each  other  at  around  280  K,  which  is 
the  glass  transition  temperature  for  the  ^^062/ PIB 
system  at  the  Larmor  frequency  of  22.6  MHz.  For 
CO2  sorbed  in  glassy  polycarbonate  and  glassy  poly- 
st3Tene,  the  dominant  relaxation  mechanism  is  CSA. 


The  drop  in  the  contribution  of  this  mechanism  at 
around  280  K  in  PIB  is  further  evidence  of  the 
change  from  glassy  behavior  to  rubbery  behavior  at 
the  frequency  of  the  NMR  experiment. 

For  comparison  with  CO2  in  another  rubbery 
polymer,  we  show  the  percent  contributions  of  each 
mechanism  to  the  spin -lattice  relaxation  for  ^®C02/ 
PDMS  system^  at  22.6  MHz  in  Figure  8.  In  the 
PDMS  case,  the  CSA  mechanism  is  a  very  minor 
contribution,  which  is  consistent  with  a  rubbery 
system.  The  dominant  mechanisms  for  CO2  in 
PDMS  are  SR  and  DD,  and  their  contributions 
change  only  slightly  with  the  temperature.  It  should 
be  pointed  out  that  the  temperature  range  of  our 
interest,  230-390  K,  is  considerably  above  the  spin- 
lattice  relaxation  time  view  of  Tg  ( 192  K)  for  ^^C02/ 
PDMS  system  at  22.6  MHz.  The  weak  temperature 


Table  V.  Contributions  of  the  Dipole-Dipole  Mechanism  to  the  Overall  Spin-Lattice  Relaxation 
of  ^®C02  Sorbed  in  Rubbery  Polyisobutene 

Temp  (K) 

22.6  MHz 

50.3  MHz 

75.4  MHz 

T,  (DD)  (s) 

Percent 
of  Total 

Ti  (DD)  (s) 

Percent 
of  Total 

T,  (DD)  (s) 

Percent 
of  Total 

233 

74.9 

1.7 

354 

0.3 

786 

0.09 

253 

15.9 

8.7 

69.1 

0.7 

149 

0.2 

273 

6.16 

39 

22.8 

4.2 

46.4 

1 

293 

4.00 

73 

11.7 

16 

21.8 

5 

313 

3.70 

87 

8.70 

38 

14.5 

14 

333 

4.13 

91 

8.20 

56 

12.4 

28 

353 

5.06 

90 

8.87 

66 

12.5 

40 

373 

6.49 

86 

10.3 

67 

13.8 

47 

393 

8.47 

77 

12.4 

62 

15.9 

47 
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Table  VI.  Contributions  of  the  Spin  Rotation  Mechanism  to  the  Overall  Spin-Lattice  Relaxation 
of  ^^C02  Sorbed  in  Rubbery  Polyisobutene 


%  Contribution  of  the  Spin  Rotation  Mechanism 


Temp  (K) 

(SR)  (s)« 

22.6  MHz 

50.3  MHz 

75.4  MHz 

233 

60040 

0 

0 

0 

253 

10260 

0 

0 

0 

273 

2602 

0.1 

0 

0 

293 

868 

0.3 

0.2 

0.1 

313 

352 

0.9 

1.0 

0.6 

333 

166 

2.3 

2.8 

2 

353 

87.3 

5.2 

6.4 

6 

373 

50.2 

11.1 

14 

13 

393 

31.0 

21.0 

25 

24 

•  Calculated  spin-lattice  relaxation  times  due  to  the  spin-rotation  mechanism. 
Ti  (SR)  is  independent  of  the  Larmor  frequency. 


dependence  is  consistent  with  the  “WLF  behavior” 
at  temperatures  well  above  Tg.  The  importance  of 
the  SR  mechanism  in  PDMS  relative  to  PIB  is  also 
an  indication  of  greater  mobility  in  the  lower  Tg  ma¬ 
terial. 

A  distribution  of  CO2  motions  in  the  silicone  rub¬ 
ber  was  observed  and  characterized  by  a  fractional 
exponential  correlation  function  of  the  Williams- 
Watts  type  with  a  fractional  exponent,  a,  of  0.58.^ 
In  contrast,  the  distribution  of  CO2  motions  in  PIB 


Figure  7.  Percent  contributions  of  CSA,  DD,  and  SR 
mechanisms  at  22.6  MHz  for  “CO2  (10  atm)  sorbed  in 
polyisobutene  as  a  function  of  temperature. 


must  be  very  narrow  since  the  value  of  the  width 
parameter  a  is  0.99.  Consequently  the  population 
distribution  of  the  self-diffusion  coefficient,  D,  for 
CO2  sorbed  in  PIB  will  also  be  very  narrow.  The 
mean  value  of  D  at  300  K  is  given  by 

<D>  =  fcV<TDD>  =  5.15  X  10'®  cmVs  (4) 

which  is  close  to  the  mutual  diffusion  coefficient  of 
5.78  X  10“®  cm^/s  reported  by  Van  Amerongen^ 


Figure  8.  Percent  contributions  of  CSA,  DD,  and  SR 
mechanisms  at  22.6  MHz  for  ^^C02  ( 10  atm)  sorbed  in 
silicone  rubber  as  a  function  of  temperature. 
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from  his  permeability  measurements  of  CO2  sorbed 
in  PIB.  Given  the  complexity  of  the  interpretation 
of  the  NMR  data  and  the  short-range  viewpoint  of 
the  NMR  experiment,  the  agreement  is  better  than 
should  be  expected. 

Good  agreement  between  the  NMR  results  and 
permeability  results  was  also  obtained  from  the  in¬ 
terpretation  of  the  PDMS  system.  This  lends  some 
credibility  to  the  connections  between  the  dynamics 
determined  from  NMR  and  from  permeability  data. 

With  regard  to  the  spin  rotation,  the  correlation 
time  TsR  for  CO2  in  PIB  is  2.9  X  10”^®  s  which  is 
smaller  by  a  factor  of  30  than  the  corresponding 
value  of  8.8  X  10”^^  s  in  PDMS.  This  is  explicable 
in  terms  of  the  number  of  collisions  suffered  by  CO2 
molecules  in  these  polymers:  it  is  much  greater  in 
PIB  than  in  PDMS,  since  the  former  is  a  much  more 
tightly  packed  polymer  than  the  latter.®  Next,  the 
value  of  rcsA  is  1.3  X  10’^°  s  in  PIB  which  may 
be  compared  with  the  corresponding  value  of  1.6 
X  10“^®  s  at  300  K  in  PDMS.  The  fact  that  the  ro¬ 
tational  correlation  time  is  almost  three  orders  of 
magnitude  greater  in  PIB  than  in  PDMS  can  also 
be  attributed  to  the  tighter  packing  of  the  former 
than  the  latter.® 

If  an  Arrhenius  fit  is  used,  the  energy  of  activation 
AH  for  CO2  in  PIB  is  found  to  be  28  kJ/mol  which 
is  considerably  greater  than  the  corresponding  value 
of  only  1.6  kJ/mol  in  PDMS.  Qualitatively  this  dif¬ 
ference  is  in  the  direction  to  be  expected,  but  quan¬ 
titatively  the  difference  is  extremely  large.  It  seems 
to  require  re-examination  of  the  meaning  of  the  ac¬ 
tivation  energy  for  diffusion  of  sorbed  CO2  in  poly¬ 
mers.  Since  a  rubbery  polymer  is  known  to  follow 
the  WLF  equation®  and  not  the  Arrhenius  type 
equation,  the  value  of  AH  estimated  is  expected  to 
change  with  the  temperature  of  observation.  In  other 
words,  AH  is  a  function  of  temperature  T,  but  more 
correctly,  it  is  a  function  of  T I  Tg  where  Tg  is  the 
glass  transition  temperature  of  CO2 /polymer  sys¬ 
tem.  If  the  temperature  T  is  300  K  and  the  Larmor 
frequency  is  22.6  MHz,  then  the  values  of  T/Tg  are 
1.07  and  1.56  for  CO2/PIB  and  CO2/PDMS,  re¬ 
spectively.  At  300  K,  the  CO2/PIB  system  is,  there¬ 
fore,  much  closer  to  Tg  than  the  CO2/PDMS  system. 
The  value  of  AH  near  glass  transition  (CO2/PIB 
case)  would  be  much  greater  than  that  at  temper¬ 
ature  far  above  the  glass  transition  (CO2/PDMS 
case )  in  agreement  with  our  observation. 

In  the  fitting  process,  the  best  result  is  obtained 
with  the  value  of  a.  nearly  equal  to  1,  though  the 
quality  of  the  fit  changes  little  as  a  is  varied  from 
0.8  to  1.0.  For  this  very  narrow  distribution,  eq.  (4) 


is  accurate  but  it  does  not  apply  to  the  PDMS  system 
where  a  broad  distribution  is  in  effect.  For  silicone 
rubber  at  300  K  we  have 

<tdd>  =  7.8X10"'^s  and 

(D>  =  1.8  X  10"®  cmVs.  (5) 

Note  that  (tdd)  for  PIB  is  a  factor  of  30  smaller 
but  ^  ^  is  a  factor  of  350  smaller.  It  is  the  effect  of 
the  istribution  of  correlation  times  in  PDMS  which 
leads  to  the  relatively  high  permeability  in  that  sys¬ 
tem. 

The  time  scale  of  segmental  motion  in  PDMS  at 
300  K  is  2  X  10"^®  s  based  on  an  interpretation  of 
proton  spin  relaxation  data  using  a  Williams-Watts 
stretched  exponential  correlation  function  with  an 
exponent  of  0.47.^®  A  preliminary  measurement  of 
carbon-13  Ti  data  on  PIB  shows  a  Ti  minimum  at 
slightly  above  room  temperature  and  at  a  Larmor 
frequency  of  75  MHz.  This  corresponds  to  a  corre¬ 
lation  time  of  about  2  X  10"®  s.  Thus  for  PIB,  tdd 
for  CO2  is  within  an  order  of  magnitude  of  the  cor¬ 
relation  time  for  segmental  motion.  This  approxi¬ 
mate  coincidence  would  indicate  a  possible  relation¬ 
ship;  namely,  the  backbone  motion  of  the  PIB  gates 
the  translational  diffusion  of  the  CO2.  In  PDMS, 
the  Tdd  for  CO2  is  6  orders  of  magnitude  larger  than 
the  correlation  time  for  segmental  motions  indicat¬ 
ing  a  substantially  different  relationship  between 
segmental  motion  and  translational  diffusion  of  the 
sorbed  CO2.  A  more  definitive  correlation  of  pol3nner 
motions  and  penetrant  motions  will  be  left  for  future 
considerations.  The  position  of  the  T 1  minimum  in 
PIB  is  very  close  to  the  shift  in  temperature  esti¬ 
mated  from  the  WLF  equation.  It  should  also  be 
noted  that  the  temperature  at  which  carbon-13 
spectra  can  no  longer  be  resolved  under  conditions 
of  simple  scalar  decoupling  is  only  also  slightly  lower 
than  the  temperature  of  the  Ti  minimum.  The 
proximity  of  the  minimum  and  the  temperature  at 
which  scalar  decoupling  fails  to  produce  resolution 
is  especially  close  in  PIB  and  reflects  the  unusual 
WLF  parameter  associated  with  this  polymer.® 

This  behavior  of  the  two  contrasting  polymer 
systems,  PIB  and  PDMS,  may  be  considered  in  the 
context  of  the  description  developed  by  Angell  et 
al.®  In  the  CO2/PDMS  system,  CO2  molecules  probe 
the  dynamics  of  the  “fragile”  polymer  with  a  broad 
distribution.  In  the  CO2/PIB  system,  on  the  other 
hand,  CO2  molecules  probe  the  “strong”  polymer 
with  a  narrow  distribution  of  correlation  times.  PIB 
is  a  densely  packed  polymer  with  a  small  total  free 
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volume,®  and  responds  gradually  to  the  temperature 
variation  and  gas  sorption.  In  contrast,  PDMS  is  a 
relatively  loosely  packed  polymer  with  a  propensity 
to  exhibit  rapid  dynamic  response  to  gas  sorption 
and  temperature  change. 
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SYNOPSIS 

The  spin-lattice  relaxation  times  as  well  as  nuclear  Overhauser  enhancements  (NOE)  are 
determined  for  the  methylene  carbon  of  polyisobutylene  ( PIB ) .  Three  different  correlation 
functions  are  used  to  treat  the  relaxation  data  to  describe  the  local  motion  of  the  PIB 
backbone  rubber.  The  first  is  the  Hall-Helfand  function  combined  with  restricted  anisotropic 
rotational  diffusion,  the  second  is  the  Dejean-Laupretre-Monnerie  function,  and  the  third 
is  the  Williams-Watts  function  combined  with  restricted  anisotropic  rotational  diffusion. 

Spin-lattice  relaxation  times  are  also  determined  on  PIB  under  pressure  of  CO2  gas  to 
check  for  changes  in  segmental  motion  in  the  presence  of  this  penetrant.  All  the  models 
give  a  satisfactory  description  for  the  local  dynamics  of  the  backbone  carbon  of  the  polymer 
chain,  yielding  comparable  values  for  the  time  scale  of  conformational  changes,  apparent 
activation  energies  and  distribution  of  correlation  times.  In  addition,  a  proton  line  shape 
experiment  is  performed  on  bulk  PIB  as  a  function  of  temperature.  The  results  are  inter¬ 
preted  in  terms  of  two  motional  components  in  the  polymer  chains:  fast  and  slow  ones. 

Applying  a  revised  version  of  Rossler^s  formalism  to  our  data,  an  average  energy  barrier  is 
found  which  is  consistent  with  the  apparent  activation  energies.  The  apparent  activation 
energy  obtained  for  PIB,  from  the  nuclear  magnetic  resonance  (NMR)  measurements, 
about  35  kJ,  is  slightly  higher  than  the  apparent  activation  energy  for  the  translational 
diffusion  of  CO2  in  PIB  from  NMR  data.  The  distribution  of  correlation  times  and  the 
correlation  times  themselves  for  segmental  motion  of  the  PIB  and  translational  diffusion 
of  CO2  in  PIB  are  very  similar,  indicating  a  close  link  between  the  two  motional  processes. 

The  distribution  of  correlation  times  is  close  to  a  single  exponential  in  PIB  consistent  with 
Angell  characterization  of  this  material  as  a  strong  liquid.  ©  1994  John  Wiley  &  Sons,  Inc, 

Keywords:  polyisobutylene  •  chain  dynamics  •  gas  diffusion 

above  the  glass  transition.  In  the  Angell  classifica¬ 
tion,^  PIB  is  a  strong  liquid.  Thus  the  slow  increase 
of  the  time  scale  of  motion  above  the  glass  transition 
is  closer  to  an  Arrhenius  temperature  dependence 
than  most  other  polymers.  Also  such  strong  systems 
typically  display  a  narrower  distribution  of  corre¬ 
lation  times.  The  segmental  motion  in  this  polymer 
is  thought  to  be  sterically  hindered  even  though  the 
glass  transition  temperature  is  quite  low.  The  poly¬ 
mer  has  a  low  permeability  to  gases  which  is  typically 
associated  with  the  low  specific  volume  in  this  sys¬ 
tem.  Molecular  packing  has  been  investigated  by 
Boyd^  using  Monte  Carlo  techniques;  these  authors 


INTRODUCTION 

As  a  polymeric  rubber,  poly  ( isobutylene )  ( PIB )  is 
quite  different  from  nearly  all  other  such  systems 
and  therefore  has  been  the  subject  of  ongoing  in¬ 
terest.  From  a  traditional  standpoint,  the  Williams- 
Landel-Ferry  ( WLF)  parameters  for  time-temper¬ 
ature  superposition  are  unusual.^  The  values  cor¬ 
respond  to  a  slow  increase  in  segmental  motion 
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also  investigated  penetrant  diffusion  using  molecular 
dynamics  simulations. 

Ih  a  previous  nuclear  magnetic  resonance  (NMR) 
study,®  the  time  scale  of  translational  diffusion  of 
CO2  in  PIB  on  a  local  scale  was  observed  to  be  close 
to  the  time  scale  of  the  segmental  motion  of  the 
polymer.  The  objective  of  this  report  is  to  more 
carefully  characterize  the  time  scales  of  segmental 
motion  in  PIB  using  NMR  spectroscopy  so  any  link 
to  the  time  scale  of  penetrant  diffusion  will  be  clearer. 
The  Kohlrausch -Williams -Watts  or  stretched  ex¬ 
ponential  was  used  to  characterize  the  distribution 
of  translational  correlation  times  for  the  CO2  and  it 
was  found  to  be  quite  narrow.  Comparable  infor¬ 
mation  on  the  segmental  motion  of  the  polymer  itself 
appeared  to  be  not  readily  available,  lending  addi¬ 
tional  motivation  for  this  study.  Based  on  the  strong 
liquid  characterization,  a  narrow  distribution  of 
segmental  correlation  times  would  however  be  ex¬ 
pected.^ 

Dejean  et  al.*^  have  used  carbon-13  Ti  and  NOE 
measurements  at  two  Larmor  frequencies  to  deter¬ 
mine  the  time  scale  of  segmental  motion  which  is  a 
well-established  technique  for  polymers  above  the 
glass  transition  temperature.  However  the  stretched 
exponential  was  not  used  in  the  interpretation. 
Rather  than  reinterpret  these  data,  new  measure¬ 
ments  were  performed  at  three  Larmor  frequencies 
to  extend  the  database.  Several  correlation  functions 
are  employed  in  this  report  to  interpret  the  new  data 
so  that  the  effect  of  the  interpretational  model  on 
the  time  scale  and  temperature  dependence  of  the 
motion  can  be  seen. 

In  addition  some  relaxation  measurements  were 
performed  on  PIB  under  a  pressure  of  CO2  compa¬ 
rable  to  that  employed  in  the  study  of  CO2  dynamics 
to  check  for  alteration  of  the  segmental  motion  of 
the  polymer  in  the  presence  of  the  gas.  Stuk®  has 
made  predictions  concerning  the  effect  of  a  gaseous 
penetrant  on  the  segmental  motion  in  a  rubbery 
polymer. 

With  these  new  measurements  and  the  associated 
interpretation,  a  more  accurate  comparison  of  the 
time  scale  of  segmental  motion  can  be  made  with 
the  time  scale  of  translational  diffusion  of  the  CO2 
in  the  polymer. 

EXPERIMENTAL 

The  polyisobutylene  sample  used  for  the  study  was 
obtained  from  Cellomer  Associates,  Inc.  of  Webster, 
NY  in  the  form  of  thick  sheets  of  high  molecular 
weight  (Cat.  #40E,  Lot  #02,  m.wt  1,000,000).  NMR 


experiments  were  carried  out  on  Bruker  SXP  (20- 
100),  AC-200  and  MSL-300  spectrometers,  which 
correspond  to  22.6  MHz,  50.3  MHz  and  75.5  MHz 
^®C  Larmor  frequencies,  respectively.  The  measure¬ 
ments  on  the  PIB  in  the  PIB-CO2  system  were  made 
on  the  AC  200. 

A  simple  180-t-90°  pulse  sequence  was  used  'to 
measure  the  spin-lattice  relaxation  time  ( Ti) .  Gated 
NOE  experiments  were  performed  on  the  AC-200 
and  on  the  SXP  (20-100).  High  power  decoupling 
was  used  for  the  measurements  on  the  MSL-300 
and  broadband  scalar  decoupling  on  the  AC-200  and 
SXP  (20-100).  Two  to  three  hours  have  been  al¬ 
lowed  for  the  sample  to  reach  a  thermal  equilibrium 
before  each  measurement. 

Proton  line  shape  experiments  were  performed 
on  the  MSL-300  using  a  solid  echo  sequence 
(90-t-90°). 


RESULTS 

Values  of  T 1  and  the- NOE  for  the  methylene  carbon 
in  PIB,  at  various  temperatures  and  Larmor  fre¬ 
quencies  are  presented  in  Figures  1  and  2.  As  can 
be  seen  from  the  figures  the  Ti  minimum  is  shifted 
to  higher  temperatures  as  the  field  strength  is 
increased®  in  agreement  with  expectations.  Figure 
1  also  shows  Tj  as  a  function  of  temperature  for  the 
methylene  carbon  in  PIB  in  the  PIB-CO2  system. 

Experimental  and  theoretical  line  shape  spectra 
are  shown  in  Figure  3.  Figure  4  shows  the  mobile 
population  of  pure  PIB  as  a  function  of  temperature. 


Figure  1.  Ti  vs.  temperature  at  22.6,  50.3,  and  75.5 
MHz.  Lines  correspond  to  the  fit  for  the  DLM  autocor¬ 
relation  function.  Fits  from  the  other  models  are  com¬ 
parable  to  the  fit  shown  for  the  DLM  model.  Squares  show 
the  Ti  vs.  temperature  at  50.3  MHz  for  the  PIB-CO2  sys¬ 
tem.  The  solid  line  corresponds  to  the  fit  for  the  Williams- 
Watts  function  combined  with  anisotropic  restricted  ro¬ 
tational  difiusion. 
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T(K) 

Figure  2.  NOE  vs.  temperature  at  22.6  and  50.3  MHz. 
Lines  correspond  to  the  fit  for  the  DLM  autocorrelation 
function.  Fits  from  other  models  are  comparable  to  the 
fit  shown  for  the  DLM  model. 


INTERPRETATION 

The  equations  relating  the  spin-lattice  relaxation 
time  and  the  spectral  densities  J(o))  for  a  methylene 
carbon-13  nucleus  are: 


—  =  +  W2 

1 1 


(hVJ,M 


^0  “  ~~ 


(CO2)  _ 

j-  0)2  -  0)H  ^ 


The  equation  for  the  NOE  is  given  by: 


^  .  Th  6J2(o)2)  ■“ 

NOE  =  1-1 - ;r-;^ — ,  :  — V  .  ;;v  "; — r 

yc  «A)(Wo)  "b  Oc/iCwc)  "T  oJ2{0)2) 

0)0  =  0)H  —  0)c  0)2  —  0)H 


The  methylene  C  —  H  distance  is  known  from  pre¬ 
vious  studies  to  be  1.09  A. 

A  number  of  expressions  for  the  spectral  density 
function  or  alternatively  the  correlation  function 
have  been  developed  for  segmental  motion  of  the 
polymer. 

One  of  them  is  the  Hall-Helfand^°  model  com¬ 
bined  with  anisotropic  restricted  rotational  diffu- 
sion.^^"^^  The  spectral  density  function  is  shown  be¬ 
low: 
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A  =  \  cos^A  -  1)^  B  =  I  (sin^2A) 

C  =  I  (sin^A)  (1) 

The  angle  A  is  between  the  internuclear  vector  and 
the  axis  of  rotation,  the  latter  being  the  carbon- 
carbon  backbone  bonds  for  this  case. 

In  this  description,  the  time  scale  of  the  segmental 
motion  is  set  by  two  parameters:  tq,  the  correlation 
time  for  single  conformational  transitions  and  ri, 
the  correlation  time  for  cooperative  or  correlated 
transitions.  Other  parameters  used,  are  the  activa¬ 
tion  energies  for  the  above  correlation  times  ( r 0 ,  ti  ) 
Eao  and  Eai,  respectively,  the  angular  amplitude  over 
which  rotational  diffusion  occurs,  Z,  and  the  rota- 
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Figure  3.  Proton  line  shapes  for  neat  PIB;  all  the  scales  are  in  ppm. 


tional  diffusion  constant  Av.  An  Arrhenius  temper¬ 
ature  dependence  was  given  to  the  correlation  times 
('^o»  'Ti)  and  the  temperature  dependence  of  the  re¬ 
stricted  diffusion (following  Gronski’s  example) 
is 

l  =  AT^'Udeg) 

The  preceding  spectral  density  function  gave  a 
satisfactory  fitting  based  on  the  following  parame¬ 
ters: 

I  =  S.IST®-®  (deg) 

Dir  =  3.26  X'lO’T  -  3.5  X  10®  (s”*) 

To  =  (T„)oexp(,Eao/RT)  Ti  =  {T^)iexp (E^i/ RT ) 


Figure  4.  Fraction  of  highly  mobile  population  for  bulk 
PIB  as  a  function  of  temperature. 
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where 

(r^)o  =  2.1  X  10“^^  s  (t«)i  =  1.0  X  10'^'  s 
Eqo  =  36.2  kJ  /mole  Ea\  =  36.7  kJ /mole 

The  reasonable  fits  are  obtained  for  activation 
energies  varying  between  29  to  37  kJ  /mole  and  the 
uncertainties  for  (too)o,  (roo)i,  are  a  factor  of  5.  The 
uncertainty  in  I  is  about  10  degrees  and  for  Dir  is 
about  30  10® 

The  data  were  also  fitted  using  a  WLF  temper¬ 
ature  dependence  for  the  correlation  times  instead 
of  an  Arrhenius  dependence: 


opment  of  this  correlation  function  was  in  response 
to  the  realization  of  the  importance  of  librational 
motion  as  a  source  of  spin-lattice  relaxation  in  rub¬ 
bery  systems  and  the  inclusion  of  librational  motion 
in  the  other  correlation  functions  reflects  this  finding 
by  Dejean,  Laupretre,  and  Monnerie.^®"^'^  The  cor¬ 
relation  time  Ti  is  associated  with  correlated  jumps 
responsible  for  orientation  diffusion  along  the  chain 
and  t2  corresponds  to  damping  which  consists  either 
of  nonpropagative  specific  motions  or  of  distortions 
of  the  chain  with  respect  to  its  most  stable  local 
conformations,  tq  is  associated  with  the  local  aniso¬ 
tropic  reorientation,  i.e.  libration.  The  spectral  den¬ 
sity  function  is  given  by: 


,  t{T)  -ciiT-To) 
°®'°r(To)  cl  +  T-To 


Ji 


_aro_\ 

1  +  0)^To) 


(2) 


where ^  Tq  =  298  K  and  C2  =  200.4  K  while  c?  and 
t(To)  were  floated.  This  time  an  acceptable  fitting 
was  obtained  for  the  following  parameters: 

to(298K)-5.4  10"®  s 
ti(298K)  =  3.4  10“®  s 

I  =  5  4^0.5 

A>  =  2.1  X  lO'^Tfs"') 
c?  =  5 

The  quality  of  the  fitting  is  the  same  as  the  one  with 
an  Arrhenius  temperature  dependence  presented 
before.  The  uncertainty  of  the  c?  is  about  0.3  and 
the  uncertainties  of  the  other  parameters  are  ap¬ 
proximately  the  same  as  before  (Arrhenius  fitting) . 

A  different  spectral  density  function  was  also  used 
to  fit  the  data.  This  is  the  Dejean-Laupretre-Mon- 
nerie  (DLM)  function^'^^^®"^’^.  The  autocorrelation 
function  is: 

G(t)  =  (1  -  a)exp(-t/r2)exp(-t/Ti)/o(t/ri) 

+  a  exp(-t/ro) 

where  Jo  is  the  modified  Bessel  function  of  order 
zero  and 


1  —  a  = 


’cos  6  —  cos®^P 
2(1  cos  0) 


where  0  is  a  cone  half-angle  in  which  random  an¬ 
isotropic  reorientation  of  the  CH  vector  occurs.^®  In 
this  model  the  local  chain  motion  consists  of  con¬ 
formational  jumps  and  bond  libration.  The  devel- 


where 

a  =  T2^  H-  2ti^T2^  -  0)^  =  --2a;(rr^  + 

In  this  following  Dejean,  Laupretre,  and  Monnerie^®"^^ 
case,  five  parameters  were  used  to  fit  the  data,  Ea, 
activation  energy,  Too  related  to  ri,  and  the  ratios 
ri/ro,  T2/ri.  The  fitting  gave  the  following  values: 

Ti  =  TooexpiEa/RT) 

Ea  -  35.3  kJ /mole 
Too  =  1.2  X  10"^^  s 
'Ti/tq  =  300 
T2/ri  =  3.1 
a  =  0.42 

In  this  model  the  activation  energy  also  lies  be¬ 
tween  29  and  36  kJ  /mole  and  the  uncertainty  for 
the  prefactor  correlation  time  is  a  factor  of  5.  For 
‘‘a”  the  uncertainty  is  0.1. 

Again,  a  WLF  temperature  dependence  can  be 
applied  to  the  correlation  time  ti  again  yielding  a 
suitable  description.  The  parameters  obtained  are 
the  following: 

c?  =  5 

ri(298K)  =  2.4  X  10’®  s 
ri/ro  =  250 
T2/ri  =  2.9 
a  =  0.42 

The  uncertainties  for  the  correlation  times  are  as 
before  and  the  uncertainty  in  c?  is  ca.  0,3. 

The  last  spectral  density  function  to  be  employed 
corresponds  to  the  Kohlrausch- Williams- Watts 
( KWW )  time  correlation  function 
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=  exp[-(t/Tp)“] 

combined  with  anisotropic  restricted  rotational  dif¬ 
fusion.  In  practice,  the  stretched  exponential  is 
written  as  a  sum  exponentials.^®  The  composite  au¬ 
tocorrelation  function  including  libration  that  will 
be  integrated  to  find  the  spectral  density  function, 
has  the  form^^’^® 


X  2 


1"/^—  cosjl  —  nir)  1  —  cos(l  -t  n7r)\® 
I  —  rnr  1  +  rnr  j 

sin(/  +  htt)'® 


/ sin(/  —  nir)  ^  si 
\  I  —  me 


I  -F  me 


<7"»(a))  +  C 


X2 

n=l  L 


cos(2Z  —  mr)  1  “  cos(2Z  + 


21  “•  mr 


21  +  mr 


^(t)  =  expi-t/TQ)^L(t) 
where  is  given 

B 

$i,(t)  =  A  +  —{  5  Ao(sin  4>)  +  ^  Ao(cos  <f)) 
+  2  [A®(sin  0) -t-i4®(cos 


+  —  {  j  i4o(sin  20)  -t-  Ao(cos  20) 


+  2  [i4®(sin  20) -t  i4®(cos  20)]e 

n-l 

where  the  A,  B,C  coefficients  have  been  defined  in 
eq.  ( 1 ) .  Angle  A  has  also  been  defined  earlier  and 
is  equal  to  71®  and^^*^^ 

^n{f{4>))  ^  J  /(0)sin^-^^j  def)  0  ^  (j)  <  I 

The  resulting  spectral  density  function  is: 

f  2B 

c/Cw)  =  j  2 A  +  [(1  “  cos  Z)^  +  sin^Z] 

C  1 

+  ^  [(1  —  cos  2Z)^  +  sin^2Z]  +  B 


+ 
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= 


To 


1  +  OjVq 

n7r\^ 

t)"’" 


.?»■(»)  - 


K  +  To 


-1 


(Xn  +  T5-^)^  +  a; 


In  this  case  six  parameters  were  adjusted;  the  ac¬ 
tivation  energy  associated  with  tq,  Too,  the  angular 
amplitude  Z  over  which  rotational  diffusion  occurs, 
the  rotational  diffusion  (two  parameters),  and  a, 
the  fractional  exponent.  The  values  that  were  found 
are: 


-  1.4  X  10"^^  s 
Ea  =  34.7  kJ/mole 
a  =  0.98 

Z  =  5.18r°-^  (deg) 

Dir  =  1.7  X  lO^T  -  6.5  X  10®  (s'^) 

It  must  be  noted  that  Xp  is  the  characteristic  cor¬ 
relation  time  and  sets  the  position  of  the  distribution 
associated  with  the  stretched  exponential  on  the 
time  axis.  The  parameter  a,  which  is  between  0 
<  a  <  1,  determines  the  width  of  the  distribution. 
The  fitting  was  done  by  determining  first  the  central 


Table  I.  Proton  Line  Shape  Parameters 


Temperature 

(K) 

Gaussian  Width 
(kHz) 

Lorentzian  Width 
(kHz) 

Fraction  Gaussian 

Fraction  Lorentzian 

200 

49.0 

20.0 

0.90 

0.10 

210 

48.0 

20.0 

0.86 

0.14 

220 

48.0 

20.0 

0.75 

0.25 

230 

44.0 

20.0 

0.65 

0.35 

240 

37.0 

18.0 

0.50 

0.50 

250 

25.0 

15.0 

0.40 

0.60 

260 

20.0 

11.0 

0.30 

0.70 

270 

12.0 

8.0 

0.25 

0.75 

299 

4.0 

3.3 

0.10 

0.90 
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value  of  Tp,  for  the  distribution.  Next  a  series  of  tq 
values  on  both  sides  of  the  central  value  was  cal¬ 
culated.  A  total  of  81  To’s  were  used  for  this  fitting. 
Then  for  every  tq,  a  relaxation  time  Ti  is  calculated 
and  in  order  to  acquire  the  total  relaxation  time  the 
individual  rates  are  weighted.^® 

The  uncertainty  for  Tqo  ,  is  again  a  factor  of  five 
and  the  activation  energy  lies  between  29  and  36 
kJ/mole;  the  uncertainty  for  a  is  0.05,  for  I  is  ap¬ 
proximately  10°  and  for  Dir  about  30  X  10®  s^^. 

A  WLF  temperature  dependence  can  again  be 
applied  to  the  correlation  time  Tp  and  the  associated 
fitting  parameters  are: 

Tp  (298  K)  =  2.1  X  10'®  s 
Cl  =  5 
a  =  0.982 
/  =  5T®®(deg) 

1.7  X  10^T(s”') 

Once  more  the  parameters  have  not  changed  very 
much  from  those  obtained  before  and  the  quality  of 
the  fitting  is  approximately  the  same.  The  uncer¬ 
tainties  are  the  same  as  before  for  all  the  parameters. 

The  last  spectral  density  function  (Williams- 
Watts  combined  with  anisotropic  rotational  diffu¬ 
sion)  was  also  used  for  fitting  the  Ti  data  for  PIB 
in  the  PIB-CO2  system.  In  this  case,  all  the  param¬ 
eters  were  kept  the  same  as  before,  and  a  was  the 
only  parameter  that  was  floated.  The  value  for  which 
the  best  fit  was  obtained,  is 


component.  Three  parameters  are  involved:  the 
width  of  the  Gaussian  component,  the  width  of  the 
Lorentzian  component,  and  the  fractional  popula¬ 
tion  of  the  Gaussian  component.  The  values  ob¬ 
tained  are  shown  in  Table  I. 

For  relating  polymer  dynamics  to  line  shapes, 
Rossler’s  formalism^^  of  dividing  molecular  motion 
into  two  categories,  fast  and  slow,  has  gained  con¬ 
siderable  popularity.  Rossler  employed  a  log-Gauss 
distribution  of  correlation  times  which  corresponds 
to  a  Gaussian  function  for  the  energy  distribution 
giE),  There  is,  however,  a  technical  problem  with 
Rossler’s  use  of  the  log-Gaussian  probability  density 
g{E)  in  the  case  of  random  variables  which  are  non¬ 
negative.  The  log  Gaussian  is  defined  for  real  vari¬ 
ables  which  are  positive  or  negative.  Activation 
energies  are  positive  only,  and  the  Gaussian  corre¬ 
sponds  to  an  unphysical  model  in  this  case.  A  simple 
alternative  distribution  is  the  gamma  density,  de¬ 
fined  for  positive  random  variables  only,  and  the 
modifications  to  Rossler’s  method  are  sketched 
briefly  below. 

We  introduce  a  gamma  density  distribution  of 
energy  barriers 


p{E)- 


qigEV-^e-"^ 

T(a) 


where  oc  and  q  are  two  adjustable  parameters  and 
r(  a )  is  a  gamma  function  of  a.  It  follows,  then,  that 
the  population  with  activation  energy  less  than  E* 
is  given  by: 


a  =  0.91 

In  this  case  one  parameter  was  used  to  fit  eight  ex¬ 
perimental  data  points.  As  an  alternative  approach, 
the  activation  energy  was  floated  while  a  was  kept 
constant  at  0.98  as  in  neat  PIB.  The  optimum  fitting 
yielded 

Ea  =  34.3  kJ /mole 

which  is  approximately  the  same  (within  the  un¬ 
certainty)  of  the  previous  fitting,  for  the  neat  PIB. 

For  all  the  fittings  a  nonlinear  least-squares  fit¬ 
ting  program  was  employed  to  treat  the  data.  The 
resulting  fits  have  root  mean  square  residuals  around 
10%  or  less,  which  is  also  the  approximate  value  for 
the  experimental  error. 

Finally  the  proton  line  shapes  of  pure  PIB,  shown 
in  Figure  3,  were  fitted  using  three  parameters.  The 
line  shapes  were  described  as  a  combination  of  a 
broad  Gaussian  component  and  a  narrow  Lorentzian 


p{E)  dE 

0 

=  iqEV-^e-'^^dE  (2a) 

r{a)  Jo 

If  we  use  Rossler’s  relation,  E*  =  RT  ln(TV tq), 
then  we  can  write  equation  ( a )  as: 

Pfast(T)  =  U[a,qRTln(ryro)]  (2b) 

where  11  in  eq.  (2b)  is  the  incomplete  gamma  func¬ 
tion,  1/to  is  the  attempt  frequency,  and  r*  is  the 
characteristic  correlation  time  determined  by  the 
width  of  the  line.  The  mean  activation  energy  is 
given  by: 

(E*)  =  a/q 

The  fitting  of  the  population  of  the  Lorentzian 
or  fast  component  of  the  proton  line  shapes  as  a 
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function  of  temperature  is  shown  in  Figure  4.  The 
parameters  obtained  are  a  =  50  and  <?  ==  6.4  X  10"^ 
cai"^  with  a  mean  activation  energy  of  (JS*  )  =  32.7 
kJ/mole.  The  above  values  were  obtained  when  r* 
==  10  s  and  tq  =  10“^^  s  as  suggested  by  Rossler. 
The  choice  of  these  characteristic  correlation  times 
is  somewhat  arbitrary  and  varying  them  over  rea¬ 
sonable  ranges  changes  the  estimate  of  the  average 
energy  from  26  to  44  kJ /mole. 


DISCUSSION 

Any  of  the  three  correlation  functions  or  equiva¬ 
lently  the  three  spectral  densities  provide  a  suitable 
basis  for  the  interpretation  of  the  and  NOE  data 
taken  on  PIB  as  a  function  of  temperature  and  Lar- 
mor  frequency.  All  of  the  interpretational  functions 
also  produce  comparable  estimates  of  the  time  scale, 
activation  energy,  and  the  effective  breadth  of  the 
distribution  of  exponential  correlation  times.  Thus 
the  view  of  segmental  motion  developed  for  PIB  is 
not  particularly  sensitive  to  the  interpretational 
model.  The  Rossler  formulation  which  is  also  based 
on  the  existence  of  a  distribution  of  correlation  times 
enables  one  to  treat  the  line  shape  data,  where 
a  correlation  function  approach  is  not  easily  avail¬ 
able,  and  thus  make  a  comparison  with  the  anal¬ 
yses.  In  this  case,  reasonable  agreement  with  the 
activation  energies  is  obtained. 

In  agreement  with  the  results  of  Dejean,  Lau- 
pretre,  and  Monnerie, librational  motion 
must  be  included  to  interpret  the  Ti  data.  The 
Gronski  description  of  librational  motion  when 
combined  with  either  the  Hall-Helfand  model  or 
the  Kohlrausch-Williams-Watts  (KWW)  modeP®'^° 


Figure  5.  /  and  Z),>  as  a  function  of  temperature;  values 

from  the  Hall-Helfand  model. 


produces  a  fit  of  the  data  which  is  comparable  to  the 
Dejean-Laupretre-Monnerie  approach.  The  Gronski 
model  assumes  rotation  about  a  single  axis  and  leads 
to  fairly  large  amplitudes  for  the  librational  motion. 
As  can  be  seen  from  Figure  5,  the  amplitude  in¬ 
creases  from  91°  at  313  K  to  102°  at  383  K.  The 
cone  angle  in  the  DLM  model  is  smaller  and  possibly 
more  realistic.  The  DLM  model  treats  libration  as 
an  activated  process  while  the  Gronski  approach  is 
just  motion  in  the  bottom  of  a  square  well.  With 
respect  to  this  latter  point,  the  view  of  librational 
motion  as  an  activated  process  is  not  so  realistic. 

The  distribution  of  exponential  correlation  times 
is  very  narrow  whichever  model  is  considered  though 
this  is  most  apparent  in  the  context  of  the  KWW 
model  where  the  exponent  approaches  unity  and 
thus  the  distribution  collapses  to  a  single  exponen¬ 
tial.  The  narrow  distribution  is  in  agreement  with 
the  strong  liquid  description  of  Angell  et  al.^  noted 
in  the  introduction  and  the  NMR  data  presented 
here  is  one  of  the  better  pieces  of  information  avail¬ 
able  in  the  literature  on  the  distribution  of  corre¬ 
lation  times  in  rubbery  PIB. 

In  relation  to  the  original  purpose  of  this  study, 
segmental  motion  in  PIB  has  a  time  scale,  width  of 
distribution,  and  activation  energy  which  are  close 
to  those  observed  for  the  translational  diffusion  of 
CO2  sorbed  in  PIB.  Table  II  shows  a  direct  com¬ 
parison  of  the  correlation  times  for  segmental  mo¬ 
tion  of  PIB  and  the  correlation  times  for  transla¬ 
tional  diffusion  of  CO2  sorbed  in  PIB  at  the  same 
temperatures.  The  coincidence  of  the  time  scales  of 
these  two  types  of  motion  does  not  always  occur  for 
gases  in  rubbery  polymers  since  in  the  only  other 
system^®  [CO2  sorbed  in  poly  ( dimethyl  siloxane)] 
where  both  motions  were  studied  by  NMR,  segmen¬ 
tal  motion  is  orders  of  magnitude  faster  than  trans¬ 
lational  diffusion.^®  It  is  tempting  to  cite  the  low 
specific  volume  in  PIB  as  the  source  of  the  connec¬ 
tion  between  the  time  scale  of  translational  diffusion 
of  a  penetrant  and  segmental  motion  in  PIB.  In 
other  polymers  where  less  efficient  packing  is  pres¬ 
ent,  a  translation  jump  of  a  penetrant  could  occur 
between  two  local  regions  of  high  free  volume  un¬ 
correlated  to  any  conformational  event  in  the  poly¬ 
mer  backbone.  In  the  well-packed  PIB,  the  pene¬ 
trant  might  move  only  when  the  chain  backbone 
moves.  This  picture  can  only  be  verified  through  the 
measurement  of  the  two  types  of  motion  in  a  series 
of  polymer/penetrant  systems.  It  is  clear,  however, 
that  PIB  behaves  quite  differently  from  other  rub¬ 
bery  polymers  both  in  terms  of  penetrant  diffusion 
as  well  as  other  behaviors. 
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Table  II.  {tp)bd  from  Ref.  (6)  and  Xp  from  the  Kohlrausch-Williams-Watts  Function 
Combined  with  Anisotropic  Restricted  Rotational  Diffusion,  As  a  Function  of 
Temperature 


Central  Correlation  Time  for 
Translational  Diffusion  of  CO2  in 
the  PIB-CO2  System 

Central  Correlation  Time 
for  Segmental  Motion  of 
PIB 

T(K) 

(tp)dd  (sec) 

Tp  (sec) 

313 

1.56  10-® 

8.65  10"® 

323 

1.08  10“® 

5.73  10"® 

333 

7.64  10'® 

3.88  10"® 

343 

5.51  10"® 

2.70  10"® 

353 

4.05  10"® 

1.91 10"® 

363 

3.02  10"® 

1.38  10"® 

373 

2.29  10"® 

1.01 10"® 

383 

1.76  10"® 

7.56  10"‘® 

A  limited  number  of  measurements  of  spin-lattice 
relaxation  times  were  also  made  on  PIB  in  the  pres¬ 
ence  of  sorbed  CO2  and  rather  small  changes  were 
observed  from  the  values  obtained  on  pure  PIB.  The 
position  of  the  Ti  minimum  is  unchanged  though 
perhaps  the  distribution  of  exponential  correlation 
times  is  slightly  broader  relative  to  the  data  on  pure 
PIB.  The  deviation  observed  at  313K  in  Figure  1  is 
not  significant  since  the  precision  of  the  Ti  mea¬ 
surement  is  less  here  due  to  the  approach  to  the 
NMR  shifted  and  the  resulting  increase  in  line 
width.  Thus  on  average,  the  motion  of  PIB  is  not 
greatly  altered  by  the  presence  of  the  relatively  small 
amount  of  sorbed  CO2,  a  couple  of  percent  by  weight. 
This  is  to  be  expected;  however,  Stuk®  did  estimate 
a  significant  shift  of  the  segmental  correlation  time 
locally  in  the  vicinity  of  the  penetrant  relative  to 
the  correlation  time  of  the  bulk.  In  that  picture  the 
correlation  time  locally  for  segmental  motion  and 
also  for  translational  diffusion  would  be  considerably 
shorter  than  the  correlation  time  for  segmental  mo¬ 
tion  in  the  bulk  polymer.  Stuk  calculates  this  in  the 
form  of  a  temperature  shift  of  45  degrees  for  the 
segmental  motion  correlation  time.  According  the 
correlation  times  listed  in  Table  II  this  temperature 
shift  corresponds  to  a  time  scale  change  of  a  factor 
of  3  or  4.  This  can  be  considered  relative  to  the  in¬ 
terpretation  given  here  by  comparing  the  correlation 
time  for  translational  diffusion  which  is  determined 
by  observing  CO2  which  has  locally  increased  seg¬ 
mental  motion  with  the  bulk  segmental  motion  cor¬ 
relation  time.  As  seen  in  Table  II,  the  correlation 
time  for  the  diffusion  of  CO2  is  within  a  factor  of  2 
of  the  correlation  time  for  segmental  motion  but  it 


is  actually  longer  rather  than  shorter.  Thus  these 
correlation  times  for  the  two  processes  generated 
from  spin-relaxation  data  do  not  support  the  Stuk 
prediction  at  first  glance.  However,  in  a  simple  model 
of  spheres  undergoing  rotational  and  translational 
diffusion  with  a  Stokes'  viscosity  dependence,  the 
translational  correlation  time  is  nine  times  longer 
than  the  rotational  correlation  time.^^  If  this  ratio 
were  applicable  in  the  current  case,  then  local  seg¬ 
mental  motion  in  the  vicinity  of  a  CO2  molecule 
would  appear  to  be  three  to  four  times  faster  in  good 
agreement  with  the  Stuk  prediction.  There  is  no 
clear  model-based  prediction  for  the  ratio  of  trans¬ 
lational  to  rotational  correlation  times  given  the 
complexity  of  the  system  and  the  correlation  func¬ 
tions  used  for  interpretation,  so  the  ratio  based  on 
diffusing  spheres  is  the  best  estimate. 

Acceptable  interpretations  of  the  segmental  mo¬ 
tion  in  PIB  or  the  motion  of  CO2  through  the  PIB 
are  produced  if  a  WLF  temperature  dependence  is 
used  in  place  of  the  Arrhenius  description.  Figure  6 
shows  a  plot  of  the  correlation  time  tq  for  segmental 
motion  as  a  function  of  temperature  that  are  ob¬ 
tained  by  fitting  the  data  with  the  two  descriptions. 
Over  the  temperature  range  of  the  measurements, 
either  an  Arrhenius  or  a  WLF  type  description  led 
to  comparable  set  of  correlation  times.  The  corre¬ 
lation  times  predicted  for  lower  temperatures  outside 
the  range  of  experimental  observation  by  this  tech¬ 
nique  are  substantially  different  and  the  WLF  de¬ 
pendence  certainly  is  more  likely  valid  over  the  wider 
range. 

The  WLF  parameters  obtained  from  the  NMR 
data  do  not  agree  with  those  obtained  from  visco- 
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1/T(K) 

Figure  6.  Comparison  of  the  temperature  dependence 
of  To  from  the  Hall-Helfand  model  combined  with  re¬ 
stricted  anisotropic  rotational  diffusion  assuming  first 
WLF  and  then  Arrhenius  behavior.  The  lower  two  curves 
compare  the  Arrhenius  and  the  WLF  temperature  depen¬ 
dence  of  To  produced  by  fitting  Ti  when  c?  was  floated 
(ci  =  5);  the  two  curves  coincide  at  high  temperatures. 
The  upper  two  curves  show  the  Arrhenius  and  the  WLF 
temperature  dependence  of  to  produced  from  fitting  Ti 
when  Cl  was  kept  constant  and  equal  to  a  literature  value  ^ 
(c?  =  8.61 ) ;  the  two  curves  do  not  coincide  at  high  temper¬ 
atures  and  it  was  impossible  to  produce  an  acceptable  fit. 


elastic  experiments.  In  the  NMR  interpretations 
developed  here,  c?  is  about  5  while  literature  values 
are  around  8,6.^  The  NMR  data  coUld  not  be  fit  using 
the  literature  values  and  Figure  6  shows  a  compar¬ 
ison  of  one  of  the  Arrhenius  fits  with  the  WLF  be¬ 
havior  given  by  c?  =  8.6  setting  the  two  dependencies 
to  a  common  value  for  the  correlation  time  at  a  tem¬ 
perature  of  300K.  Dejean  et  al.^  reported  a  similar 
inability  to  reconcile  the  temperature  dependence 
of  the  correlation  times  from  NMR  with  the  ac¬ 
cepted  values  of  the  WLF  parameters  for  PIB.  For 
other  polymers,  NMR  data  lead  to  the  same  values 
for  the  WLF  parameters  as  are  obtained  by  other 
techniques.^®  Here  again  PIB  may  be  an  exception 
to  the  rule. 

If  the  information  from  NMR  on  the  temperature 
dependence  of  the  translational  diffusion  of  CO2  is 
compared  with  literature  values,  inconsistencies 
again  arise.  According  to  the  results  of  computer 
simulation  for  the  diffusion  of  methane  in  PIB  by 
Pant^^  in  the  temperature  range  350-600K,  an  ac¬ 
tivation  energy  for  diffusion  of  39  kJ/mole  was  ob¬ 
tained.  While  experimental  values  for  methane  in 
PIB  were  not  available,  activation  energies  for  ox¬ 
ygen,  nitrogen,  and  carbon  dioxide  obtained  from 
permeation  data  are  in  the  range  of  50.2  ±  0.4  kJ / 
mole.  The  activation  energy  of  28.2  kJ/mole  for  the 
diffusion  of  CO2  in  PIB  obtained  from  an  NMR 


study®  is  closer  to  Pant’s  simulation  result  of  39  kJ / 
mole  than  that  of  50.2  kJ  /mole  obtained  by  gas  per¬ 
meation  experiments.  Again  the  values  are  not  very 
consistent  and  it  is  not  clear  how  to  mesh  the  results 
from  different  techniques  into  an  overall  under¬ 
standing. 
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Introduction 

In  a  previous  NMR  studyCU,  the  time  scale  of  translational 
diffusion  of  CO2  in  PIB  on  a  local  level  was  observed  to  be  close  to  the 
time  scale  of  the  segmental  motion  of  the  polymer.  The  objective  of  this 
report  is  to  more  carefully  characterize  the  time  scales  of  segmental  motion 
in  PIB  using  NMR  spectroscopy  so  that  link  to  time  scale  of  penetrant 
diffusion  may  be  clarified.  The  Williams  -  Watts  or  stretched  exponential 
was  used  to  characterize  the  distribution  of  translational  correlation  times 
for  the  CO2  and  it  was  found  to  be  quite  narrow.  Comparable  information 
on  the  segmental  motion  of  the  polymer  itself  appeared  not  to  be  readily 
available  lending  additional  motivation  for  this  study.  Based  on  the  strong 
liquid  characterization,  a  narrow  distribution  of  segmental  correlation 
times  would  however  be  expected(2). 

Experimental 

The  PIB  sample  used  for  this  study  was  obtained  from  Cellomer 
Associates,  Inc.  of  Webster,  NY  in  the  form  of  thick  sheets  of  high 
molecular  weight  (Cat.  #40E,  Lot  #02,  mol.  wl.  1,000,000).  NMR 
experiments  were  carried  out  on  Bruker  SXP  (20  -  100),  AC  -  200  and 
MSL  -  3(X)  spectrometers,  which  correspond  to  22.6  MHz,  50.3  MHz 
and  75.5  MHz  i3c  Larmor  frequencies  respectively.  The  measurements 
on  the  PIB  -  CO2  system  were  made  on  the  AC  -  200. 

Results 

Values  of  Ti  for  the  methylene  carbon  in  neat  PIB  and  PIB  -  CO2 
system,  at  various  temperatures  and  Larmor  frequencies  are  given  in  Table 
I.  The  NOE  as  a  function  of  temperature  and  Larmor  frequency  is  shown 
in  figure  1. 


j  2  It  40  fj 
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The  equation  for  the  NOE  is  given  by: 

Yu  6  J,(co,)  -  J„(o)„) 

NOE  =  1  +-I!i - ^  ^  °  ° - 

Yc  +  3  J,((0q)  +  6  J2(C02) 

where  cOq  =  -  03c  and  cOo  =  coh  +  (Oc-  The  methylene  C  -  H  distance  is 

known  from  previous  studies  to  be  1.09  A. 

A  number  of  expressions  for  the  spectral  density  function  or 
alternatively  the  correlation  function  have  been  developed  for  segmental 
motion  of  the  polymer.  Three  correlation  functions  arc  used  in  this  study 
to  treat  the  relaxation  data  and  describe  the  local  motion  of  the  PIB 
backbone.  The  first  is  the  Hall  -  Helfand  function  combined  with 
anisotropic  libration,  the  second  is  the  Dejean  -  LauprStre  -  Monnerie 
function  and  the  third  is  the  Williams  -  Watts  function  combined  with 
libration: 

<I>(t)~exp[-(t/tp)o 

In  practice,  the  stretched  exponential  is  written  as  a  sum  exponentials.  The 
composite  autocorrelation  function  including  libration  that  will  be 
integrated  to  find  the  spectral  density  function,  has  the  form 

0(t)  =  cxp(-t/to)<I>L(t) 

where  C>L(t)  is  given  by 

®  -  A  +  ~  { i  >^(sin4.)  + 1  aJ(cos^)  +  X I  A*(sin4>)  +  A^(cos4i)  1  e  . 

+  7  ( J  Aj(sin2^)  +  >^(cos2<|»)  +  X  ( A^(sln2(>)  +  A^(cos24»)  ]  e  ^  ) 
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The  resulting  spectral  density  function  is: 

J(co)«{2A  +  -^[(1  -cosC)^+sin^I  +  -~[(1  -  cos2f)*  +  sin*2C  ] )  J°((o)  + 


+  B  ^  1  -005(1 +  n7c)^2  ^  ^  sin{C-»  ^ 

n-i  C-njc  C  +  nn  C-n«  t  +  njc 


1 -005(2^11 «)  ^  1-cos(2t .1)11)  8  sln(2t-nB)  ^  sln(2C.f  n«)j2 


Interpretation 

The  equations  relating  the  spin  -  lattice  relaxation  time  and  the 
spectral  densities  J(co)  for  carbon  - 13  nucleus  are: 


- Tz — 2 


In  this  case  six  parameters  were  adjusted;  the  activation  energy 
associated  with  Tq,  x«,  the  angular  amplitude  I  over  which  restricted 
rotational  diffusion  occurs,  the  rotational  diffusion  (two  parameters)  and 


a,  the  fractional  exponent  The  values  are: 


It.  =  1.4  10‘14  sec 
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62.40 

81.50 


83.20 

86.85 


133.60 

149.90 


83.35 


Ea  =  34.7  kJ/ mole 

a  =  0.98 

t  =  5.18  1^*5 

Dir  =  1.7  107  T- 6.5  108 

Discussion 

Any  of  the  three  correlation  functions  provide  a  suitable  basis  for 
the  interpretation  of  the  Tj  and  NOE  data  taken  on  PIB  as  a  function  of 
temperature  and  Larmor  frequency.  All  of  the  interpretational  functions 
also  produce  comparable  estimates  of  the  time  scale,  activation  energy  and 
the  breadth  of  the  distribution  of  exponential  correlation  times.  Thus  the 
view  of  segmental  motion  developed  for  PIB  is  not  particularly  sensitive 
to  the  interpretational  model 

The  distribution  of  exponential  correlation  times  is  very  narrow 
whichever  model  is  considered,  though  this  is  most  apparent  in  the 
context  of  the  Williams  -  Watts  model  where  the  exponent  approaches 
unity  and  thus  the  distribution  collapses  to  a  single  exponential.  The 
narrow  distribution  is  in  agreement  with  the  strong  liquid  description  of 
Angell  and  the  NMR  data  presented  here  is  one  of  the  better  pieces  of 
information  available  in  the  literature  on  the  distribution  of  correlation 
times  in  rubbery  PIB. 

In  relation  to  the  original  purpose  of  this  study,  segmental  motion 
in  PIB  has  a  time  scale,  width  of  distribution  and  activation  energy  which 
is  close  to  that  observed  for  the  translational  diffusion  of  CO2  sorbed  in 
PIB.  Table  11  shows  a  direct  comparison  of  the  correlation  times  for 
segmental  motion  of  PIB  and  the  correlation  times  for  translational 
diffusion  of  CO2  sorbed  in  PIB  at  the  same  temperatures.  The  coincidence 
of  the  time  scales  of  these  two  types  of  motion  does  not  always  occur  for 
gases  in  rubbery  polymers,  since  in  the  only  other  system  (CO2  sorbed  in 
PDMS)  where  both  motions  were  studied  by  NMR,  segmental  motion  is 
orders  of  magnitude  faster  than  the  translational  diffusion^). 
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TABLE  I:  Methylene  carbon  Ti  values  for  neat  PIB  and  for  PIB  /  CO2 


system  as  a 
frequencies. 

function  of 

temperature  at 

various  Larmor 

T(K) 

Ti  (ms) 

Ti  (ms) 

Ti  (msec) 

Ti  (msec) 

22  MHz 

50.3  MHz 

75.4  MHz 

50.3  MHz 

Neat  PIB 

Neat  PIB 

Neat  PIB 

PIB/CO2 

313 

44.50 

135.66 

213.31 

103.16 

323 

41.37 

98.79 

157.99 

86.56 

333 

41.10 

71.06 

127.67 

79.58 

343 

42.47 

68.04 

113.39 

73.64 

353 

53.29 

71.04 

112.98 

74.95 

363 

64.40 

74.48 

117.90 

79.50 

373 

383 


TABLE  H:  (Xp)DD  from  reference  (1)  and  Tp  from  the  Williams  -  Watts 
function  combined  with  restricted  anisotropic  rotational 
diffusion  as  a  function  of  temperature. 


T(K) 

(tp)DD  (sec) 

Tp(scc) 

313 

1.56  10*8 

8.65  10-9 

323 

1.08  10-8 

5.73  10-9 

333 

7.64  10-9 

3.88  10-9 

343 

5.51  10-9 

2.70  10-9 

353 

4.05  10-9 

1.91  10-9 

363 

3.02  10-9 

1.38  10-9 

373 

2.29  10-9 

1.01  10-9 

383 

1.76  10-9 

7.56  10-10 
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Figure  1.  NOE  vs  Temperature  at  22.6  MHz  and  50.3  MHz.  Lines 
correspond  to  the  fit  for  the  DLM  autocorrelation  function. 
Fits  from  other  models  are  comparable  to  the  fit  shown  for  the 

DLM  piodel. 
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Abstract 

The  self  -  diffusion  coefficients  of  toluene  in  polyisobutylene  (PIB)  solutions  were 
determined  using  the  Pulsed  Field  Gradient  Nuclear  Magnetic  Resonance  technique.  The  volume 
fraction  of  toluene  in  the  polymer  was  varied  from  0.045  up  to  0.712  and  the  temperature  was 
varied  from  225  K  up  to  368  K.  The  concentration  dependence  of  the  data  was  interpreted  using 
the  Fujita  Free  Volume  theory  and  the  temperature  dependence  was  interpreted  with  the  WLF 
equation.  These  models  describe  separately  the  concentration  and  temperature  dependencies  of 
the  toluene  self  -  diffusion  coefficients  very  well  and  the  resulting  free  volume  parameters  are  in  a 
good  agreement  with  the  ones  extracted  from  the  analysis  of  viscosity  data  on  the  same  system. 


Keywords:  Polyisobutylene,  Toluene,  Pulse  Field  Gradient  Nuclear  Magnetic  Resonance, 
Diffusion,  Fujha  Free  Volume  Theory. 


Introduction 


The  diffusion  of  small  molecules  into  polymers  has  been  the  subject  of  many 
investigations,  because  of  its  importance  in  industrial,  medical  and  biological  applications.  Some 
examples  are  the  membrane  separation  of  gases,  controlled  drug  delivery,  control  of 
polymerization,  production  of  barrier  materials  for  the  packaging  industry  and  chemical  defense. 
Diffusion  of  penetrants  in  polyisobutylene  is  of  special  interest  since  this  material  displays  low 
permeability  for  a  polymeric  rubber  well  above  the  glass  transition.  This  unusual  property  of 
polyisobutylene  has  been  considered  in  computer  simulations(^)  of  the  rubber  and  in  a  detailed 
characterization  of  segmental  motion(2)  in  the  pure  rubber  by  NMR. 


There  are  a  variety  of  techniques  for  measuring  the  self  -  diffusion  coefficients  of  a 
penetrant  into  a  polymer  matrix.  These  include  sorption,  permeation,  light  scattering,  radioactive 
tracing,  forced  Rayleigh  scattering  and  Nuclear  Magnetic  Resonance  (NMR).  The  Pulsed  Field 
Gradient  NMR  Technique  is  valuable  because  it  does  not  perturb  the  system.  With  this  technique 
a  direct  measure  of  the  self  -  diffusion  coefficient  of  the  penetrant  is  achieved  by  observing  the 
molecules  microscopically,  while  other  methods  (sorption  method  for  example)  indirectly 
determine  the  self  -  diffusion  coefficient  from  macroscopic  measurements.  Also  the  experiment 
can  give  information  about  the  geometry  of  the  surroundings  of  the  penetrant  molecules  in  the 
case  of  heterogeneous  systems;  the  dimension  of  the  boundaries  that  the  molecules  may 
encounter  and  the  mean  square  distance  that  the  penetrant  has  travelled  during  the  observation 
of  diffusion.  Other  advantages  of  the  technique  are  the  ability  to  use  small  sample  volumes,  the 
precision  of  the  measurements  and  the  ability  to  be  applied  to  both  liquid  and  solid  samples.  The 
only  limitation  of  the  technique  is  that  the  lowest  diffusion  coefficient  that  can  be  commonly 
determined  is  around  1 0'®  cm^  /  sec  to  1 0'^  ®  cm^  /  sec. 

In  addition  to  the  advantages  of  Pulsed  Gradient  technique  itself,  with  NMR  the 
reorientational  dynamics  of  the  both  the  polymeric  and  penetrant  components  of  the  solution  can 
be  measured  independently  on  the  same  sample.  Thus  both  translational  and  rotational  motion  of 
the  penetrant,  and  segmental  motion  of  the  polymer  can  be  determined  in  the  same  samples. 
With  modern  NMR  instrumentation,  these  measurements  can  be  made  over  a  wide  range  of 
temperatures  and  concentrations  which  lead  to  a  fairly  complete  description  of  local  dynamics.  It 
is  the  goal  of  the  series  of  papers  given  here  to  present  the  results  of  such  a  fairly  complete  set  of 
NMR  measurements  on  a  single  system  beginning  with  translational  diffusion  of  the  penetrant. 

Diffusion  in  polyisobutylene  (PIB),  has  been  studied  before  using  different  techniques(^  ' 
^).  In  references  (3)  and  (4)  the  study  of  the  diffusion  coefficients  of  six  hydrocarbons  (propane,  n 
-  butane,  isobutane,  n  ■  pentane,  isopentane  and  neopentane)  in  PIB  was  reported  using  the 
sorption  method.  The  self  -  diffusion  of  benzene  in  PIB  has  been  studied  in  reference  (5)  using 
the  pulsed  field  gradient  technique  over  a  limited  range  of  temperatures  and  concentrations.  In 
reference  (6)  the  diffusion  of  radioactively  tagged  cetane  in  PIB  has  been  studied  and  finally  in 
reference  (7)  the  diffusion  of  toluene  in  butyl  rubber  has  been  studied  using  the  sorption  method. 
In  all  cases  the  data  were  analyzed  using  the  Fujita  free  volume  theory.  Viscosity  data  on 
polyisobutylene  solutions  are  also  available(® '  In  references  (8-10)  the  intrinsic  viscosities 
of  PIB  solutions  in  various  solvents  have  been  studied  and  their  relations  to  polymer  chain 
structure  and  to  thermodynamic  parameters  governing  the  interaction  between  the  polymer  and 
the  solvent  have  been  made. 


Free  volume  theories  have  been  used  widely  for  the  interpretation  of  the  concentration 
dependence  of  the  self  -  diffusion  of  many  small  molecules  into  polymer  matrices(''2  - 15).  One 
frequently  used  theory,  developed  almost  thirty  years  ago(^2)^  is  the  Fujita  theory  which  will  be 
employed  in  this  paper  to  summarize  penetrant  motion.  In  the  second  report  of  this  series,  the 
same  penetrant  data  will  be  reanalyzed  with  a  more  recent  free  volume  theory,  that  of  Vrentas 
and  Duda.  These  two  free  volume  theories  differ  in  several  respects  and  it  will  be  of  value  to 
compare  them  in  the  context  of  a  rather  complete  data  set  not  only  including  information  of 
translational  motion  of  the  penetrant  but  also  reorientational  motion  of  both  components:  polymer 
and  penetrant.  The  molecular  level  information  available  from  NMR  will  provide  insight  into  some 
of  the  concepts  introduced  into  these  theories  which  are  at  least  at  first  glance  based  on  very 
local  descriptions  of  polymer  and  solvent  dynamics. 

The  Fujita  theory  begins  with  the  empirical  result  of  Doolittle(1 5  ■  1 7)  who  showed  that 
the  viscosity  q  of  a  liquid  of  low  molecular  weight  is  given  by: 

lnTi  =  lnA  +  B/f  (1) 

where  A  and  B  are  empirical  constants  for  a  given  liquid  and  f  is  the  fractional  free  volume, 
defined  to  be  equal  to  the  ratio  of  the  volume  of  the  space  not  occupied  by  the  constituent 
molecules,  vf,  over  the  total  volume  of  the  solution,  v: 


(2) 

V 

V  =  VQ  +  Vf 

(3) 

where  vq  is  the  volume  occupied  by  the  molecules  in  accordance  with  Van  der  Waals  radii  and 
vibrational  motions.  Equation  (1)  measures  the  volume  of  a  molecule  relative  to  the  volume  of  a 
void  next  to  it.  When  the  volume  of  the  void  is  larger  than  the  volume  of  the  molecule  the  viscosity 
is  low.  When  the  volume  of  the  void  is  smaller  than  the  volume  of  the  molecule  the  viscosity  is 
high.  Another  assumption  that  is  made  is  that  the  Doolittle  relation  applies  to  each  jumping  unit  in 
polymer  -  solvent  systems(1®).  Jumping  units  are  imagined  to  be  of  the  size  of  small  penetrant 
molecules  or  repeat  units  in  the  case  of  the  polymer.  Also  the  theory  relates  the  mobility  of  the 
polymer  and  the  solvent  jumping  units  to  the  free  volume  per  unit  volume  of  solution. 

In  this  study  the  temperature  and  concentration  dependencies  of  the  self  -  diffusion  of 
toluene  in  PIB  generated  from  the  Pulsed  field  gradient  technique  will  be  examined  and  then  a 
comparison  with  viscosity  data  from  macroscopic  techniques  will  be  attempted.  The  concentration 
dependence  of  the  data  will  be  analyzed  with  the  Fujita  theory(^2)  while  the  temperature 


dependence  will  be  analyzed  using  the  WLF  equation  for  diffusion('^Q),  which  is  also  commonly 
placed  in  the  context  of  free  volume  theories.  The  well  known  form  of  the  WLF  equation  is 


logic 
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Co 


-c?(T-To) 

c^+T-To 


(4) 


where  C  is  the  translational  friction  coefficient  of  a  penetrant  molecule  into  the  polymer  matrix,  Co 
is  the  value  of  C  at  Tq  and  ci  D  and  C2^  are  the  WLF  coefficients  at  the  reference  temperature  Tq. 
The  diffusion  coefficient  ,D  .can  be  introduced  into  this  equation  by  using  the  well  known 
expression : 


D  = 


kT 

C 


(5) 


where  k  is  the  Boltzmann's  constant.  The  WLF  equation  can  thus  be  applied  to  diffusion. 


Experimental 

The  polyisobutylene  sample  used  for  the  study  was  obtained  from  Cellomer  Associates, 
Inc.  of  Webster,  NY  in  the  form  of  thick  sheets  of  high  molecular  weight  (Cat.  #40E,  Lot  #02,  mol. 
wt  1,000,000).  HPLC  Grade  toluene  was  obtained  from  RsherChemical  Fisher  Scientific.  The 
diffusion  measurements  were  performed  on  a  Bruker  MSL  -  300  spectrometer  equipped  with  a 
High  Resolution  Diffusion  Unit  and  gradient  probe,  with  a  proton  frequency  of  300.1  MHz.  For  the 
diffusion  measurements  the  phenyl  proton  peak  of  toluene  was  observed. 

The  well  known  Stejskal  -  Tanner  pulse  sequence(2®)  %  !  2  -  G(8)  -  7t  -  G(5)  -  (echo)  was 
used  to  measure  the  diffusion  coefficient  of  toluene  in  PIB  solutions.  The  stimulated  echo(21) 
pulse  sequence  7t  /  2  -  G(5)  -  n  /  2  -  t  -  7t  /  2  -  G(5)  -  (echo)  was  also  used  for  measuring  the 
diffusion  coefficient,  employing  longer  diffusion  time  observations.  The  strength  of  the  magnetic 
field  gradients  was  varied  from  14  gauss  /  cm  up  to  140  gauss  /  cm;  the  associated  length  5,  from 
500  psec  up  to  7  msec;  and  the  distance  between  the  two  gradient  pulses  A,  from  2  msec  up  to 
80  msec.  Those  values  of  A  correspond  to  diffusion  distances  from  approximately  1 .1  pm  up  to 
9.1  pm. 

Intfirpretation 

The  temperature  dependence  of  toluene  in  PIB  was  analyzed  with  the  following  form  of 
the  WLF  equation("*^): 
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(6) 


where  ci  ^  =  Bj  /  2.303  f(<|)s),  C2^  =  f(<l)s)  /  a((t)s).  f(<t>s)  's  the  fractional  free  volume  of  the  solution 
at  the  reference  temperature  To.  <l)s  is  the  volume  fraction  of  toluene  in  PIB  and  a((|)s)  is  the 
fractional  free  volume  expansion  factor  and  is  the  self  diffusion  coefficient  at  Tq.  Equation  (6) 

can  be  rewritten  as: 


ln( 


DTq 

DtJ 


)  = 


^(T-To) 

%^(%^  +  ^(T-To)) 


(7) 


The  parameters  that  were  floated  in  order  to  fit  the  data  were  a((|>s)  /  Bd  and  f(<t)s)  /  BcJ.  Bd  is 
analogous  to  the  B  in  equation  (1)  and  will  be  defined  precisely  later.  The  fits  along  with  the 
experimental  self  -  diffusion  coefficients  at  different  concentrations  and  temperatures,  are  shown 
in  Figure  1 ,  and  the  reference  temperature  To  was  chosen  to  be  298  K.  The  parameters  obtained 
are  shown  in  Table  1 .  After  a((j)s)/Bd  and  f((|)s)/Bd  were  obtained,  the  WLF  coefficients  ci  ^  and 
C2*^  were  calculated  (from  c-|  ^  =  Bd  /  2.303  f((|)s)  and  C2^  =  f(<t>s) !  «  (^s))  snd  they  are  also 
shown  in  Table  1 . 


To  fit  the  data  at  a  constant  temperature,  as  a  function  of  concentration,  the  Fujita  free 
volume  theory('*2)  was  applied.  As  mentioned  above,  the  theory  begins  with  the  result  of 
Doolittle(‘*6  - 18)  ^^ho  showed  that  the  viscosity  ti  of  ordinary  liquids  of  low  molecular  weight  is 
givervby  equation  (1). 

To  relate  this  to  diffusion,  the  mobility,  m,  of  a  molecule  is  considered  to  be  the  velocity 
with  which  a  molecule  translates  under  the  action  of  unit  force  and  is  inversely  proportional  to  the 
molecular  friction  coefficient  and  the  viscosity  ti  (Stokes  law).  Thus  equation  (1)  can  be  written  as: 

ln(m)  =  C-B/f  (8) 

where  C  is  a  constant.  The  mobility  of  a  molecule  in  a  medium  depends  on  the  probability  that  the 
molecule  will  find  in  its  vicinity  a  hole  which  is  large  enough  to  allow  displacement  into  it.  The 
constant  B  is  considered  to  be  a  measure  of  the  size  of  a  given  hole.  Let  Bd  be  the  value  of  B  for 
the  minimum  hole  which  is  needed  to  allow  a  given  molecule  to  undergo  such  a 
displacement(12),(22)  Therefore  equation  (8)  will  become: 


ln(m)=A-Bd/f 


(9) 


The  definition  of  mobility  in  terms  of  the  self  -  diffusion  coefficient  of  a  penetrant  in  a  polymer 
matrix  is(23); 

D  =  RTm  (10) 

where  R  is  the  gas  constant,  T  is  the  absolute  temperature  and  D  is  the  self  -  diffusion  coefficient. 

In  general  the  free  volume  of  a  given  polymer  -  penetrant  solution  depends  on 
temperature  and  concentration  and  it  is  more  appropriate  to  denote  the  fractional  free  volume  as 
^0",  <t)s)  where  again  <j»s  is  the  volume  fraction  of  the  penetrant.  It  has  been  shown(24)  that  if  the 
increase  in  free  volume  by  the  addition  of  the  diiuent  is  proportional  to  the  volume  of  the  added 
penetrant,  then: 


Vf  =  (vf)2+Y(T)vi  (11) 

V=:V1+V2  (12) 


where  vf  is  the  average  free  volume  of  the  solution,  (vf  )2  is  the  average  free  volume  of  the  pure 
polymer,  y(T)  is  a  proportionality  factor  that  may  be  compared  with  the  fractional  free  volume  of 
the  diluent,  v  is  the  total  volume  of  the  solution,  vi  is  the  volume  of  the  penetrant  and  V2  is  the 
volume  of  the  amorphous  polymer.  By  dividing  equation  (11)  with  equation  (12)  the  fractional  free 
volume  of  the  solution  f(T,  is  found  to  be  a  linear  function  of  the  volume  fraction  of  the 
penetrant  (l>s: 


^0*,  <l)s)  =  ^pCT,  0)  +  ( fs(T)  -  fp(T ,  0) )  <|)s 


(13) 


where  fp(T,  0)  is  the  fractional  free  volume  of  pure  polymer  and  fs(T)  is  a  parameter  that 
represents  the  fractional  free  volume  of  the  solvent  trapped  in  the  polymer  matrix.  fs(T)  -  fp(T,  0) 
represents  the  effectiveness  of  the  penetrant  for  increasing  the  free  volume  when  it  is  dispersed 
in  the  given  polymer.  Equation  (13)  can  be  rewritten  as: 


f(T.(i)s)_fp(T)  f3(T)  fp(T) 

Bd  B,  ^  B,  B, 


(14) 


If  equations  (13)  and  (9)  are  substituted  into  equation  (10)  considering  that  as  <l)s  ->  0,  D  ->  Dq  we 
obtain: 
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Where  <j)s  is  the  penetrant  volume  fraction,  Dq  is  the  diffusion  coefficient  of  toluene  at  the  limit  of 
zero  penetrant  concentration,  fp  is  the  average  fractional  free  volume  of  pure  polymer  and  fs  is 
the  fractional  free  volume  of  toluene  in  the  solution. 

For  each  temperature,  the  diffusion  coefficient  of  toluene  as  a  function  of  concentration 
was  fitted  using  equation  (15).  The  parameters  varied  are  Dq,  (fp  /  Bd)  and  (fg  /  Bd).  The  reason 
for  using  (fp  /  Bd)  and  (fg  /  Bd)  as  fitting  parameters  and  not  fp  and  fg  is  to  avoid  making  any 
assumptions  about  the  parameter  Bd.  The  parameters  (fp  /  Bd)  and  (fg  /  Bd)  do  not  have  physical 
significance  whereas  fp,  fg  and  Bd  do.  The  parameters  produced  by  the  fit  are  shown  in  Table  2 
and  the  lines  corresponding  to  the  fit  in  Figures  2  and  3. 


From  the  diffusion  coefficient  of  toluene  at  the  limit  of  zero  penetrant  concentration  In 
Table  2  and  the  WLF  equation  (6),  WLF  parameters  can  be  obtained  similar  to  those  shown  In 
Table  1 .  Those  are  shown  in  the  Table  3  below. 


Discussion 


A  comparison  can  be  made  between  the  fractional  free  volume  determined  from  the 
temperature  dependence  of  the  translational  diffusion  of  the  toluene  with  that  determined  from  the 
concentration  dependence  of  the  same  quantity.  Such  a  comparison  is  made  in  Figure  4  where 
the  fractional  free  volume  of  the  solution,  f(T,  <|>g),  is  presented  as  a  function  of  the  volume  fraction 
of  toluene  in  PIB.  The  points  come  from  Table  1  and  they  were  obtained  from  the  WLF  equation 
(7)  where  the  temperature  dependence  of  the  toluene  self  -  diffusion  was  analyzed,  except  the 
filled  circle  that  represents  the  fractional  free  volume  of  pure  toluene  (at  293  K)  obtained  from 
viscosity  data(25).  The  solid  line  was  obtained  from  equation  (14)  using  the  values  of  Table  2  for 
fpCO  and  fsCO  (obtained  from  equation  (15)  where  the  concentration  dependence  of  the  toluene 
self  -  diffusion  was  analyzed).  The  error  bars  represent  an  uncertainty  of  approximately  20%  for 
f(T.  <|)s)  in  Table  1,  which  comes  from  quality  of  the  fitting  (namely  the  fitting  of  the  temperature 
dependence  of  the  self  -  diffusion  of  toluene  in  PIB).  As  can  be  seen  from  Figure  4,  the  values  of 
the  fractional  free  volume  of  the  solution  that  are  obtained  from  the  analysis  of  the  temperature 

and  concentration  dependencies  of  the  self  -  diffusion  of  toluene  in  PIB,  are,  in  fact,  close  to  each 
other. 
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The  thermal  expansion  coefficients  of  the  fractional  free  volume  obtained  from  this  work 
can  be  compared  with  typical  literature  values  obtained  from  viscosity  data  as  is  shown  in  Figure 
5.  The  points  come  from  Table  1  and  they  were  obtained  from  equation  (7)  when  the  temperature 
dependence  of  the  toluene  self  •  diffusion  was  analyzed.  The  solid  line  was  obtained  from 
equation  (16): 

tt(<t>s)  =  (1  *  <l>s)  ®p  +  ‘)>s  tts  ^ 

with  the  literature  values('^^'  of  ap  =  4.4  10'^  deg"^  and  as  =  1.710  ^  deg  (at  293  K), 
where  subscripts  p  and  s  refer  to  polymer  and  solvent  respectively.  The  error  bars  represent  an 
uncertainty  of  approximately  20%  for  a(<|)s)  in  Table  1 ,  which  comes  from  quality  of  the  fitting  of 
the  temperature  dependence  of  the  self  -  diffusion  of  toluene  in  PIB.  Reassuringly,  again  the  two 
sets  of  data  are  in  good  agreement  with  each  other. 

A  free  volume  analysis  of  the  viscosity  of  the  toluene  -  PIB  solution  was  done  by 
Fujita(25).  In  that  study  the  fractional  free  volume  of  the  solution,  f(293  K,  (l»s),  and  the  fractional 
free  volume  of  toluene  trapped  in  the  polymer  matrix,  fs(293  K),  are  reported  as  a  function  of 
concentration  at  293  K.  In  reference  (23)  the  author  uses  the  following  equation,  for  the  viscosity 
data  analysis,  with  the  assumption  that  Bp  =  1 : 

f(293K,4..)  (,(293K)  f.(293K)  f.(293K) 

—  B,  B„ 

where  Bp  represents  the  minimum  void  necessary  for  one  flow  unit  to  be  displaced  (i.e.  the 
analog  of  Bd  for  viscous  flow).  The  values  that  are  derived  are(25): 

fp  (293  K)  =  0.048  fs(293  K)  -  fp  (293  K)  =  0.052  (1 8) 

In  order  to  obtain  a  meaningful  comparison  with  this  viscosity  data  it  would  be  desirable 
to  calculate  the  quantities  f(293K,  ({>s).  fp(293  K)  and  fs(293  K)  -  fp(293  K)  from  the  diffusion 
results .  From  Table  2  the  values  of  fs  /  Bd  and  fp  /  Bd  can  be  extrapolated  to  293  K  following  the 
form  f  (T)  =  f(To)  +  a  (T  -  To)  yielding  : 

fs  (293  K)  /  Bd  =  0.263  fp  (293  K)  /  Bd  =  0.084  (1 9) 

Using  equation  (14)  forT  =  293  K  and  the  values  shown  in  equation  (19)  the  parameter  Bd,  which 
is  a  measure  of  the  minimum  void  size  required  for  penetrant  (toluene)  displacement,  is  varied  in 
order  to  match  the  viscosity  data.  The  best  match  of  the  free  volume  curves  versus  volume 
fraction  toluene  is  shown  in  Figure  6.  The  values  derived  from  this  procedure  for  the  minimum 


necessary  void  size  for  diffusion  are  between  0.4  and  0.45.  This  is  a  very  approximate  estimate 
for  Bd  since  an  extensive  viscosity  data  set  as  a  function  of  concentration  that  might  be  compared 
with  the  diffusion  data  presented  here,  was  not  available.  For  Bd  between  0.4  and  0.45  the 
intercept  and  slope  of  the  dotted  line  in  Figure  6  for  diffusion  data  corresponding  to  equation  (14) 
yields ; 


fp  (293  K)  between  0.0336  and  0.0378 
fs  (293  K)  -  fp  (293  K)  between  0.072  and  0.081 

These  values  are  to  be  compared  with  those  given  in  equation  (18).  The  differences  between  the 
values  obtained  from  the  diffusion  data  and  the  viscosity  data  may  be  attributed  to  Bp  which  is 
assumed  to  be  equal  to  1  in  the  viscosity  data  analysis,  the  uncertainties  of  the  free  volume 
parameters,  the  accuracy  of  the  viscosity  and  diffusion  data  and  the  uncertainties  in  the 
knowledge  of  the  ap  and  fg  (the  fractional  free  volume  of  PIB  at  the  glass  transition  temperature) 
forPIB. 


It  would  have  been  better  to  calculate  Bd  from  the  slope  of  the  following  equation(^  2); 

(20) 


B. 


ti(T.<|)s) 


with  Bp  set  equal  1 ,  but  extensive  viscosity  data  as  a  function  of  toluene  concentration  were  not 
available.  If  Bp  is  not  assumed  to  be  equal  to  1 ,  then  the  ratio  Bd  /  Bp  is  estimated. 

From  the  temperature  dependence  of  fp  /  Bd  a  value  for  the  free  volume  expansion 
coefficient  af  can  be  obtained  (given  the  determination  of  Bd  as  between  0.4  and  0.45)  utilizing 
the  following  equation: 

fp(T)  =  fp(T„)+a,(T-T„)  (21) 

The  value  for  af  is  found  to  be  between  3.38x1 0"^  and  3.8x10''^  deg'"*  which  is  close  to  the 
values  reported  elsewhere  from  viscosity  data(^®)  which  range  from  4.4  10‘4  deg'^  to  2.5  10"^ 
deg'"* .  Figure  7  shows  the  temperature  dependence  of  fp  for  Bd  =  0.4  and  for  Bd  =  0.45. 

Figure  8  shows  a  comparison  of  the  free  volume  of  neat  toluene  as  a  function  of 
temperature  obtained  from  viscosity  data(25)  and  the  free  volume  of  toluene  trapped  in  PIB 
obtained  from  the  diffusion  data  presented  here  (Table  2).  It  shows  that  the  fractional  free  volume 


of  toluene  decreases  when  the  latter  is  trapped  in  entangled  PIB  in  agreement  with  the  conclusion 
in  reference  (25). 

In  Figure  9  the  fractional  free  volume  of  the  PIB  component  in  the  solution  (extracted  from 
the  diffusion  data  presented  in  this  study;  Table  2)  is  compared  with  the  free  volume  parameters 
obtained  from  other  studies  of  penetrants  in  polyisobutylene(5  *  '•9).  In  this  case  the  minimum 

void  size  was  assumed  to  be  1 .  The  limited  fractional  free  volume  data  from  other  diffusion 

studies(5  ■  are  closer  to  the  values  determined  in  this  work  than  the  values  determined  from 
viscosity(^  . 

Finally,  looking  at  the  enthalpy  of  activation  for  diffusion,  an  Arrhenius  temperature 
dependence  was  given  to  the  self  -  diffusion  coefficient  of  toluene; 

D  =  Di  exp(-AHD/RT)  (22) 

Plotting  InD  as  a  function  of  1  /  T(K)  the  enthalpy  activation  energy  for  diffusion  is  calculated.  In 
this  procedure  the  non-linearity  which  can  be  seen  in  Rgure  1  and  is  reflective  of  non-Arrhenius 
behavior  is  ignored.  The  values  obtained,  for  different  concentrations,  are  shown  in  Table  4.  The 
value,  13.5  kJ  /  mole,  for  pure  toluene  is  close  to  the  value  (10.9  kJ  /  mole)  reported 

elsewhere(26). 


When  InD  is  not  a  linear  function  of  (1  /  T),  a  temperature  dependent  "apparent" 
activation  energy  can  be  calculated  from  the  WLF  equation,  in  terms  of  the  shift  factor  aj, 
formally  as(®>  ^); 


d(i) 


(23) 


If  Inaj  is  given  by  a  WLF  equation(l®),  then  a  temperature  dependent  apparent  activation  energy 
controlled  by  free  volume  considerations  can  be  obtained(®>  ^  ®): 


2.303RC°C°T 

(C^  +  T-To)^ 


Since  the  WLF  coefficients  are  available  from  the  information  In  Table  1  and  for  T  =  Tq  =  298  K, 
Cl  =  Bd  /  2.303  f((|)s),  C2^  =  f(<l)s)  /  a  (<j)s)  equation  (24  becomes: 
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R(298)2(|l) 

^  <25) 

Bd 

Thus  from  equation  (25)  and  Table  1 ,  we  obtain  Table  6.  A  comparison  of  the  activation  energies 
shown  in  Tables  4  and  5  is  shown  in  Rgure  10. 

Finally  using  equation  (23)  with  Tq  =  298  K  and  the  WLF  coefficients  ciD,  02^,  taken 
from  Table  1 ,  the  activation  energy  of  each  concentration  as  a  function  of  temperature  is 
calculated.  The  results  are  shown  in  Figure  11.  Essentially,  the  activation  energies  determined 
from  either  approach  are  within  experimental  error  of  each  other  though  those  obtained  from 
equation  (25)  more  accurately  reflect  the  trends,  particularly  at  higher  polymer  contents.  It  Is  clear 
from  Figure  1 1  that  the  temperature  dependence  of  the  activation  energy  increases  with  polymer 
content  and  the  apparent  activation  energy  increases  rapidly  with  decreasing  temperature  as 
expected  from  WLF  predictions.  In  contrast  the  behavior  at  high  solvent  content  is  essentially 
temperature  independent  and  more  characteristically  Arrhenius  as  expected. 

Conclusions 

It  is  clear  from  the  above  analysis  that  the  concentration  dependence  of  the  self  - 
diffusion  coefficient  of  toluene  in  PIB  is  well  described  by  the  FujHa  fractional  free  volume  theory. 
Also  it  is  clear  that  the  temperature  dependence  of  the  self  -  diffusion  coefficient  of  toluene  In  PIB 
IS  well  described  by  the  WLF  equation  for  diffusion.  In  addition  it  has  been  established  that  the 
fractional  free  volumes  and  the  fractional  free  volume  thermal  expansion  coefficients  of  the 
solutions  calculated  from  the  temperature  dependence  of  the  self  -  diffusion  coefficient  of  toluene 
in  PIB  (equation  (6))  are  close  to  the  values  obtained  from  the  analysis  of  the  concentration 
dependence  of  the  self  -  diffusion  coefficient  of  toluene  in  PIB  (equation  (7));  a  fact  that  lends 
credibility  to  the  validity  of  the  interpretation. 

Also  the  fractional  free  volumes  of  the  solutions  at  293  K  obtained  from  the  self  -  diffusion 
data  are  close  to  the  values  obtained  from  the  viscosity  data.  This  establishes  the  utility  of  the 
free  volume  approach  in  relation  to  molecular  mobility  of  a  penetrant  in  a  polymer.  Viscosity  is 
macroscopic  measure,  while  the  diffusion  measurements  presented  here  are  microscopic 
measurements  looking  specifically  at  the  penetrant  molecule.  Two  different  experiments  probing 
two  different  levels  of  dynamics  give  approximately  the  same  fractional  free  volume  for  the 
solutions,  when  the  appropriate  value  for  the  minimum  void  size  necessary  for  diffusion  relative  to 
viscosity  was  chosen. 


Another  result  that  has  been  confirmed,  is  the  decrease  of  the  fractional  free  volume  of 
toluene  when  the  latter  is  trapped  in  entangled  PIB.  The  same  conclusion  has  been  reached  from 
viscosity  data(25),  with  a  completely  different  approach.  Also  the  fractional  free  volume 
parameters  for  pure  PIB  obtained  from  this  work  are  close  to  the  values  reported  in  the  literature 
and  extracted  from  other  diffusion  or  viscoelastic  measurements.  This  confirms  the  correctness  of 
the  magnitude  of  the  free  volume  of  the  particular  polymer. 

From  the  preceding  analysis  it  is  obvious  that  Bd,  the  parameter  that  represents  the 
minimum  void  size  necessary  for  a  penetrant  molecule  to  diffuse,  is  important  if  results  obtained 
from  different  methods  that  observe  different  kind  of  motions  are  to  be  compared.  This  is  to  be 
expected  since  the  size  of  a  moving  unit  is  not  the  same  for  different  kind  of  motions.  The 
unfortunate  fact  is  that  it  is  not  possible  to  obtain  the  dimension  of  the  void  when  the  parameter 
Bd  is  known. 
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Table  1 


Parameters  obtained  from  the  temperature  dependence  fit  of  the  toluene  self  -  diffusion 

coefficient,  using  equation  (6). 


<l>s 

f((j)s)/Bd  a(<|)s)/Bd(deg-1 

)  ciD 

cgD  (K) 

0.045 

0.0904 

5.8310-4 

4.803 

155.060 

0.104 

0.1050 

6.8410-4 

4.135 

153.509 

0.205 

0.1270 

7.2510-4 

3.419 

175.172 

0.327 

0.1550 

8.4410-4 

2.801 

183.649 

0.409 

0.1900 

9.6010-4 

2.285 

197.917 

0.712 

0.2700 

1.1010-3 

1.608 

245.455 

1.000 

0.3700 

1.7010-3 

1.174 

217.647 

Table  2 

Parameters 

obtained  from  the  concentration  dependence  fit  of  the  toluene  self  -  diffusion 
coefficient,  using  equation  (15). 

T(K) 

fs  /  B(j 

fp  /  Bd 

Do  (cm3  /  sec) 

298 

0.273 

0.08735 

9.0810-9 

310 

0.300 

0.09780 

2.5510-3 

320 

0.333 

0.10700 

5.1510-3 

329 

0.364 

0.11630 

8.7010-3 

339 

0.381 

0.12500 

1.5510-7 

349 

0.398 

0.13200 

2.3010-7 

358 

0.423 

0.14100 

3.4010-7 

368 

0.449 

0.14450 

4.1010-7 
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Table  3 

Parameters  obtained  from  the  temperature  dependence  fit  of  the  toluene  self  -  diffusion 
coefficient  at  zero  toluene  concentration  (Dq,  shown  in  Table  2),  using  equation  (6). 


<|>s 

^(<l>s)/Bd 

a(<t)s)  /  Bd  (deg-1) 

c-i^ 

C2D(K) 

0.000 

0.0818 

5.3410-4 

5.308 

153.184 

Table  4 

Apparent  activation  energy  for  diffusion  of  toluene  in  PIB  as  a  function  of  concentration,  obtained 

from  the  InD  vs  1  /  T  plot  of  the  self  -  diffusion  coefficient  of  toluene  in  PIB  (Arrhenius  temperature 
dependence). 


<}»s 

AHd  (kJ  /  mole) 

0.000 

50.4 

0.045 

44.4 

0.104 

40.4 

0.205 

34.1 

0.327 

28.6 

0.409 

22.5 

0.712 

12.4 

1.000 

13.5 

16 


Table  5 

Apparent  activation  energy  for  diffusion  of  toluene  in  PIB,  controlled  by  free  volume,  as  a  function 
of  concentration,  obtained  from  equation  (25)  and  Table  1 . 


<t's 

Ed  (kJ  /  mole) 

0.000 

58.9 

0.045 

52.6 

0.104 

-45.8 

0.205 

33.1 

0.327 

25.9 

0.409 

19.6 

0.712 

11.1 

1.000 

9.16 

Figure  Captions 


Figure  1 :  Self  -  diffusion  coefficient  of  toluene  in  PIB  solutions  as  a  function  of  temperature  at 
different  concentrations.  The  solid  lines  represent  fits  to  the  WLF  diffusion  equation 
presented  in  the  text.  The  zero  penetrant  (0%)  concentration  points  are  from  Table  2. 

Figure  2.  Self  -  diffusion  coefficient  of  toluene  in  PIB  solutions  as  a  function  of  concentration  at 
different  temperatures.  Solid  lines  represent  fits  to  Fujita  free  volume  theory 
presented  in  the  text. 

Figure  3:  Self  -  diffusion  coefficient  of  toluene  in  PIB  solutions  as  a  function  of  concentration  at 
different  temperatures.  Solid  lines  represent  fits  to  Fujita  free  volume  theory 
presented  in  the  text. 

Rgure  4:  Comparison  of  the  fractional  free  volume  of  the  toluene  -  PIB  solution,  obtained  from  the 
temperature  (points)  and  concentration  (line)  dependence  of  the  diffusion  coefficient. 
Solid  line  represents  values  obtained  from  equation  (14).  Riled  circle  represents  the 
fractional  free  volume  of  pure  toluene  at  293  k(25). 

Figure  5:  Thermal  expansion  coefficient  of  the  solution  fractional  free  volume  as  a  function  of  the 
volume  fraction  of  toluene  in  polyisobutylene.  Points  represent  values  from  Table  1 
and  the  solid  line  represents  values  obtained  from  equation  (16)  with  ap  =  4.4  10"^ 
deg*"*  and  as  =  1 .7 1 0*3  deg*"* . 

Figure  6;  Comparison  of  fractional  free  volume  of  the  toluene  -  PIB  solution  as  a  function  of 
concentration  obtained  from  diffusion  data  for  Bd  =  0.45  (dashed  line)  with  the  values 
obtained  from  viscosity  data  (solid  line)  from  reference  (25),  for  Bp  =  1.  The  dashed 
line  was  calculated  from  equation  (14)  for  T  =  293  K. 

Figure  7:  Fractional  Free  Volume  of  the  polyisobutylene  component  in  the  solution  as  a  function 
of  Temperature  using  different  values  for  Bd:  Solid  lines  represent  curve  fitting 
according  to  equation  (21). 

Figure  8:  Comparison  of  the  fractional  free  volume  of  neat  toluene(^2)  j|^g  fractional  free 
volume  of  toluene  trapped  in  PIB,  as  a  function  of  temperature. 

Figure  9.  Comparison  of  the  fractional  free  volume  of  the  polyisobutylene  component  in  the 
solution  with  reported  values  in  the  literature  as  a  function  of  temperature. 
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Figure  10:  Activation  energy  for  diffusion  as  a  function  of  the  volume  fraction  of  toluene  in  PIB. 

Open  circles  represent  values  obtained  from  an  Arrhenius  temperature  dependence 
of  the  self  -  diffusion  coefficient  of  toluene  in  PIB.  Riled  circles  were  calculated  from 
equation  (25). 

Figure  1 1 :  Apparent  activation  energy  for  diffusion,  controlled  by  free  volume,  as  a  function  of 
temperature  for  each  concentration,  obtained  from  equation  (23). 
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Abstract 


The  self  -  diffusion  coefficients  of  toluene  in  polyisobutyiene  have  been  analyzed  using 
the  Vrentas  -  Duda  free  volume  diffusion  model.  The  diffusion  coefficients  were  determined  at 
different  temperatures  and  concentrations,  using  the  pulsed  field  gradient  nuclear  magnetic 
resonance  technique.  The  data  were  satisfactorily  described  by  the  model  and  the  size  of  the 
polymer  jumping  unit  was  extracted.  Comparisons  were  made  with  the  Fujita  free  volume  theory 
and  the  Fujita  free  volume  parameters  were  extracted  from  the  Vrentas  -  Duda  free  volume 
parameters.  From  the  diffusion  data  that  now  available,  it  can  be  concluded  that  for  most 
polymers  the  jumping  unit  is  about  1 .5  times  the  polymer  monomer  molecular  weight.  The 
activation  energy  of  the  toluene  diffusion  in  polyisobutyiene  is  compared  with  the  activation 
energies  of  other  penetrants  in  the  same  polymer.  The  diffusion  data  presented  in  this  work  show 
that  the  energy  per  mole  required  to  overcome  the  attractive  forces  which  constrain  a  diffusing 
species  to  its  neighbors  should  be  considered  to  be  zero,  in  order  to  be  able  to  extract  the  free 
volume  parameters  (from  viscosity  and  diffusion  data)  with  an  acceptable  uncertainty. 


Keywords:  Polyisobutyiene,  Toluene,  Pulsed  Field  Gradient  Nuclear  Magnetic  Resonance, 
Diffusion,  Vrentas  -  Duda  Theory,  Fujita  Theory. 


introduction 


In  the  previous  paperd).  the  self  -  diffusion  coefficients  of  toluene  in  polyisobutylene 
(PIB),  as  measured  by  the  pulsed  field  gradient  NMR  technique,  have  been  presented  as  a 
function  of  temperature  and  concentration  and  analysed  using  the  Fujita  free  volume  theory(2).  ln 
that  analysis,  the  free  volume  parameters  were  compared  with  the  ones  obtained  from  the  study 
of  viscosity  data(3)  of  the  same  system  (toluene  -  polyisobutylene).  Also  comparisons  were  made 
with  the  polymer  free  volume  parameters  reported  in  the  literature  which  were  obtained  from 
viscoelastic(^)  and  other  diffusion  measurements. 

In  this  study  the  self  -  diffusion  coefficient  data(1)  of  toluene  in  PIB  are  analysed  using  the 

Vrentas  -  Duda  free  volume  theory(5-1 0).  This  theory  is  known  to  be  able  to  describe  the 

penetrant  self  -  diffusion  coefficient  in  both  rubbery  and  glassy  polymers.  The  theory  introduces  a 

number  of  parameters  which  have  physical  significance  and  furthermore  they  are  not  simply 

adjusted  in  order  to  achieve  a  good  empirical  fit  of  the  diffusion  data.  Most  of  the  parameters  can 

be  determined  independently  from  sources  other  than  diffusion  studies,  such  as  viscosity  and 

C  Ti  data.  Thus  it  may  be  possible  to  predict(11)  general  diffusion  behavior  without  the 

existence  of  any  direct  experimental  data.  However,  quantitative  predictions  of  D  are  not  always 

possible  and  usually  one  to  three  adjustable  parameters  are  needed  to  fit  diffusion  data.  In 

addition  Zielinski  and  Duda(11)  have  pointed  out  the  utility  of  carbon-13  spin-lattice  relaxation 

experiments  in  providing  values  for  the  parameters  in  the  Vrentas-Duda  theory  and  accordingly  a 

detailed  NMR  study  on  pure  PIB  was  completed  prior  to  this  study  of  penetrant  translational 
diffusion. 

The  Vrentas  -  Duda  theory  is  based  on  the  Cohen  -  Turnbull  result(12)  showing  that  the 
mobility  of  a  molecule  in  a  pure  liquid  is  related  to  the  average  free  volume,  vf.  per  molecule 
through  the  following  equation: 

lnm  =  E-yv*/vf 

where  y  is  a  numerical  factor  between  0.5  and  1.  also  known  as  overlap  factor,  that  is  introduced 
because  a  given  free  volume  may  be  available  to  more  than  one  molecule  and  v*  Is  the  minimum 
volume  of  a  hole  that  a  molecule  must  find  in  its  vicinity  in  order  for  it  to  move.  According  to 
Cohen  and  Turnbull  molecules  are  resident  In  cages  bound  by  their  neighbors  for  the  majority  of 
the  time.  When  a  density  fluctuation  occurs  and  a  hole  of  enough  size  opens  within  the  molecular 
cage,  the  molecule  jumps  into  the  hole.  DIffusional  motion  is  considered  to  be  translation  of  a 
molecule  across  the  space  within  its  cage  and  bulk  diffusion  then  occurs  as  a  result  of 
redistribution  of  the  free  volume  within  the  liquid.  The  above  equation  was  derived  under  the 


assumption  that  redistribution  of  free  volume  (holes)  at  constant  volume  requires  no  change  in 
energy.  Equation  (1)  demonstrates  that  when  the  minimum  volume  of  a  hole  that  a  molecule  must 
find  for  a  jump  to  occur  is  larger  than  the  average  free  volume  per  molecule  then  the  mobility  of 
the  penetrant  is  small.  Conversely,  when  v*  is  smaller  than  Vf  then  the  molecule  can  jump  to  a 
neighbouring  position  and  the  mobility  Is  high. 

Vrentas  -  Duda(5)  formulated  their  theory  based  on  the  following  assumptions  as 
stated(®>  ^3). 

(i)  Vf'  is  the  free  volume  equal  to  the  space  not  occupied  by  the  actual  molecules. 

(ii)  The  free  volume  available  to  each  jumping  unit:  solvent  molecule  or  polymer  segment,  is  equal 
to  Vf*  /  N  where  N  is  the  total  number  of  solvent  molecules  plus  polymer  segments  and  Vf"  is 
the  total  free  volume  of  the  solution. 

(iii)  The  critical  hole  volume,  v  ,  of  a  jumping  unit  is  equal  to  its  intrinsic  volume  at  zero  Kelvin. 

The  partial  specific  volumes  of  the  solvent  and  polymer  are  independent  of  composition  i.e. 
there  is  no  volume  change  on  mixing  and  thermal  expansion  coefficients  are  given  by  their 
average  values  over  the  temperature  range  considered. 

(v)  As  the  temperature  is  increased,  the  increase  in  volume  comes  from  the  homogeneous 
expansion  of  the  material  due  to  the  increasing  amplitude  of  the  enharmonic  vibrations  with 
temperature  and  from  the  formation  of  holes  or  vacancies  which  are  distributed 
discontinuously  throughout  the  material  at  any  instant.  The  free  space  that  is  not  occupied  by 
the  constituent  molecules  and  is  not  associated  with  any  holes  or  vacancies  is  called 
interstitial  free  volume.  The  energy  for  redistribution  of  the  interstitial  free  volume  is  very  large 
and  this  type  of  free  volume  therefore  must  be  distributed  uniformly  among  the  molecules. 
The  free  volume  associated  with  the  discontinuous  distribution  of  holes  in  the  liquid  is  called 
the  hole  free  volume.  This  free  volume  can  be  redistributed  with  no  increase  in  energy. 

To  summarize  the  volume  of  a  liquid  can  be  divided  into  two  parts:  the  volume  occupied 
by  the  molecules  of  the  liquid  and  the  empty  space  between  the  molecules,  designated  as  the 
free  volume.  Due  to  thermal  fluctuations,  part  of  the  free  volume  is  continuously  redistributed, 
causing  fluctuations  in  the  local  density.  This  free  volume  is  called  hole  free  volume.  A  molecule 
will  translate  through  these  fluctuating  voids.  According  to  the  theory  in  order  for  a  molecule  to 
undergo  a  diffusing  jump  two  things  must  happen:  (i)  a  hole  of  sufficient  size  must  appear  next  to 
the  molecule  and  (ii)  the  molecule  must  possess  enough  energy  to  jump  Into  the  void.  Thus  the 


diffusion  of  a  molecule  is  controlled  by  the  probability  that  a  hole  of  enough  size  will  form  in  Its 
vicinity  and  by  the  probability  that  the  molecule  will  have  enough  energy  to  jump  into  this  hole.  For 
polymer  -  solvent  systems  above  the  glass  transition  temperature,  the  hole  free  volume  is  small 
and  therefore  the  diffusion  of  the  penetrant  is  mostly  controlled  by  the  probability  of  finding  a  hole 

of  enough  size  in  the  vicinity  of  the  penetrant  molecule  and  by  extension  the  size  of  that  migrating 
molecule. 

In  this  work  the  parameters  of  the  Vrentas  -  Duda  theory  for  the  toluene  -  PIB  system  are 
determined  from  toluene  self  -  diffusion  data  and  compared  with  the  ones  obtained  for  other 
polymers.  Also  the  Fujita  free  volume  theory  Is  known(‘*4)  to  be  identical  to  the  Vrentas  -  Duda 
free  volume  theory  under  certain  conditions,  in  accord  with  expectations.  The  polymer  jumping 
unit  of  PIB  is  compared  with  the  jumping  units  of  other  polymers  as  a  function  of  their  glass 
transition  temperatures,  using  different  polymer  -  solvent  systems.  The  apparent  activation 
energies  for  diffusion  of  different  penetrants  in  PIB  are  also  considered. 

Experimental 


The  materials  used  for  this  study  as  well  as  the  experimental  pulse  sequence  and 
conditions  are  described  in  reference  1. 

Interpretation 


The  self  -  diffusion  of  toluene  In  PIB  was  analysed  via  the  Vrentas  -  Duda  theoryC^'IO). 

According  to  this  theory  the  self  diffusion  coefficient  of  a  penetrant  into  a  rubbery  polymer  is  given 
by(5-10). 


Di  =  Doexp(— )exp[ 


-Kv;+o)p4v;) 


]  (2) 


where  Vj  is  the  specific  critical  hole  free  volume  of  component  i  required  for  a  jump  (i  can  be 
equal  to  1  or  2;  1  for  penetrant  or  solvent  and  2  for  polymer),  coi  is  the  mass  fraction  of  component 
i.  Tgi  is  the  glass  transition  temperature  of  component  i.  Do  is  a  constant  pre-exponential  factor,  E 
is  the  energy  per  mole  that  a  molecule  needs  to  overcome  attractive  forces  which  constrain  it  to 
its  neighbors  and  y  is  an  overlap  factor  (between  1/2  and  1)  which  is  introduced  because  the 
same  free  volume  is  available  to  more  than  one  molecule.  Kn  and  K21  are  free-volume 
parameters  for  the  solvent,  while  K12  and  K22  are  free  volume  parameters  for  the  polymer. 


The  free  volume  parameters  K-]  -i ,  K21 ,  K-j 2i  K22  are  defined  as  follows(^); 


5 


Ki,  =  V?(T,,)[a,-(1-f°)a,J 

K  -  ^H1 _ 

^12  ~  ^2  ("^92  )[®2  ~  (^  ~  ^H2  )®c2  ] 

K  _ _ _ 

^  a2-(1-4)a,, 

fC3  _ 

V?(TgO 


K,,K22 

V^CTga) 


(3) 

(4) 

(5) 

(6) 

(7) 

(8) 


where  aj  Is  the  thermal  expansion  coefficient  for  the  equilibrium  liquid  component  i,  occi  is  the 
thermal  expansion  coefficient  for  the  component  i  for  the  sum  of  the  specific  occupied  volume  and 
the  specific  interstitial  volume,  Vj®  is  the  specific  volume  of  the  pure  equilibrium  component  I  and 
fHi®  is  the  fractional  hole  free  volume  of  the  component  i  at  Tgj. 

Finally,  ^  is  the  ratio  of  molar  volume  of  a  solvent  jumping  unit  to  the  molar  volume  of  a 
polymer  jumping  unit.  The  parameter  ^  can  be  a  measure  of  the  molar  volume  of  the  solvent 
jumping  unit.  When  the  entire  solvent  molecule  performs  a  jump  during  the  process  of  the  self  - 
diffusion,  then  the  molar  volume  of  the  solvent  jumping  unit  is  equal  to  the  molar  volume  of  the 
liquid  solvent.  In  some  cases  the  solvent  jumping  unit  constitutes  only  a  portion  of  the  molecule 
and  hence  the  molar  volume  of  the  jumping  unit  is  smaller  than  the  molar  volume  of  the  liquid 
solvent.  The  solvent  jumping  unit  is  always  smaller  than  or  equal  to  the  molar  volume  of  the  liquid 
solvent.  Also  the  size  of  the  solvent  jumping  unit  that  undergoes  segment  wise  diffusion,  depends 
on  the  characteristics  of  the  polymer  through  which  it  is  diffusing(''8). 

There  are  nine^independent  parameters  in  equation  (2),  Dq,  E.  Kn/y.  K21  -  Tgi, 
Ki2/y,  K22  -  Tg2,  V.J  ,  V2  that  need  to  be  determined.  However,  most  of  them  can  be 

determined  from  viscosity  data  or  sources  other  than  diffusion  data  and  in  favorable  cases  all  of 
them  could  conceivably  be  set  from  other  data. 

The  two  critical  volumes  V.j  ,  and  V2  ,  represent  the  minimum  specific  hole  free  volumes 
required  to  allow  a  solvent  and  polymer  molecule  respectively  to  take  a  diffusive  jump  and  are 
estimated  as  the  specific  volumes  of  the  solvent  and  polymer  at  temperature  T  =  0  K.  Group 


contribution  methods  developed  by  Sugden(''6)  and  Bittz(17)  are  used  for  the  estimation  of  the 

molar  volumes  at  0  K.  Those  methods  were  summarised  by  Haward  in  reference  (18).  For  the 
system  under  study  these  parameters  are(''  I);  *  =  0.917  (cm3  /  gr)  and  Vg*  =  1 .004  (cm3  /  g^). 

The  polymer  free  volume  parameters  Ki2/y,  K22  -  Tg2,  can  be  estimated  either  from 
polymer  viscosity  data  or  from  the  WLF  parameters  through  the  equations(3>  "I); 

^  =  2.303(C?,)(C^)  (9) 

^12 

K22  =  (CI2 )  (l  0) 

where  the  “g"  on  the  exponent  means  that  the  reference  temperature  for  the  WLF  parameters  is 
the  glass  transition  temperature.  NMR  studies  of  the  130  spin  ■  lattice  relaxation  time  as  a 
function  of  temperature  and  Larmor  frequency  can  be  usefuK^^)  for  determining  the  WLF 

coefficients  and  for  PIB  quite  a  few  different  sets  of  WLF  coefficients  have  been  reported  in  the 
literature: 

(1)  From  reference  (4),  from  viscoelastic  measurements  Ci29  =  16.07,  C229  =  107.4  K  and  we 
obtain: 


Ki2/y  =  2.526  10-4  (cm3 /grK),  K22  -  Tg2  =  -  97.6  K. 

(2)  From  reference  (11)  from  viscosity  data  C-\2^  =  16.63,  022^  =  104.4  K  and: 

Ki2/y=  2.51  10-4  (cm3/  grK),  K22  -  Tg2  =  - 100.6  K. 

(3)  From  a  detailed  NMR  study(20)  of  the  ^30  spin  .  lattice  relaxation  times  for  neat  PIB  as  a 
function  of  temperature  and  Larmor  frequency  we  obtain  (for  Ci29  =  9.33,  C229  =  107.4  K): 

Ki 2/y  =  4.351  1 0-4  (cm3  /  gr  K),  K22  -  Tg2  =  ■  97.4  K. 

The  solvent  free  volume  parameters  Kii/y  and  K21  -  Tgi,  can  be  estimated  from 
equations  similar  to  equations  (9)  and  (10)  from  WLF  coefficients.  Viscosity  data  alone  for  the 
solvent  can  provide  all  four  parameters(1 1)  (Kn/y.  K21  -  Tgi ,  Dq  and  E),  however  the  fitting  in 
this  case  is  uncritical  and  cannot  yield  a  unique  set  of  parameter  values.  The  only  way  to 
overcome  this  problem  is  to  assume  E  =  0.  With  this  assumption  a  different  set  of  parameters  is 
obtained  from  those  assuming  E  0,  and  these  parameters  give  better  predictions  for  the 


diffusion  coefficient  of  a  penetrant  into  a  polymer(2^).  The  values  of  the  parameters  that  are 
reported  in  the  literature  are: 

(1)  From  reference  (19),  for  toluene  from  viscosity  data,  C-|  i9  =  10.75,  C2i9  =  23.6  K  which  give 
(Tgi=117K): 

K1  i/y  =  1 .569  1 0-3  (cm3  /  gr  k),  K21  -  Tgi  =  -  93.4  K. 

(2)  From  reference  (19),  for  toluene  from  NMR  data,  CuS  =  6.64,  C2l9  =  34.51  K  which  give 
(Tgi  =117K): 

Ki  i/y  =  1 .738  1 0*3  (cm3  /  gr  K),  K21  -  Tgi  =  -  82.49  K. 

(3)  From  reference  (22),  from  viscosity  data  for  toluene: 

Ki  i/y=  2.21  10-3  (cm3  /  gr  k)_  K21  -  Tgi  =  - 103  K. 

(4)  From  reference  (9),  from  viscosity  data,  for  toluene: 

Ki  i/y  =  1 .57  1 0-3  (cm3  /  gr  k),  K21  ■  Tgi  =  -  90.5  K. 

(5)  From  reference  (23)  from  viscosity  data  for  toluene,  Ci  1 9  =  1 5.8,  C21 9  =  1 1 .4  K  which  give 

(Tgi  =114  K): 

Ki  i/y = 2.21 1  10-3  (cm3 /grK),  K21 -Tgi  =-102.6  K. 

(6)  From  reference  (11),  from  viscosity  data: 

Ki  i/y  =  1 .45  10-3  (cm3  /  gr  k),  K21  -  Tgi  =  -  86.32  K. 

The  differences  in  the  free  volume  parameters  are  due  to  experimental  errors  in  the 
viscosity  data  and  to  the  uncertainty  in  the  fitting  parameters. 

The  pre-exponential  factor  Do,  is  considered  to  be  solvent  dependent(^  In  order  to 
verify  if  this  is  true,  the  Dq  value  obtained  from  the  analysis  of  the  viscosity  data  of  pure  toluene  is 
used  to  describe  the  diffusion  of  toluene  in  different  polymers.  It  has  been  foundC’*'*)  that  a 
specific  value  of  Dq  adequately  describes  the  transport  of  toluene  in  various  polymers 
(polystyrene,  polyvinyl  acetate,  polymethyl  methacrylate).  A  recent  report(^^)  shows  that  the 
value  of  Do  associated  with  the  free  volume  parameters  Ki  i/y  =  1 .45  10-3  (cm3  /  gr  K)  and  K21  - 
Tgi  =  -  86.32  K,  is  4.82  10-^  cm^  /  sec.  It  must  be  noted  that  if  a  different  set  of  values  for  Kn/y 
and  K21  -  Tgi  are  used,  then  a  different  number  for  Do  should  be  used  also. 


The  last  parameter  left  to  be  determined,  is  the  ratio  of  the  molar  volume  of  a  solvent 
jumping  unit  to  the  molar  volume  of  a  polymer  jumping  unit.  For  the  case  of  small  penetrants  / 
solvents  such  as  benzene,  methane  etc.  that  move  as  single  units,  ^  is  defined  as  : 

(11) 

where  Vi  0(o)  fe  the  solvent  molar  volume,  Vgj  is  the  polymer  jumping  unit  volume.  Mi  is  the 

molecular  weight  of  the  entire  solvent  and  M2j  is  the  molecular  weight  of  the  polymer  jumping 

unit.  It  has  been  shown(24)  that  (y  /  Ki  2)  depends  linearly  on  the  solvent  molar  volume  (at  0 

K)  and  the  resulting  straight  line  passes  through  the  origin  when  the  solvent  moves  as  a  single 
unit: 


=  PV?(0) 


(12) 


Equation  (12)  suggests  (using  the  definition  of  ^  above)  that  the  size  of  the  polymer  jumping  unit 
doesn't  depend  on  the  solvent  and  is  polymer  dependent.  Therefore,  if  p  Is  known  for  a  particular 
polymer,  the  value  of  ^  for  any  solvent  in  that  polymer  can  be  calculated  accurately,  if  the  solvent 
diffuses  as  a  single  unit.  In  order  to  estimate  p  for  a  polymer,  a  lot  of  diffusion  data  of  different 
penetrants  must  be  studied  first,  in  order  to  extract  the  values  of  ^  and  then  make  the  plot  of  (y 
V2  §  /  K12)  vs  Vi  0(0)  (the  slope  will  be  p).  After  p  is  known,  predictions  of  the  diffusion  coefficient 
of  a  penetrant  in  the  particular  polymer  can  be  made.  For  polyisobutylene  the  value  of  p  is  not 

known  and  thus  the  parameter  ^  will  be  the  only  parameter  that  will  be  floated  in  order  to  calculate 
the  diffusion  coefficient  of  toluene  in  polyisobutylene. 

In  order  to  fit  the  diffusion  data  of  toluene  in  PIB  as  a  function  of  concentration,  the 
parameters  shown  in  Table  1  were  used.  Since  all  the  toluene  and  polymer  free  volume 
parameters  mentioned  above  are  close  to  each  other,  the  actual  choice  of  a  parameter  set  is 
somewhat  arbitrary  and  doesn't  significantly  change  the  quality  of  the  fit. 

As  a  result  of  the  fit,  ^  was  found  to  be  equal  to  0.892  ±  0.08.  The  experimental  self  - 
diffusion  coefficients  of  toluene  in  PIB  along  with  the  fitting  curves  (solid  lines)  are  shown  as  a 
function  of  temperature  and  concentration  in  Rgures  1  (a)  and  1  (b). 

Discussion 


The  filling  of  the  concentration  dependence  ol  the  self  -  diffusion  of  toluene  In  PIB  has 
deterntined  a  value  for  i  which  is  0.892  ±0.08.  This  value  corresponds  to  a  polymer  iumping  unit 
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of  molecular  weight  94.3  as  derived  from  equation  (11)  when  the  toluene  rrxjves  as  a  single  unit. 
This  value  is  close  to  the  molecular  weight  of  toluene  (M-|  =  92.13).  It  is  thus  seen  that  the  size  of 
the  polymer  jumping  unit  is  approximately  the  same  as  the  size  of  the  penetrant  jumping  unit-.  In 
view  of  this  result,  it  would  be  appropriate  to  expect  that  both  polymer  and  penetrant  jumping 
units  would  require  approximately  the  same  minimum  hole  size  necessary  for  a  diffusion  jump.  It 
has  been  shown(5)  previously  that  both  Fujita  and  Vrentas  -  Duda  free  volume  theories  can  be 
used  when  the  molecular  weights  of  the  two  jumping  units  are  the  same.  As  indicated  in  reference 
(14),  the  two  theories  are  identical  when  the  following  condition  is  true: 


V°  V? 


(13) 


and  when  Bd  (minimum  hole  size  required  for  a  given  molecule  to  permit  a  displacement, 
according  to  Fujita  theory(2))  is  defined  as: 


Bd 


(14) 


This  can  be  done  by  rewriting  equation  (1)  as: 


v;f2, 


ln(  D,  )-  " 

D.(0)  [l+,|,i( 


V?Y' 


f 


(15) 


where  Di  (0)  is  the  self  -  diffusion  coefficient  of  the  penetrant  at  the  limit  cai  =  0,  Vp  is  the  specific 
volume  of  pure  component  i  at  the  temperature  of  interest,  Vj*  is  the  specific  critical  hole  free 
volume  of  component  i  required  for  a  jump,  <(11  is  the  volume  fraction  of  the  solvent  and  fj  is  the 
fractional  hole  -  free  volume  of  pure  component  i  at  the  temperature  of  interest  and  the  fj's  are 
given  by: 


(K,,  +  T-T,,)K,, 
V? 


(K22  +  T-T,2)K,2 


(16) 


Using  equations  (13)  and  (14),  equation  (15)  can  be  written  as: 
Inf— )  =  B()(fi  —  f2  Hi 

^D,(0)^  l  +  faCfi-faHi 


(17) 


Equation  (1 7)  is  identical  to  the  equation  derived  by  Fujita(2).  For  ^  =  0.892,  V2*  =  1 .004  cm^  /  gr. 
Vi  =  0.917  cm3  /  gr,  V2O  =  1.091  cm^  /  gr,  and  Vi^  =  1.154  cm^  /  gr  we  find  that  the  left  hand 
side  of  equation  (13)  is  equal  to  0.821  and  the  right  hand  side  of  equation  (13)  is  equal  to  0.794. 
Therefore  we  see  that  since  the  two  ratios  are  close,  either  of  the  two  theories  can  be  used  for  the 
toluene  -  polyisobutylene  system.  Thus  It  is  no  surprise  that  both  approaches  accurately  describe 
the  temperature  and  concentration  dependence  of  the  self  -  diffusion  constant. 


Also  as  it  has  been  pointed  out(^^)  the  free  volume  parameters  of  the  Fujita  theory  can 
be  determined  from  the  free  volume  parameters  of  the  Vrentas  -  Duda  theory  if  they  are  known. 
This  Is  done  with  the  help  of  the  following  equations: 

B,  =  B,  f7 


where  fi  and  f2  are  given  by  equation  (16)  and 

v;fp 

V°Y  V?Y 


(19) 


Vj3  is  the  specific  volume  of  pure  component  i,  fp(T)  is  the  average  fractional  free  volume  of  pure 
polymer  (as  defined  in  Fujita  theory(2)),  fs(T)  the  fractional  free  volume  of  solvent  (toluene)  in  the 
solution  at  the  limit  of  zero  toluene  concentration  (as  defined  by  Fujita  for  the  toluene  -  PIB 
system(3))  and  Bj  a  constant  that  represents  the  minimum  hole  size  required  for  a  given 
molecule  to  permit  a  displacement,  (as  defined  in  Fujita  theoty(2)). 


Using  equations  (16),  (18)  and  (19)  and  the  parameters  from  Table  1,  the  Fujita 
parameters  can  be  predicted.  Using  T  =  298  K,  and  the  free  volume  parameters  (Table  1)  we 
obtain: 


fp(298K) 


=  0.096; 


fs(298K) 

Bp 


0.32  . 


(20) 


These  values  are  close  to  the  ones  (for  T  =  298  K,  fp(T)  /  Bd  =  0.08735  and  fs(T)  /  Bd  =  0.273 ) 
obtained  from  the  actual  fit  using  the  Fujita  theory  (Paper  I  of  this  series).  This  verifies  the 
correspondence  of  the  two  theories  when  the  conditions  described  by  equations  (13)  and  (14)  are 
true  and  demonstrates  that  it  is  possible  to  determine  the  Fujita  free  volume  parameters  if  the 
Vrentas  -  Duda  free  volume  parameters  are  known. 


The  results  obtained  here  for  ^  from  the  PIB  -  toluene  system  can  be  compared  with  the 
data  on  a  number  of  other  polymers.  If  equation  (11)  is  multiplied  and  divided  by  M2mono.  the 
molecular  weight  of  the  polymer  monomer  unit,  it  can  be  rewritten  as: 

V?(0) 

Ma, 

In  Figure  2  a  plot  of  ^  vs  V-|  ®(0)  /  V2  M2mono  is  shown  for  different  polymers.  The  values  of 
Vi  0(0)  /  V2*  M2mono  were  calculated  from  Tables  2  and  3  from  reference  (1 1)  where  the  Vi  0(0), 
^'t2mono>  M2j.  V2  are  given  for  different  polymers  and  solvents.  For  each  polymer  shown  in 
Rgure  2,  the  parameters  M2monoi  ^2)  and  V2  are  given  in  Table  3  in  reference  (11).  Then  from 
Table  2  (in  reference  (1 1))  the  values  of  Vi  0(o)  are  taken  for  different  solvents.  Those  solvents 
are:  acetone,  benzene,  carbon  tetrachloride,  chloroform,  ethyl  benzene,  n-hexane,  methyl 
acetate,  methylene  chloride,  n-octane,  n-pentane,  tetrahydrofuran,  toluene,  o-xylene  and 
cyclohexane.  Thus  for  every  polymer  the  value  of  ^  is  calculated  from  equation  (12),  for  every 
solvent  mentioned  above,  since  the  values  of  p,  for  some  polymers,  are  known  and  reported  in 
Table  3  in  reference  11.  The  values  of  ^  for  polybutadiene  were  calculated  using  Table  2  from 
reference  (11)  and  the  free  volume  parameters  reported  in  reference  (25).  The  values  of  ^  for 
polyisobutylene  were  calculated  using  Table  2  from  reference  (2)  and  the  free  volume  parameters 
showed  in  the  Table  1  above.  The  values  reported  here  for  PIB  and  polybutadiene  are  estimates 
from  the  currently  available  data  obtained  from  a  single  polymer  -  penetrant  system;  and  thus  the 
uncertainty  is  large.  Better  estimates  for  p  and  M2j  could  be  made  when  additional  diffusion 
data  are  available. 

The  parameter  p  is  accurately  known  for  polystyrene  (PS),  poly(methyl  acrylate)  (PMA), 
poly(methyl  methacrylate)  (PMMA),  poly(ethyl  methacrylate)  (PEMA)  and  poly(vinyl  acetate) 
(PVAc)  because  of  extensive  prior  studies(1 "'). 

Figure  2,  plotted  according  to  equation  (21),  shows  that  most  of  the  data  have  a  slope 
about  M2mono  !  ^2\  =  2/3  which  indicates  that  the  size  of  the  polymer  jumping  unit  is  1 .5  times 
as  large  as  the  monomer  segment.  PMMA  seems  to  be  an  exception;  its  jumping  unit  is  almost 
two  times  larger  than  the  monomer  unit.  For  polyisobutylene  M2j  /  M2mono  is  1 .68  which  is  close 
to  the  typical  value  of  1 .5.  From  equation  (1 1)  it  can  also  be  concluded  that  for  the  same  solvent, 
as  the  polymer  jumping  unit  becomes  larger,  the  value  of  %  becomes  smaller.  This  implies  faster 
diffusion  for  a  larger  polymer  jumping  unit  which  seems  inconsistent  though  the  complex 
dependency  of  D  on  a  number  of  parameters  given  by  equation  (2)  makes  it  difficult  to  predict 
trends  for  a  single  parameter. 


In  Figure  3  the  molecular  weight  of  the  polymer  jumping  unit  is  plotted  as  a  function  of  the 
polymer  monomer  molecular  weight  .  The  solid  line  has  a  slope  of  1.5.  This  plot  shows  more 
clearly  that  the  polymer  jumping  unit  is  approximately  1.5  times  polymer  monomer  molecular 
weight.  Polyisobutylene  is  consistent  with  the  other  polymers  in  this  plot.  Of  course  more  data 
must  be  collected  for  more  polymers  in  different  solvents  in  order  to  be  able  to  conclude  that  the 
molecular  weight  of  the  polymer  jumping  unit  is  in  general  1 .5  times  the  molecular  weight  of  the 
polymer  monomer  unh. 


To  try  to  relate  the  size  of  the  polymer  jumping  unit  with  the  chain  stiffness('*  1 ),  a  plot  of 
the  polymer  glass  transition  temperature  vs  the  molar  volume  of  the  polymer  jumping  unit  is  made 
(Rgure  4).  The  data  are  derived  from  equation  (11)  using  the  parameters  given  in  Tables  2  and  3 
in  reference  (11),  while  toluene  was  chosen  arbitrarily  as  the  penetrant.  Again  the  data  for  the 
polybutadiene  were  taken  from  reference  (25).  The  validity  of  the  relationship  of  V2j  with  the 

polymer  glass  transition  temperature  reported  in  reference  (1 1) 

V2j  (cm3  /  mole)  =  0.6224  Tg2  (K)  -  86.95  (22) 

is  not  evidenced  when  two  new  data  points  were  added;  that  of  PIB  and  that  of  polybutadiene 
(parameters  taken  from  reference  (25)).  Therefore  it  seems  that  additional  diffusion  data  are 
required  in  order  to  have  a  better  understanding  about  the  relationship  between  the  volume  of  the 
polymer  jumping  unit  and  the  polymer  glass  transition  temperature. 


A  constant  temperature  independent  activation  energy  for  diffusion  can  be  defined 
relative  to  the  Arrhenius  equation: 


D  =  Doexp(-E*/RT)  (23) 

by  taking  the  derivative  of  InD  with  respect  to  (1  /  T)  and  multiplying  by  (-R).  That  is,  in  general 
one  can  define  24). 


Eo=-R^  =  RT2^'"^ 


3(:p) 


dT 


(24) 


where  when  InD  varies  linearly  with  (1  /  T),  then  Ed  =  E*.  However  when  InD  does  not  vary 

linearly  with  (1  /  T),  then  the  last  equation  formally  defines  a  temperature  dependent  “activation 
energy". 


Equation  (2)  can  be  rewritten  at  the  limit  of  zero  penetrant  concentration  (ooi  =  0  and  02 
=  1)  and  for  E  =  0: 


lnD  =  lnDo 


K,2(K22-T^+T) 


(25) 


Therefore  the  apparent  activation  energy  for  diffusion  in  the  limit  of  zero  penetrant 
concentration,  which  is  obtained  by  equations  (24)  and  (25),  is  related  to  the  free  volume 
parameters  by  the  expression^'^  ’  ^4, 26). 


RT^(^) 

_ ^2 

(K^  +  T-Tg^r 


(26) 


Using  the  parameters  shown  in  Tabie  1 ,  the  apparent  activation  energy  for  diffusion  at  T  =  298  K 
is  estimated  to  be  about  37.8  kJ  /  mole.  The  activation  energy  as  a  function  of  temperature  is 
shown  in  Figure  5.  The  values  obtained  are  close  to  the  values  reported  in  reference  (20)  for 
segmental  motion  of  PIB,  but  lower  than  the  values  obtained  from  the  Fujita  analysis  of  penetrant 
diffusion(l)  (paper  1  in  this  series). 

From  equation  (26),  if  the  apparent  activation  energy  for  diffusion  is  known,  then  the 
value  of  the  expression  ('ifV2*^/  K12)  can  be  extracted  and  then  ^  obtained.  From  the  Fujita 
interpretation,  it  was  found  that  the  apparent  activation  energy  at  the  limit  of  zero  toluene 
concentration  is  58.9  kJ  /  mole  at  298  K.  Using  this  determination,  the  value  of  (7V2*§  /  K12)  is 
found  to  be  about  3204  K.  This  value  is  compared  with  other  values  obtained  for  other  penetrants 
in  polyisobutylene(26)  in  Rgure  6.  The  values  for  /  K12)  and  the  corresponding  for  the 
different  polymer  -  penetrant  systems,  are  shown  in  Table  2.  The  data  were  taken  from  reference 
(26),  except  for  the  point  that  corresponds  to  the  toluene  -  PIB  system.  The  ^ 's  were  extracted 
using  V2*  =  1 .004  (cm^  /  gr)  and  Ki  2/y  =  2.51  1 0*4  cm^  /  gr  K  (which  was  obtained  from  viscosity 
data)  for  all  the  penetrant  -  polymer  systems.  For  the  toluene  -  PIB  system,  in  order  to  be 
consistent  with  the  rest  of  the  data  presented  in  the  graph,  the  value  of  Ki2/y = 2.51  1 0*4  cm^  /  gr 
K  was  used  in  the  calculation  of  (yV2*^/  K-|2)  instead  of  K-|2^  =  4.351  10‘4  cm^/  gr  obtained 
from  the  NMR  data.  If  Ki  2/y  =  4.351  1 0*4  cm^  /  gr  is  considered  then  ^  is  found  to  be  1 .38  for  the 
toluene  -  PIB  system.  This  last  estimate  of  ^  is  very  approximate  and  has  a  large  error  because  of 
the  error  in  the  estimation  of  the  activation  energy  and  the  uncertainty  in  the  Ki2/y  parameter. 
The  value  ^  =  0.892  is  more  accurate  since  it  was  extracted  directly  from  the  diffusion  data  and 
implies  that  toluene  is  moving  as  a  single  unit  in  accord  with  recent  conclusions(^  ^  •  ^4, 19,  21) 
From  Figure  6  it  also  can  be  seen  that  the  eight  solvents  (filled  circles)  that  may  be  expected  to 


move  as  single  units  deviate  from  a  straight  line  in  contradiction  to  equation  (12).  This  deviation 
may  be  due  to  experimental  error  since  primitive  diffusion  data  from  different  laboratories  are 
being  used.  Another  possibility  is  that  the  energy,  E  per  mole,  that  a  molecule  needs  in  order  to 
overcome  the  attractive  forces  which  constrain  it  to  its  neighbors  makes  a  significant  contribution 
to  the  apparent  activation  energy,  Eq,  for  diffusion  and  therefore  equation  (26)  should  be 
modified  as(2®); 


Rr(^) 


Ed=E  + 


K 


12 


(27) 


Assuming  that  E  =  6000  cal  /  g  mole,  a  recalculation  of  {fJ2%  I K12)  trom  equation  (27)  shows 
that  the  points  under  question  can  be  reasonably  represented  by  a  straight  line  according  to 
equation  (12),  as  it  is  demonstrated  in  reference  (26). 

However,  the  assumption  that  E  is  different  from  zero,  introduces  one  more  parameter 
into  the  fitting  process.  This  would  increase  the  uncertainty  of  the  fitting  parameters  by  a 
significant  percentage.  Also,  this  would  have  an  effect  on  the  values  of  the  parameters  Do,  Ki  1  / 
y  and  K21  -  Tg-|  which  are  obtained  from  the  analysis  of  viscosity  and  specific  volume  data  for  the 
pure  solvent.  If  E  is  included,  a  four  parameter  regression  (Do,  E,  Ki  1  /  y  and  K21  -  Tgi)  would 
be  required  to  fit  the  temperature  dependence  of  the  viscosity  and  specific  volume  data(21).  This 
would  make  it  very  difficult  to  produce  fits  which  converge  to  a  single  set  of  parameters  describing 
the  viscosity  and  specific  volume  data(2'l).  Thus  the  assumption  that  E  is  different  from  zero 
would  increase  the  uncertainty  in  Do,  Ki  t  /  y  and  K21  •  Tgi  .  Also  it  has  been  demonstrated(21) 
that  the  Vrentas  -  Duda  theory  is  by  far  more  successful  in  predicting  the  self  -  diffusion  coefficient 
of  the  penetrant  in  different  polymer  -  solvent  systems  (toluene  -  polystyrene;  ehtylbenzene  - 
polystyrene:  benzene  -  polystyrene)  when  the  assumption  of  E  =  0  is  employed.  Thus  it  is  our 
feeling  that  the  deviation  from  a  straight  line  of  the  eight  points  (filled  circles  in  Figure  6)  is  due  to 
experimental  error  and  not  due  to  the  neglect  of  E  in  equation  (27). 

Another  characteristic  that  is  seen  in  Rgure  6  is  the  existence  of  a  ceiling  value  of  (yV2*^ 

/  Ki  2)  and  therefore  of  Eq  for  large  penetrant  molecules.  This  seems  to  be  in  contradiction  with 
the  free  volume  theory  as  expressed  by  equation  (26).  For  a  given  polymer,  if  the  size  of  the 
penetrant  increases,  %  increases  and  thus  Ed  increases.  Thus  free  volume  theory  predicts  that 
Ed  should  increase  indefinitely  as  the  size  of  the  penetrant  molecule  increases.  However  the 
apparent  activation  energy  Ed  does  not  necessarily  increase  with  an  increase  in  the  size  of  an 
entire  penetrant  molecule  if  the  penetrant  molecule  moves  in  a  segment  wise  manner,  i.e.  if  the 


penetrant  jumping  unit  is  a  portion  of  the  molecular  unit.  Therefore  the  ceiling  value  of  (7V2  ^  / 
K12)  in  Figure  6  for  the  penetrants  n-dodecane,  n-hexadecane,  1,1-diphenylethane  and  2,2,5- 
trimethyl  hexane  may  be  an  indication  of  segment  wise  motion  for  these  penetrants.  This  would 
not  be  very  surprising  because  of  the  size  and  structure  of  these  molecules. 

Conclusions 

The  Vrentas-Duda  free  volume  interpretation  of  the  diffusion  of  toluene  in  PIB  is  of 
comparable  quality  to  the  Fujita  interpretation  presented  in  the  preceding  paper.  Indeed,  this 
should  be  so  since  the  two  theories  coincide  when  the  size  of  the  jumping  unit  of  the  polymer 
matches  the  size  of  the  penetrant  as  it  does  in  this  system.  It  is  informative  to  consider  the 
parameters  of  the  Vrentas-Duda  theory  for  this  system  relative  to  other  systems  as  was 
presented  in  various  ways  in  the  Discussion  section.  Nevertheless,  both  free  volume 
interpretations  seem  equally  appropriate  and  informative  for  this  system. 

It  would  be  desirable  if  the  ideas  from  free  volume  concepts  were  able  to  shed  some  light 
on  the  rather  unique  behaviour  of  PIB  relative  to  other  polymers.  In  particular  the  unusual  slow 
diffusion  exhibited  by  PIB  with  respect  to  a  variety  of  penetrants  might  be  expected  to  be  reflected 
in  the  value  of  certain  free  volume  parameters.  This  is  seen  in  the  fact  that  the  WLF  parameters: 
ci  and  C2,  are  considerably  different  from  most  other  polymers,  which  can  be  adequately 
represented  by  the  so-called  universal  values(^).  The  parameter  a,  the  fractional  free  volume 
expansion  factor,  is  considerably  smaller  than  is  typical  for  other  polymers  and  does  give  some 
physical  insight  into  the  origin  of  the  lower  free  volume  at  temperatures  above  the  glass  transition. 
This  parameter  is  derived  directly  from  the  Fujita  analysis('l).  The  Vrentas-Duda  parameters, 
though  useful  in  their  potential  for  predictability,  do  not  appear  to  yield  any  physical  insight  into 
uniqueness  of  PIB. 
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Table  1 

Values  of  parameters  for  PIB  /  Toluene  systems. 


Vi*(cm3/gr) 

0.917 

V2*  (cm3  /  gr) 

1.004 

K-|-|  /y(cm3/grK) 

1.45  10-3 

Ki2/Y(cm3/grK) 

4.351  10-4 

K21-Tgi(K) 

-86.32 

K22-Tg2(K) 

-97.6 

0.892  ±  0.08 

Do  (cm2  /  sec) 

4.82  10-4 

E  (kcal  /  mole) 

0 

Ml 

92.13 

V-fO  (cm3/gr  mole) 

84.48 

Tgl(K) 

117 

M2mono 

56.04 

0129 

9.33 

0229  (K) 

107.4 

Tg2(K) 

205 

Table  2 


Parameter  ^  obtained  from  equation  (26)  for  different  penetrants  in  PIB. 


Penetrant 

7V2*^/Ki2(K) 

1 .  Hydrogen 

1730 

0.433 

2.  Oxygen 

2460 

0.616 

3.  Nitrogen 

2660 

0.667 

4.  Carbon  dioxide 

2660 

0.667 

5.  n  -  Butane 

3820 

0.958 

6.  iso  -  Butane 

4010 

1.00 

7.  n  -  Pentane 

3660 

0.918 

8.  iso  -  Pentane 

4150 

1.04 

9.  Neopentane 

4120 

1.03 

10.  n  -  Dodecane 

3620 

0.907 

11.  n  -  Hexadecane 

3620 

0.908 

12. 1,1-Diphenylethane 

3780 

0.946 

13.  Cyclohexane 

3580 

0.897 

14.  2,2,5-Trimethylhexane 

3910 

0.981 

15.  Toluene 

3204 

0.801 

figure  Captions 


figure  1(a)  and  1(b):  Diffusion  coefficient  of  toluene  in  PIB  as  a  function  of  concentration  at 

different  temperatures.  Solid  lines  represent  diffusion  coefficients 
calculated  from  equation  (2). . 

Figure  2;  A  comparison  of  M2j  and  M2mono  as  demonstrated  by  equation  (20)  for  various 
polymers:  polystyrene  (PS),  poly(methyl  acrylate)  (PMA),  poly(methyl  methacrylate) 
(PMMA),  poly(ethyl  methacrylate)  (PEMA)  and  poly(vinyl  acetate)  (PVAc), 
polyisobutylene  (PIB),  polybutadiene  (PB).  The  solid  line  indicates  that  the  molecular 
weight  of  the  polymer  jumping  unit  is  1 .5  times  as  large  as  the  molecular  weight  of  the 
monomer  unit.  Various  points  of  the  same  marker  mean  different  solvents.  Those 
solvents  are:  acetone,  benzene,  carbon  tetrachloride,  chloroform,  ethylbenzene,  n- 
hexane,  methyl  acetate,  methylene  chloride,  n-octane,  n-pentane,  tetrahydrofuran, 
toluene,  o-xylene  and  cyclohexane. 

Figure  3:  The  molecular  weight  of  the  polymer  jumping  unit  M2j  as  a  function  of  M2mono  the 
molecular  weight  of  the  polymer  monomer  unit,  polystyrene  (PS),  poly(methyl 
acrylate)  (PMA),  poly(methyl  methacrylate)  (PMMA),  poly(ethyl  methacrylate)  (PEM) 

and  poly(vinyl  acetate)  (PVAc),  polyisobutylene  (PIB),  polybutadiene  (PB).  The  solid 
line  has  a  slope  1.5. 

Figure  4:  Polymer  glass  transition  temperature  vs  the  molar  volume  of  the  polymer  jumping  unit; 

polystyrene  (PS),  poly(methyl  acrylate)  (PMA),  poly(methyl  methacrylate)  (PMMA), 
and  poly(vinyl  acetate)  (PVAc),  polyisobutylene  (PIB),  polybutadiene  (PB).  cis  - 1,  4  - 
poly(isoprene)  (PI),  poly(styrene  -  butadiene)  (SB). 

Figure  5:  The  activation  energy  for  diffusion  for  the  PIB  -  toluene  system  as  a  function  of 
temperature  in  the  limit  of  zero  penetrant  concentration. 

Figure  6:  (yV2*§  /  K12)  vs  Molar  volume  of  solvent  at  0  K  for  polyisobutylene  solvent  systems. 

Numbers  designating  solvents  correspond  to  those  in  Table  3.  Solid  circles  represent 
solvents  which  might  be  expected  to  move  as  single  units  and  open  circles  represent 
solvents  that  may  indicate  a  segmentwise  motion.  Data  and  graph  were  taken  from 

reference  (20).  The  rhombus  represents  toluene  using  parameters  were  taken  from 
this  work. 


D  (cm^/sec) 


.  22 


A  Nuclear  Magnetic  Resonance  Study  of  Toluene  >  Polyisobutylene  Solutions.  (Ill)  Spin- 
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Abstract 


The  spin  -  lattice  relaxation  times  are  determined  for  the  methylene  carbon  of 
polyisobutylene  (PIB),  as  well  as  for  the  ortho  carbon  of  toluene  in  toluene  -  polyisobutylene 
solutions.  The  Hall  -  Helfand  correlation  function  combined  with  restricted  anisotropic  rotational 
diffusion  was  used  to  treat  the  T-j  data  of  the  methylene  carbon  of  PIB.  A  simple  exponential 
correlation  function  was  used  to  describe  the  local  motion  of  toluene  in  the  solutions  which  falls  in 
the  extreme  narrowing  limit  for  the  solutions  studied.  Both  models  described  satisfactorily  the 
temperature  and  field  dependence  of  the  spin  -  lattice  relation  times.  From  the  temperature 
dependence  of  the  correlation  times  for  the  polymer  segmental  motion,  the  free  volume  of  the 
solution  at  each  concentration  is  extracted  and  compared  with  the  values  obtained  from  previous 
studies  of  the  translational  motion  of  the  toluene  penetrant.  The  free  volume  values  extracted 
from  the  T 1  data  for  the  methylene  carbon  of  PIB  and  the  self  -  diffusion  data  for  the  toluene  were 
found  to  be  in  substantial  agreement.  The  interrelationship  of  the  timescale  of  segmental  motion 
of  the  polymer  and  the  translational  diffusion  of  the  toluene  was  also  examined  and  it  was  found 
that  the  two  types  of  motion  seem  to  be  correlated  in  high  polymer  concentrated  solutions.  The 
toluene  reorientational  motion  was  found  to  be  much  faster  than  both  the  polymer  segmental 
motion  and  the  toluene  translational  diffusion  leading  to  the  conclusion  that  the  toluene 
reorientational  motion  is  uncoupled  from  these  two  motions. 

Keywords:  Segmental  motion.  Chain  Dynamics,  Polyisobutylene,  Toluene,  Hall  -  Helfand,  NMR, 
Free  Volume. 


Introduction 


Transport  through  polymers  is  a  subject  that  has  attracted  interest  during  the  past  three 
decades  because  of  the  applications  in  medicine,  biology  and  industry.  The  increased  number  of 
challenges  that  practitioners  are  faced  in  industry  and  applied  chemistry  has  Increased  the  need 
to  find  and  understand  more  about  the  mechanisms  that  take  place  in  a  polymer  when  a 
penetrant  molecule  diffuses  through  the  polymer  matrix. 

Numerous  diffusion  studies  have  been  reported  in  the  literature(1  '5)  and  it  has  been 
found  and  derTX)nstrated(®  by  many  researchers  that  free  volume  governs  the  transport  of  a 
penetrant  through  a  polymer  matrix.  Numerous  studies(10*13)  have  also  described  the  behavior 
of  a  polymer  In  solution  i.e.  when  the  latter  is  in  contact  with  a  penetrant.  A  variety  of  models(‘*0- 
''5)  that  describe  the  segmental  motion  of  a  polymer  in  the  bulk,  well  above  the  glass  transition 
temperature,  have  been  reported  in  the  literature.  In  general,  the  same  models  can  describe  the 
segmental  motion  of  the  polymer  in  solutioni"*®). 

In  order  to  obtain  a  better  understanding  about  the  mechanisms  of  transport  in  rubbery 
polymers  it  is  necessary  to  study  the  overall  translational  and  rotational  motions  of  the  penetrant 
as  well  as  the  local  segmental  motion  of  the  polymer.  To  accomplish  this  for  the  toluene-PIB 
system,  the  self  -  diffusion  coefficients  of  toluene  in  polyisobutylene  (PIB)  solutions  were 
determined  and  analyzed  as  a  function  of  temperature  and  concentration  using  the  two  most 
common  free  volume  theories,  namely  the  Fu]ita(®)  and  the  Vrentas  -  Duda(^'®)  theories.  In  the 
two  previous  papers(‘l7''>®)  in  this  series.  In  those  two  papers  the  free  volume  parameters  of  the 
toluene  -  PIB  solutions  were  determined  from  diffusion  and  compared  with  the  ones  obtained  from 
viscosity  data  of  the  same  solutions.  Also,  the  Fujita  theory  was  compared(^®)  with  the  Vrentas  - 
Duda  theory. 

In  this  present  work  the  local  segmental  motion  of  the  polymer  and  the  reorientational 
motion  of  the  penetrant  /  solvent  are  determined  using  Ti  NMR  measurements.  Previous 
studies(14,  19)  shown  that  segmental  motion  and  diffusion  of  a  penetrant  in  PIB  may  be 
correlated,  which  is  not  generally  the  case  for  other  polymers.  This  was  attributed  to  the  very 
efficient  molecular  packing  of  the  polymer(2®'22)  fhia  p,jgp|  density  of  PIB  relative  to  many  other 
polymers  implies  that  chain  reorientation  must  accompany  translational  motion  of  the  penetrant. 
The  objective  of  this  work  is  to  compare  the  time  scale  of  segmental  motion  of  the  polymer  with 
the  self  -  diffusion  coefficient  of  the  toluene  in  the  PIB  -  toluene  solutions  as  a  function  of 
concentration  to  see  if  they  are  indeed  correlated.  Segmental  motion  in  pure  PIB  has  been 


examined  in  detail  in  a  recent  report(^^)  and  that  information  will  serve  as  a  guide  in  the 
interpretation  of  the  concentration  dependent  results  to  be  presented  here. 

The  molecular  tumbling  of  liquid  toluene  has  been  studied(23-24)  earlier  via  spin  -  lattice 
relaxation  time  measurements  and  it  has  been  found  that  the  relaxation  time  of  the  ortho  is 
dominated  by  the  intramolecular  dipole  -  dipole  interaction  with  the  attached  proton.  Because  of 
the  r'®  dependence  of  the  dipole  -  dipole  interaction,  intermolecular  interactions  for  carbons 
directly  bound  to  a  hydrogen  atom(24)  ^an  be  neglected.  Other  mechanisms,  such  as  spin 
rotation  and  chemical  shift  anisotropy,  have  little  contribution  to  the  spin  -  lattice  relaxation  time.  It 
is  secondary  goal  to  examine  toluene  in  the  PIB  -  toluene  system  and  see  if  it  behaves  in  a  similar 
way  to  the  pure  liquid  state  and  try  to  obtain  additional  information  about  the  solution  free  volume. 

A  comparison  of  the  time  scale  of  reorientational  motion  of  the  solvent  with  the  segmental 
motion  of  the  polymer  is  required  to  determine  the  applicability  of  the  Kramers  relationship(25)  jn 
dilute  solution.  Zhu  et  al.(26)  have  shown  that  this  relation  fails  when  the  rotational  motion  of  the 
solvent  is  comparable  or  slower  than  the  time  scale  of  segmental  motion. 

Experimental 

The  polyisobutylene  sample  used  for  the  study  was  obtained  from  Cellomer  Associates, 
Inc.  of  Webster,  NY  in  the  form  of  thick  sheets  of  high  molecular  weight  (Cat.  #40E,  Lot  #02,  mol. 
wt.  1,000,000).  HPLC  Grade  toluene  was  obtained  from  FisherChemical  Fisher  Scientific.  NMR 
experiments  were  carried  out  on  a  Bruker  AC-200,  MSL-300  and  Varian  Unity  500  spectrometers, 
which  correspond  to  50.33  MHz,  75.5  MHz  and  125.7  MHz  ^®C  Larmor  frequencies  respectively. 
The  measurements  of  the  spin  -  lattice  relaxation  time  for  the  ortho  carbon  of  toluene  were 
performed  on  the  AC-200  and  the  MSL-300. 

A  simple  1 80®  -  x  -  90°  pulse  sequence  was  used  for  the  spin  -  lattice  relaxation  time 
measurements.  Two  to  three  hours  have  been  allowed  for  the  sample  to  reach  thermal 
equilibrium  before  each  measurement. 

R9?uI1s 


Spin  -lattice  relaxation  times  for  the  methylene  carbon  of  PIB  and  the  ortho  phenyl  carbon 
of  toluene  are  presented  in  Figures  1  -8  as  a  function  of  concentration  and  temperature. 


Interpretation 


The  equations  that  relate  the  spectral  density  function  and  the  spin  -  lattice  relaxation 
time  are  given  by: 

:i  =  W„+2W„  +  W,  (1, 

'  1 

where 

(2) 

,3, 

=  (4) 

I  'Urj 

The  methylene  C  -  H  distance  is  set  at  1.09  A. 

The  spectral  density  function  that  was  used  to  interpret  the  T-j  data,  is  the  Hall  -  Helfand 
function(15)  combined  with  restricted  anisotropic  rotational  diffusion(14.  27-30).  jhe  spectral 
density  function  is  shown  below: 

J(co)  =  {A  -I-  ^[(1  -  cos£f  +sm^£]  +  -^[(1  -  cos  2£f  +  sin^  2£]}}^^  (co)  + 

-L.— ,  l-cos(^  +  n;c)  2  sin(.^-n;t)  ,  sin(^  +  iMc),2,A.  ,  •. 

,  l-cos(2^-ni7r),2  .  ,sin(2^-n;r)  .  sm(2^-f  n;t),2,5L  .  . 

2^-nK  *~2Un,i  >  ^ 

where 


J“  (m)  =  2{[tS '  ( T^;  +  Tr')  -  0)2]’  +  [2x„>f  )->«  cos[i  tan-'  , 

2  XoXXoJ-hxP)-0)2 


J  ■  (CO)  =  2{[(To'  +  A„)(Toi  +  Ti ^  +  Xj-co2]2  +  [2(To;  +X„)co]2}-i/'» 

2(to/  +  ^n)0) 


®cos[— tan 


1.  -1 


(Xo^+^n)('Coi+'Cr'+^„)-0) 


] 


where 


-r-1  -  -r-1  J.  . 
^01  “  ^0  > 


(5) 


and 

A  =  4(3cos^A-1)2:  B  =  -(sin2  2A);  C  =  -(sin'^A)  (6) 

4  4  4 

The  angle  A  is  between  the  intemuclear  vector  and  the  axis  of  rotation,  the  latter  being  the  carbon 
-  carbon  backbone  bonds  for  this  case. 

According  to  this  model  two  correlation  times  describe  the  segmental  motion  of  the 
polymer.  One  of  them  is  xq  which  is  the  correlation  time  for  single  conformational  transitions  and 
the  other  is  xi  the  correlation  time  for  cooperative  or  correlated  transitions.  Each  of  these  two 
correlation  times  was  given  an  Arrhenius  temperature  dependence  with  activation  energy  EaO 
and  Eai  respectively.  Also  t  is  the  amplitude  over  which  restricted  rotational  diffusion  occurs  and 
Dir  is  the  rotational  diffusion  constant.  The  form  of  the  temperature  dependence  of  £  and  Djr  is 
shown  below: 


£  =  A  T®-®  deg 


(7) 


Djr  =  B  T  -  C  sec*"*  (8) 

When  the  Arrhenius  temperature  dependence  was  used,  seven  parameters  are  required, 
namely:  (Xoo)o.  (1:00)1.  EaO.  Eai.  A,  B.  C  and: 


TO  =  (Too)o  exp(EaO  /  R  T) 

(9) 

Ti  =  (Too)i  exp{Eai  /  R  T) 

(10) 

However,  it  is  not  necessary  to  float  all  seven  parameters  to  fit  the  experimental  data  as 
several  are  essentially  predetermined.  A  previous  study  of  PIB('*3)  in  solution  showed  that  the 
libration  does  not  change  from  the  bulk  state  to  solution.  Therefore  the  parameters  £  and  Dir  were 

kept  constant  and  equal  to  the  values  obtained  from  the  study(^4)  of  neat  PIB.  Also  during  the 
fitting  process  the  (xoo)o  and  {xoo)i,  did  not  seem  to  change  much  and  thus  they  were  kept 
constant  and  equal  to  the  values  obtained  from  the  study(14)  of  neat  PIB.  Therefore  only  two 
parameters  (EaO  and  Eai)  were  left  to  float  in  order  to  fit  the  spin  ■  lattice  relaxation  data  of  the 
methylene  carbon  of  PIB  in  toluene  -  PIB  solutions  as  a  function  of  temperature  and 
concentration.  The  results  are  shown  in  Table  1  while  the  rest  of  the  parameters  have  the  values: 
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(Xoo)o  =  2.1  lO'l^sec;  (Too)i  =  1.0  lO'^^sec  (11) 

^  =  5.15  tO-5  deg;  Djr  =  3.26  10^  T- 3.5  10^  sec''*  (12) 


The  curves  associated  with  the  fitting  are  compared  with  the  data  in  Figures  1-6. 


When  a  WLF  temperature  dependence  was  given  to  the  correlation  times,  there  are  six 
parameters  to  be  established:  (to)o.  (xi)0.  Djr.  c-i®  and  C2°  (ci®  and  C2°  are  the  WLF 

coefficients).  The  reference  temperature  was  chosen  to  be  To  =  298  K  and  accordingly; 


In 


JEo_ 

(■^o)o 


-2. 303c?  (T- To) 

C2  +  T-T0 


(13) 


T,  _-2.303c?(T-TJ 

(T,)o  C2°  +  T-T 


c®  -  ^0 

’  2.303f((l)J 

c0_n4>s) 


(15) 

(16) 


where  f(<j)s)  is  the  fractional  free  volume  of  the  solution  and  a  is  the  thermal  expansion  coefficient 
of  the  fractional  free  volume.  Bq  represents  the  minimum  void  size  necessary  for  a  conformational 
transition. 

Again  I  and  Djr  were  kept  constant  and  equal  to  the  values  obtained  for  the  neat  Pise'll) 
corresponding  to  the  WLF  temperature  dependence: 

i  =  5.4  T®-5  deg  and  Djr  =  2.1  10^  T  sec’"* 

Also  02°  was  set  equal  to  C2^  the  WLF  coefficient  obtained  from  the  studyC*^)  of  the  temperature 
dependence  of  the  self  -  diffusion  coefficient  of  toluene  in  PIB  -  toluene  solutions.  The  parameters 
that  were  floated  to  fit  the  data  are  (xo)o>  ("^1  )o  snd  c-|  The  values  obtained  are  shown  in  Table 
2.  Comparable  fittings  to  those  shown  in  Figures  1  -  6  were  obtained  from  the  analysis  of  the  Ti 
data  when  a  WLF  temperature  dependence  was  given  to  the  correlation  times. 

From  the  WLF  coefficients  ci  ^  shown  in  Table  2  and  equation  (15)  we  can  extract  values 
for  the  fractional  free  volume  of  the  solution.  The  results  are  shown  in  Table  3,  where  the 
fractional  free  volume  of  the  solution  (for  different  concentrations)  obtained  from  the  J-\  data  is 
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compared  with  the  corresponding  values  obtained  from  the  self  -  diffusion  coefficient  data 
analysis(17).  in  both  cases  the  reference  temperature  is  298  K.  In  Table  3,  Bd  represents  the 
minimum  void  size  necessary  for  a  molecule  to  undergo  a  displacement.  In  general  Bj  is  different 
from  Bo  because  they  are  associated  with  different  modes  of  motion:  Bj  is  associated  with 
translational  diffusion  while  Bq  is  associated  with  conformational  transitions. 

As  can  be  seen  from  Rgures  7  and  8  there  is  no  field  strength  dependence  in  the  T -j  data 
for  toluene.  Therefore,  the  extreme  narrowing  limit  equation  is  used  for  the  interpretation  of  the  Ti 
data  for  the  ortho  carbon  of  toluene  in  PIB  -  toluene  solution  and  is  shown  below : 

(17) 

h  *  CH 

where  S  =  1  /  2,  x  is  the  correlation  time  for  molecular  tumbling  and  rcH  =  1 .08  A  is  the  distance 
between  the  proton  and  the  carbon  on  the  toluene  phenyl  ring. 

Again,  at  first  an  Arrhenius  temperature  dependence  was  given  to  the  correlation  time: 

X  =  too  exp(E  /  R  T)  (18) 

where  E  is  the  activation  energy  for  rotation  and  Xoo  is  a  prefactor.  Two  parameters  were  floated  E 
and  Xoo.  During  the  fitting  procedure  it  was  found  that  Xoo  is  around  1.0x10'^^  sec  and  changes 
little  as  a  function  of  concentration.  Thus  the  value  of  Xoo  was  set  equal  to  1.0  10"^^  sec  and  E 
was  the  only  parameter  that  was  floated  as  concentration  changed.  The  parameters  obtained  are 
shown  in  Table  4.  The  fitting  curves  are  represented  by  the  solid  lines  in  Figures  7  and  8. 


Next  a  WLF  temperature  dependence  was  assigned  to  the  correlation  time  x: 


X  -2.3Q3c?(T-To) 
To  c°+T-To 


(19) 


where  xq  is  the  value  of  the  correlation  time  x  at  the  reference  temperature  To  which  again  was 
chosen  to  be  298  K.  ci  ®  and  C2®  are  the  WLF  coefficients  at  the  reference  temperature: 


C?  = 


B 


2.303f((t)J 


P  _  H^s) 


c;  = 


(20) 


a 


(21) 


where  f((|)s)  is  the  frectional  free  volume  of  the  solution  and  a  is  the  thermal  expansion  coefficient 
of  the  fractional  free  volume  of  the  solution.  B  represents  the  minimal  void  size  necessary  for 
rotation.  For  the  fitting  process  of  the  Ti  data  C2®  was  set  equal  to  C2*^  the  WLF  coefficient 
obtained  from  the  temperature  dependence  of  the  self  -  diffusion  of  toluene  in  PIB  -  toluene 
solutionsCi  7).  The  parameters  that  were  floated  in  order  to  fit  the  data  are  ci  0  and  tq.  the  values 
obtained  are  shown  in  Table  5.  The  fittings  obtained  are  comparable  to  those  shown  in  Rgures  7 
and  8  for  the  Arrhenius  temperature  dependence. 

Discussion 

In  this  work  three  different  kind  of  motions  were  considered:  the  self  -  diffusion  of  toluene 
in  high  molecular  weight  PIB,  the  molecular  rotation  of  toluene  and  the  segmental  motion  of  the 
polymer.  In  general  the  c-j®  WLF  coefficient  obtained  from  these  three  types  of  motion  should  be 
different,  because  of  the  different  hole  sizes  required  for  each  type  of  motion  to  occur,  i.e.,  B’s  (in 
general  ci^  =  B/2.303f((|»s)  and  C2®  =  f(<|)s)  /  a).  For  the  case  of  diffusion  we  have  Bd,  for  the 
toluene  rotation  B  and  for  the  segmental  motion  of  the  polymer  Bq. 

In  Table  3,  (f(<()s)  /  Bq)  obtained  from  segmental  motion  and  (f((j)s)  /  Bd)  obtained  from 
translational  diffusion  are  compared.  From  Table  3,  it  can  be  seen  that  the  two  parameters  are 
close  to  each  other,  considering  their  uncertainties:  of  the  order  of  ±  0.015  at  low  concentration, 
up  to  ±  0.05  at  high  toluene  concentrations.  The  agreement  is  very  good  noting  that  the  values 
are  extracted  from  two  completely  different  sets  of  data  that  describe  different  components  of  the 
solution.  The  toluene  self  -  diffusion  data  were  collected  by  observing  directly  the  toluene 
molecule  and  the  Ti  data  were  collected  by  observing  the  methylene  carbon  of  PIB  in  the  toluene 
-  PIB  solutions.  The  fractional  free  volume  of  the  solution  was  extracted  using  two  completely 
different  methods  that  observe  two  different  components  of  the  solution.  Therefore,  from  Table  3, 
since  the  fractional  free  volume  of  the  solution,  f(<j)s)  and  the  thermal  expansion  coefficient  of  the 

free  volume,  a,  are  unique  and  independent  of  the  type  of  motion  that  is  examined,  it  can  be 
concluded  that  Bd  =  Bq. 

This  conclusion  is  supported  by  a  previous  study(^^)  which  showed  that  the  segmental 
motion  of  PIB  is  correlated  with  the  translational  diffusion  of  CO2.  The  fact  that  Bd  «  Bq  implies 
that  the  same  relationship  holds  for  the  toluene  -  PIB  system,  when  a  toluene  molecule  is 
surrounded  by  polymer  segments.  Also  the  value  of  obtained  from  the  Vrentas  -  Duda 
theoryO®)  is  0.892  which  corresponds  to  a  polymer  jumping  unit  with  molecular  weight  94.3.  This 
is  very  close  to  the  molecular  weight  of  toluene  which  is  92.  The  fact  that  the  two  nearly  coincide 
with  each  other  means  that  the  two  jumping  units  will  need  the  same  minimal  free  volume  for 


translation  (for  toluene)  and  a  conformational  transition  (for  the  polymer).  This  also  indicates  that 
Bd  should  equal  Bq- 


Another  way  to  check  to  see  if  the  two  motions  are  correlated,  is  to  directly  compare  the 
correlation  time  for  translational  diffusion  XQ  and  the  correlation  time  for  segmental  motion 
(correlation  time  for  single  conformational  transition).  The  first  can  be  obtained  from(2^): 


(22) 


where  D  is  the  self  -  diffusion  coefficient  of  toluene,  b  is  the  distance  of  closest  approach  which  is 
between  2.2  A  -  3.3  A.  Using  2.7  A  as  the  distance  of  closest  approach  and  the  self  -  diffusion  of 
toluene  in  PIB  -  toluene  solution,(^^)  the  correlation  time  for  translational  diffusion  is  calculated 
and  shown  in  Table  6.  The  correlation  time  for  a  conformational  transition,  xq,  is  calculated  from 
equations  (9),  (1 1)  and  Table  1.  xq  is  also  shown  in  Table  6.  Looking  at  Table  6  it  can  be  seen 
that  the  two  correlation  times  are  close  to  each  other  at  0.045  toluene  volume  fraction  which 
means  at  this  low  penetrant  concentration  the  two  motions  (segmental  and  translational  diffusion) 
are  correlated.  As  we  go  to  higher  concentrations  the  difference  between  the  xq  correlation  time 
and  the  xq  becomes  larger  up  to  two  orders  of  magnitude  at  0.712  toluene  volume  fraction  (see 
Table  7).  This  might  be  expected  since  as  we  go  towards  higher  concentrations  a  toluene 
molecule  has  fewer  polymer  molecules  in  its  neighbourhood  and  it  is  not  necessary  for  a  polymer 
segment  to  move  in  order  a  toluene  molecule  to  translate. 

Thus  we  see  that  the  two  processes  (diffusion  and  conformational  transition)  require 
similar  minimum  free  volumes  in  order  to  occur.  The  diffusion  of  small  molecules  in  PIB  is 
affected  by  the  free  volume  which  in  turn  influences  the  molecular  mobility  of  the  polymer.  Also 
the  time  scale  of  the  two  motions  are  correlated  at  low  concentrations  of  penetrant. 

The  activation  energies  EaO  and  Eai  as  a  function  of  the  volume  fraction  of  toluene  in 
PIB  are  shown  in  Figure  9.  The  WLF  coefficient  c-]®  obtained  from  Ti  data  analysis  and  ci^ 
obtained  from  the  self  -  diffusion  data  analysis  are  shown  in  Rgure  10.  The  trends  here  are  as 
expected  and  the  agreement  between  the  two  WLF  fittings  is  encouraging. 


The  spin  -  lattice  relaxation  data  on  the  ortho  carbon  of  toluene  do  not  provide  as  much 
information  about  the  motion  of  toluene  as  the  corresponding  data  on  PIB.  However,  at  all 
concentrations  the  correlation  time  for  toluene  reorientation  is  two  or  more  orders  of  magnitude 
faster  than  segmental  motion  and  translational  diffusion.  Thus  the  reorientational  motion  of  the 
toluene  is  uncoupled  from  the  other  two  motions.  The  WLF  analysis  of  the  temperature 


dependence  also  leads  to  different  values  of  ci  0  and  thus  B  further  supporting  the  separation  of 
these  motions.  A  comparison  of  f((|)s)  /  with  f((j)s)  /  B  is  shown  in  Table  8.  Toluene  reorientation 
would  appear  to  be  a  much  more  rapid  motion  involving  a  smaller  hole  size.  The  separation  of  the 
time  scale  of  solvent  reorientation  from  segmental  motion  at  low  PIB  concentration  also  points  to 
a  likely  Kramers  dependence  of  segmental  motion  on  solvent  viscosity  in  dilute  solutions  though 
only  a  single  solvent  is  considered  in  this  study. 
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Table  1 


Activation  energies  for  conformational  transitions  as  a  function  of  concentration  of  toluene  in  PIB. 


Volume  Fraction  of  toluene  in 

PIB. 

EaO  (kJ  /  mole) 

Eai  (kJ  /  mole) 

0.045 

34.0 

34.5 

0.104 

33.4 

33.7 

0.205 

32.7 

32.0 

0.327 

31.9 

30.0 

0.409 

31.7 

28.4 

0.712 

30.4 

25.7 

Table  2 

Fitting  parameters  of  the  spin  -  lattice  relaxation  data  when  a  WLF  temperature  dependence 

given  to  the  correlation  times. 

Volume  Fraction 

of  toluene  in  PIB 

(to)o  (sec) 

(Tl)o(sec) 

c-|9 

C2°  (K) 

0.045 

3.00  10-8 

1.2010-8 

4.30 

155.06 

0.104 

1.80  10-8 

8.4010-9 

3.80 

153.50 

0.205 

8.85  10-9 

4.25  10-9 

3.40 

175.17 

0.327 

6.23  10-9 

1.8810-9 

2.75 

183.65 

0.409 

3.40  10-9 

9.7010-‘'0 

2.65 

197.92 

0.712 

3.00  10-9 

4.00  10*'' 9 

1.85 

245.45 

13 


Table  3 

Fractional  free  volume  comparison  of  the  PIB  -  toluene  solutions  obtained  from  the  diffusion  data 
analysis(1 7)  and  from  the  T i  data  presented  in  this  work. 


Volume  fraction  of  toluene  in 

f(<|)s)  /  Bd 

1(<l>s) !  Bq 

P\B^s 

From  Diffusion  data 

FromT-i  data 

0.045 

0.090 

0.101 

0.104 

0.105 

0.114 

0.205 

0.127 

0.127 

0.327 

0.155 

0.158 

0.409 

0.190 

0.164 

0.712 

0.270 

0.235 

Table  4 

Fitting  parameters  of  the  ortho  carbon  Ti  data  when  an  Arrhenius  temperature  dependence 

Volume  fraction  of  toluene  in 

given  to  the  correlation  time. 

Too  (sec) 

E(kJ/mole) 

PIB  <j)s 

0.045 

1.010-''4 

22.73 

0.104 

1.010-14 

21.84 

0.205 

1.010-14 

20.76 

0.327 

1.010-14 

18.50 

0.409 

1.010*14 

16.77 

0.712 

1.010-14 

14.00 
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Table  5 

Fitting  parameters  of  the  ortho  carbon  T i  data  when  a  WLF  temperature  dependence  was  given 

to  the  correlation  time. 


Volume  fraction  of 
toluene  in  PIB  <|)s 

TO  (sec)  at  To  298K 

ci0at298K 

02°  (K)  at  298 

0.045 

9.72  10-''1 

2.45 

155.06 

0.104 

6.8010-'' I 

2.35 

153.50 

0.205 

^loio-n 

2.25 

175.17 

0.327 

1.6710-''1 

2.10 

183.65 

0.409 

8.40  10-''2 

2.00 

197.92 

0.712 

2.90  10- ■'2 

1.81 

245.45 
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Table  6 

Comparison  of  the  correlation  times  for  diffusion,  xd.  and  for  single  conformational  transitions  tq 

for  0.045  toluene  volume  fraction. 


•  (cm2  /  sec)  Toluene 
Volume  Fraction  = 

0.045 

T(K) 

to  (sec) 

XD  (sec) 

2.5510-8 

298 

1.91510-8 

2.859  10-8 

6.6610-8 

310 

1.12610-8 

1.09510-8 

1.2610-7 

320 

7.453  10-9 

5.78610-9 

2.0410-7 

329 

5.254  10-9 

3.57410-9 

3.11  10-7 

339 

3.642  10-9 

2.344  10-9 

4.5510-7 

348 

2.666  10-9 

1.60210-9 

6.31  10-7 

357 

1.98210-9 

1.155  10-9 

8.1410-7 

367 

1.451  10-9 

8.956  lO-"*  8 

Table  7 


16  ■ 


Comparison  ol  the  correlation  times  for  diffusion,  tQi  and  for  single  conformational  transitions  tq 

for  0.712  toluene  volume  fraction. 


D  (cm2  /  sec)  Toluene 
Volume  Fraction  = 

0.712 

T(K) 

to  (sec) 

tD  (sec) 

3.1  10-6 

249 

5.008  10-6 

2.352  10-16 

4.0210-6 

259 

2.841  10-6 

1.813  10-16 

5.0910-6 

269 

1.681  10-8 

1.432  10-16 

5.79  10-6 

280 

9.854  10-9 

1.259  10*16 

8.0010-6 

298 

4.478  10-9 

9.11310-11 

9.9010-6 

310 

2.785  10-9 

7.364  10-11 

1.19  10-6 

320 

1.92610-9 

6.339  10*11 

1.41  10-5 

329 

1.409  1  0-9 

5.17010-11 

1.4610-6 

339 

1.01510-9 

4.993  10*11 

1.81  10-6 

348 

7.681  10-16 

4.028  10*11 

1.9410-6 

357 

5.893  10-16 

3.758  10*11 
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Table  8 


Fractional  free  volume  comparison  of  the  PIB  -  toluene  solutions  obtained  from  the  diffusion  data 
analysis('*  7)  and  from  the  T i  data  of  the  ortho  carbon  - 1 3  of  toluene  presented  in  this  work. 


ction  of  toluene  in 

PIB  (t)s 

f((j>s)  /  Bd 

From  Diffusion  data 

f(<j)s)/B 

FromTi  data 

0.045 

0.090 

0.177 

0.104 

0.105 

0.185 

0.205 

0.127 

0.193 

0.327 

0.155 

0.207 

0.409 

0.190 

0.217 

0.712 

0.270 

0.240 

Houre  Captions 

Figure  1:  T-)  of  the  methylene  carbon  of  PIB  in  toluene  -  PIB  solutions  as  a  function  of 
temperature  and  Larmor  frequency.  The  volume  fraction  of  toluene  in  PIB  is  0.045. 
An  Arrhenius  temperature  dependence  was  given  to  the  correlation  times.  The  solid 
line  corresponds  to  the  fit  for  the  Hall  -  Helfand  function  combined  with  anisotropic 
restricted  rotational  diffusion. 

Figure  2:  Ti  of  the  methylene  carbon  of  PIB  in  toluene  -  PIB  solutions  as  a  function  of 
temperature  and  Larmor  frequency.  The  volume  fraction  of  toluene  in  PIB  is  0.104. 
An  Arrhenius  temperature  dependence  was  given  to  the  correlation  times.  The  solid 
line  corresponds  to  the  fit  for  the  Hall  -  Helfand  function  combined  with  anisotropic 
restricted  rotational  diffusion. 

Figure  3;  T-|  of  the  methylene  carbon  of  PIB  in  toluene  -  PIB  solutions  as  a  function  of 
temperature  and  Larmor  frequency.  The  volume  fraction  of  toluene  in  PIB  is  0.205. 


An  Arrhenius  temperature  dependence  was  given  to  the  correlation  times.  The  solid 
line  corresponds  to  the  fit  for  the  Hall  -  Helfand  function  combined  with  anisotropic 
restricted  rotational  diffusion. 

Figure  4:  Ti  of  the  methylene  carbon  of  PIB  in  toluene  -  PIB  solutions  as  a  function  of 
temperature  and  Larmor  frequency.  The  volume  fraction  of  toluene  in  PIB  Is  0.327. 
An  Arrhenius  temperature  dependence  was  given  to  the  correlation  times.  The  solid 
line  corresponds  to  the  fit  for  the  Hall  -  Helfand  function  combined  with  anisotropic 
restricted  rotational  diffusion. 

Figure  5:  Ti  of  the  methylene  carbon  of  PIB  in  toluene  -  PIB  solutions  as  a  function  of 
temperature  and  Larmor  frequency.  The  volume  fraction  of  toluene  in  PIB  Is  0.409. 
An  Arrhenius  temperature  dependence  was  given  to  the  correlation  times.  The  solid 
line  corresponds  to  the  fit  for  the  Hall  -  Helfand  function  combined  with  anisotropic 
restricted  rotational  diffusion. 

Figure  6:  Ti  of  the  methylene  carbon  of  PIB  in  toluene  -  PIB  solutions  as  a  function  of 
temperature  and  Larmor  frequency.  The  volume  fraction  of  toluene  in  PIB  is  0.712. 
An  Arrhenius  temperature  dependence  was  given  to  the  correlation  times.  The  solid 
line  corresponds  to  the  fit  for  the  Hall  -  Helfand  function  combined  with  anisotropic 
restricted  rotational  diffusion. 

Figure  7:  Ti  of  the  ortho  Qf  toluene  in  PIB  -  toluene  solutions  as  a  function  of  temperature 
concentration  and  Larmor  frequency.  The  numbers  shown  next  to  the  curves 
represent  the  volume  fraction  of  toluene  in  the  solution.  An  Arrhenius  temperature 
dependence  was  given  to  the  correlation  time  x.  The  solid  lines  correspond  to  the  Ti 
values  calculated  from  equation  1 7. 

Figure  8:  T-|  of  the  ortho  of  toluene  in  PIB  -  toluene  solutions  as  a  function  of  temperature 
concentration  and  Larmor  frequency.  The  numbers  shown  next  to  the  curves 
represent  the  volume  fraction  of  toluene  in  the  solution.  An  Arrhenius  temperature 
dependence  was  given  to  the  correlation  time  x.  The  solid  lines  correspond  to  the  Ti 
values  calculated  from  equation  1 7. 

Figure  9;  Activation  energies  as  a  function  of  concentration  obtained  from  the  Arrhenius 
temperature  dependence  of  the  correlation  times. 


Rgure  10:  WLF  coefficients  as  a  function  of  concentration. 


A  NMR  Study  of  Penetrant  Diffusion  and  Polymer  Segmental 
Dynamics  in  Rubber. 

A.  Bandis,  P.T.  Inglefield*.  A.A.  Jones  and  W.-Y.  Wen.  Chemistry 
Department,  Clark  University,  Worcester,  MA,  01610. 

Introduction 

.  Diffusion  thcorics<i-3)  utilizing  free  volume  concepts  are  known  to 
accurately  model  transport  in  both  rubberyC^)  and  glassy  polymers(5).  A 
number  of  studies  have  been  made  concerning  diffusion  in  polyisobutylcne 
(PIB)(6*9),  PIB  is  a  special  polymer  because  of  its  low  permeability  to  gases 
and  liquids.  The  low  permeability  may  come  from  the  fact  that  PIB  has  a 
low  specific  volumc(*®*i2).  Previous  studies  have  shown  that  segmental 
motion  and  diffusion  in  PIB  are  conelaicdtis-M)  while  this  is  not  generally 
the  case  for  other  polymers.  In  the  present  work  the  pulsed  field 
gradient(PFG)  NMR  technique  was  used  to  directly  measure  the  self¬ 
diffusion  coefficient  of  toluene  in  PIB  as  a  function  of  concentration  and 
temperature.  The  spin-lattice  relaxation  time(Ti)  of  the  methylene 
carbon  of  PIB  was  also  determined  to  simultaneously  monitor  the  segmental 
motion  of  the  polymer  in  PIB/Toluene  solutions. 

Experimental 

The  NMR  measurements  were  made  on  a  Brukcr  MSL300 
spectrometer.  In  the  PEG  technique  an  applied  magnetic  field  gradient  was 
varied  from  14  gauss/cm  up  to  140  gauss/cm.  The  Ti*s  of  the  methylene 
carbon  of  PIB  were  determined  from  about  5wt  %  up  to  70wt  %  of  toluene 
in  PIB. 

Results  andinterpretarion 

In  figure  1  the  self  diffusion  coefficient  of  toluene  in  PIB  is  shown 
as  a  function  of  temperature  and  concentration.  In  figure  2  the  13C  Tfs  of 
the  backbone  methylene  of  PIB  are  shown  as  a  function  of  temperature  and 
Larmor  frequency  for  70wt  %  toluene  in  PIB. 

The  temperature  dependence  of  the  self-diffusion  coefficients  of 
toluene  in  PIB  were  fitted  using  the  familiar  WLF  equationtt^): 

DT„  ^(T-T.) 

where  is  the  volume  fraction  of  toluene  in  PIB,  f(<t»s)  is  the  fractional  free 
volume  of  the  solution  at  the  reference  temperature  Tq,  a((t)*)  is  the 


T(K) 

Fig.  1.  Temperature  Dependence  of  the  toluene  self  •  diffusion  coefficient  in 
PIB.  The  solid  lines  represent  fittings  from  the  WLF  equation  presented  in 
the  text 


Fig.  2.  Spin  -  Lattice  relaxation  times  of  the  methylene  carbon  of  PIB  as  a 
function  of  temperature  and  Larmor  frequency.  The  concentration  of  toluene 
is  about  70%  wt  The  solid  lines  represent  fitting  from  the  Hall-Helfand 
model  combined  with  a  restricted  anisotropic  rotntional  diffusion.  An 
Arriienius  temperature  dependence  was  given  to  the  correlation  limes. 


The  data  were  also  fit  as  a  function  of  concentration  at  a  constant 
temperature  using  the  Fujita  theory(0: 

n  Bff  fU 

■  D_  B/f.-f,)i|>.  Bj  B, 

"d  l^d  "d  "d 

where  the  parameters  are  the  same  as  above  except  that  fp  is  the  average 
fractional  free  volume  of  pure  polymer,  f,  is  the  fractional  free  volume  of 


fractional  free  volume  expansion  factor  and  Ba  is  a  constant  which 
corresponds  to  the  minimum  hole  required  for  a  given  molecule  to  diffuse. 
The  parameters  f((l)s)/Bd  and  obtained  from  data  fitting  are: 


<t>s 

f(«5)/Ba 

a(4>s)/Bd  (deg*i) 

CjD 

0 

0.0818 

5.34  10*4 

5.308 

153.18 

0.045 

0.0904 

5.83  10-4 

4.803 

155.06 

0.104 

0.105 

6.84  10*4 

4.135 

153.51 

0.205 

0.127 

7.25  10-4 

3.419 

175.17 

0.327 

0.155 

8.44  10-4 

2.801 

183.65 

0.409 

0.190 

9.60  10-4 

2.285 

197.92 

0.712 

0.270 

1.10  10-3 

1.608 

245.45 

1.000 

0.396 

1.70  10-3 

1.097 

232.94 

toluene  in  the  solution  at  the  limit  of  zero  toluene  concentration(i6)  and  Dq  is 
the  diffusion  coefficient  of  toluene  at  the  limit  of  zero  pcnetrant(toIuene) 
concentration.  The  parameters  varied  are  Do,  (fp/Ba)  and  (fs/Ba).  The  values 
for  the  parameters  are: 


T(K) 

(f./Ba) 

(fp/Bd) 

Do  (cm2 /sec) 

298 

0.273 

0.08735 

9.08  10-9 

310 

0.300 

0.0978 

2.55  10-8 

320 

0.333 

0.1070 

5.15  10-8 

329 

0.364 

0.1163 

8.70  10-8 

339 

0.381 

0.1250 

1.55  10-7 

349 

0.398 

0.1320 

2.30  10-7 

358 

0.423 

0.1410 

3.40  10-7 

368 

0.449 

0.1445 

4.10  10-7 
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The  data  were  also  fined  using  the  Vrentas  -  Duda  theoryP): 


D,  =Doexp(^)exp[ 


_ -(CD.V‘  +CDi4Vi) 

o),(^)(K„  -T„  +  T)  +  to,(^^XK^  -T.,  +  T) 
Y  Y 


Where  Di  is  the  self  diffusion  coefficient  of  the  penetrant  and  the  other 
parameters  are  defined  elsewherefiT).  The  values  of  the  parameters  are 
given  in  reference  17  except  the  values  of  Ki2/y  and  K22  -  Tg2  that  were 
obtained  from  an  NMR  study  of  neat  PIBd^l;  their  values  are  4.351  10-4 
cirf/gr  K  and  -97.6  K  respectively.  The  only  parameter  that  was  floated  is  ^ 


which  is  the  ratio  of  molar  volume  of  a  solvent  jumping  unit  over  the  molar 
volume  of  a  polymer  jumping  unit: 

V?(0)^  MX 
^  M,^v; 

where  ViO(O)  is  the  solvent  molar  volume,  V2j  is  the  jumping  unit  volume, 
Ml  is  the  molecular  weight  of  the  solvent  and  M2j  is  the  molecular  weight 
of  the  polymer  jumping  unit  and  Vi*,  V2*  are  the  specific  critical  hole  finec 


volumes  of  pure  toluene  and  PIB  respectively.  The  value  of  ^  is  about 
0.892  ±  0.08  which  corresponds  to  a  polymer  jumping  unit  of  molecular 
weight  94.3.  A  fitting  of  the  diffusion  coefficient  as  a  function  of 
concentration  at  two  different  temperatures  for  the  two  different  theories  is 


shown  in  figure  3. 

The  Ti  data  of  the  methylene  carbon  of  PIB  in  the  toluene  solutions 
were  fitted  using  the  Hall  -  Helfand  model  combined  with  librationi^).  The 
data  were  fitted  using  both  a  WLF  approach  with  values  obtained  from  the 
diffusion  data  and  an  Arrhenius  temperature  dependence  for  the  correlation 
times. 

One  of  the  fittings  is  shown  in  figure  2.  The  parameters  obtained  for 
the  Arrhenius  temperature  dependence  show  that  the  activation  energies  for 
the  Hall-Helfand  processes  arc  close  and  decrease  slightly  with  increasing 
toluene  content  with  values  from  34  to  26KJ/mole.  Parameters  that  were 
constant  at  all  concentrations  arc: 


to  -  (too)o  exp(Eao/RT),  ti  =  (t^)i  exp(Eai/RT).  t  =  A  TO.5,  Di,  =  B  T  -  C 
(t^)o  =  2.1  10-^^  sec,  (t^)i  =  1.0  10-14  sec,  t  =  5.15  TO-5  deg.  Dir  =  3.26 


10?  T  -  3.5  109  sec-i.  Where  to  is  the  correlation  time  for  single 


conformational  transitions,  1 1  the  correlation  time  for  cooperative 

transitions,  EaO,  Eal.  are  the  activation  energies  for  to  and  ti  respectively,  C 
the  angular  amplitude  over  which  rotational  diffusion  occurs  and  Djr  the 
rotational  diffusion  coefficient. 


Fig.3  Fitting  of  the  diffusion  coefficient  of  toluene  in  PIB  as  a  function  of 
concentration  at  two  temperatures  using  Fujita  and  Vrentas  -  Duda  theory. 


Discussion 

It  is  clear  from  the  above  analysis  that  the  self-diffusion  of  toluene 
in  PIB  is  well  described  by  free  volume  theories.In  this  work  two  types  of 
motion  are  investigated:  self-diffusion  of  toluene  in  high  molecular  weight 
PIB  and  segmental  motion  of  the  polymer.  In  general  the  ci^WLF 
coefficient  obtained  from  these  two  types  of  motion  should  be  different, 
because  of  the  different  hole  sizes,  ie,  B’s  (in  general  Cj®  =  B/2.303f  and 

C2O  =  f/a).  For  the  case  of  diffusion  we  would  have  Ba,  while  for 
segmental  motion.  B.  It  seems  from  this  study  that  Ba  is  equal  to  B  since 
the  diffusion  and  the  Ti  data  can  be  fitted  using  the  same  WLF  coefficients. 

since  f  (fractional  free  volume  of  the  solution)  and  a  (thermal  expansion 
coefficient  of  the  fractional  free  volume)  arc  unique  and  independent  of  the 
type  of  motion  that  is  examined.  This  is  supported  by  a  previous  study(i3) 
that  showed  that  the  segmental  motion  of  PIB  is  correlated  with  the 
translational  diffusion  of  CO2.  The  fact  that  Bd  =  B  shows  that  the  same 

thing  happens  for  the  toluenc/PIB  system.  Also  the  value  of  obtained 
from  the  Vrentas  -  Duda  theory  is  0.892  which  corresponds  to  a  polymer 
jumping  unit  with  molecular  weight  94.3.  This  is  very  close  to  the 
molecular  weight  of  toluene,  which  is  92.  The  fact  that  the  two  are  very 
close  to  each  other  means  that  the  two  jumping  units  will  need  the  same 
minimal  free  volume  for  translation  (for  toluene)  and  a  conformational 
transition  (for  the  polymer).  This  indicates  that  Bd  =  B.  It  seems  that  the 
diffusion  of  small  foreign  molecules  in  PIB  is  affected  by  the  free  volume 
which  in  turn  influences  the  molecular  mobility  of  the  polymer. 
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Abstract 

:czzr.‘~ 

,  '  * concentration  dependence  nf  tho 

bath  the  Puika  and  the  VentrasJOuda  free  volume  thiries  a^m  T 
interpreted  with  the  WLF  equation  THp  w  .  *®'"P®fature  dependence  was 

temperature  dependencies  of  the  toluene  lerdtsT""'' 

free  volume  parameters  are  in  a  good  agreement  resulting 

viscosity  data  on  the  same  system  Comoa  •  ^  extracted  from  the  analysis  of 

theories  and  the  Fuifta  free  volume  ParameterruXL^l'^^^^  ‘  V 

volume  parameters.  xtracted  from  the  Vrentas  -  Duda  free 

polyisobutylene  (PieririeTrf^^  determined  for  the  methylene  carbon  of 

solutions.  The  Hall  -  Helfand  correct  f  ’  P^'V-butylene 

diffusion  was  used  to  treat  the  data  oTth  "  combined  with  restricted  anisotropic  rotational 

polymer  segmenfal  melon  aIoL  '  ' 

-lnm.io„e.olue„e:Llrs:~i“ 

Studied.  Both  models  described  satisfart  ’i  ♦»,  '■®'"®  "^’’^wing  limit  for  the  solutions 

'a-«..,on..oc.Promr.o:r:re:~ 

—htr:  ~  -  ~-:rc.:'^^ 

-eev.raeedfmm.her,. .for, heme, 


““z,rrj:;:rr=— — 

TKrs.irMmeCwes  of  moto  seem  .Ob.  correlated  in  high  po^mer<»ncemra., on 

lworr,oiuene  reoden.e.^r«i  mCion  was  found  fo  be  much  faster  dtan  both  the  ^mer 
segmental  motion  and  the  toluene  banslational  d«fusion  leading  to  the  conclusion  W  the  toluene 
reorienlational  motion  is  uncoupled  from  Urese  two  moMns. 

rcnuwnrds:  Polyisobutylene,  Toluene,  Pulse  Field  Gradient  Nuclear  Magnetic  Resonance, 
n-rffusion.  Free  Volume  Theoiy.  Segmental  Motion. 


Introduction 

The  dmusion  of  small  molecules  into  polymers  has  been  the  subiec.  of  many 
,...bons,becauseo„^^ 

segmental  motlonP)  in  the  pure  rubber  by  NMR. 

There  are  a  variety  of  techniques  for  measuring  the  self  -  diffusion  coefficients  of  a 
I  „  matrix  These  incfude  sorption,  permeation,  light  scattering,  radioactive 
penetrant  in  a  polymer  matrix.  Tn  /mmp\  tkp  Pulsed  Field 

rr-r:==ri": 

precision  of  the  measurements  and  the  ability  to  be  applied  to  both  liqukf  and  solid  ^ 

Ly  limitation  of  the  technique  is  that  the  lowest  diffusion  coeHicien.  lha.  can  be 
determined  is  around  tO'P  cm2  /  sec  to  1 0'l  °  cm2  /  sec. 
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In  addition  to  the  advantages  of  the  Pulsed  Gradient  technique  Hsell,  with  NMR  the 
reonentatK.nal  dynamics  of  the  both  the  polymeric  and  penetrant  components  of  the  solution  can 
be  measured  rndependently  on  the  same  sample.  Thus  both  translational  and  rotattonal  motlo;,  of 
the  penetrant,  and  segmental  motion  of  the  polymer  can  be  determined  in  the  same  samples 
With  modern  NMR  instrumentation,  these  measurements  can  be  made  over  a  wide  range  cl 
temperatures  and  concenlrations  which  lead  to  a  fairly  complete  description  of  local  dynamics  It 
la  the  goal  of  this  work  to  present  the  results  of  such  a  complete  set  of  NMR  measurements  on  a 
single  system  beginning  with  translational  diffusion  of  the  penetrant. 

Diffusion  in  polyisobulylene  (PIB),  has  been  studied  before  using  different  techniquesP  ■ 

.  In  references  (3)  and  (4)  the  study  of  the  diffusion  coefficienls  of  six  hydrocarbons  (propane,  n 
-  butane,  Isobulane,  n  -  pentane,  isopentane  and  neopentane)  in  PIB  was  reported  using  the 
sorption  method.  The  sc#  -  diffusion  of  benzene  in  PIB  has  been  studied  in  reference  (5)  using 
the  pulsed  field  gradient  technique  over  a  limited  range  of  temperatures  and  concentrattons.  In 
reference  (6)  the  diffusion  of  radioactively  lagged  cetane  In  PIB  has  been  studied  and  finally  in 
reference  (7)  the  diffusion  of  toluene  In  butyl  rubber  has  been  studied  using  the  sorption  method. 

In  all  cases  the  data  were  analyzed  using  the  Fullta  tree  volume  theory.  Viscosity  data  on 
poviscbutylene  solutions  are  also  availableP ' '  D.  In  references  (8  -  to)  the  intrinsic  viscosities 
of  PIB  solutions  in  various  solvents  have  been  studied  and  their  relations  to  polymer  chain 

structure  and  to  thermodynamic  parameters  governing  the  interaction  between  the  polymer  and 
the  solvent  have  been  made. 


Free  volume  theories  have  been  used  widely  for  the  interpretation  of  the  concentration 
dependence  of  the  self  -  diffusion  of  many  small  molecules  into  polymer  matrlces('2  - 15)  one 
frequently  used  theory,  developed  almost  thirty  years  ago(12),  is  the  Fujita  theory  and  a  more 
recent  free  volume  theory  is  that  of  Vrentas  and  Dudal13  •  IS),  phese  two  free  volume  theories 
differ  in  several  respects  and  it  will  be  of  value  to  compare  them  In  the  context  of  a  rather 
complete  data  set,  not  only  including  information  of  translational  motion  of  the  penetrant  but  also 
reonentational  motion  of  both  components:  polymer  and  penetrant.  The  molecular  level 
information  available  from  NMR  will  provide  insight  into  some  of  the  concepts  introduced  into 

these  theories  which  are  at  least  at  first  glance  based  on  very  local  descriptions  of  polymer  and 
solvent  dynamics. 


The  Fujita  theoiy  begins  with  the  empirical  result  of  Doolittle(16  - 17)  showed  that 
the  viscosity  ti  of  a  liquid  of  low  molecular  weight  is  given  by; 


lmi  =  lnA  +  B/f 


0) 


where  A  and  B  are  empirical  constants  for  a  given  liquid  and  f  is  the  fractional  free  volume, 
defined  to  be  equal  to  the  ratio  of  the  volume  of  the  space  not  occupied  by  the  constituent 
molecules,  vf,  over  the  total  volume  of  the  solution,  v: 


(2) 

V 

V  =  VO  +  Vf 

(3) 

where  vq  is  the  volume  occupied  by  the  molecules  in  accordance  with  Van  der  Waals  radii  and 
vibrational  motions.  Equation  (1)  measures  the  volume  of  a  molecule  relative  to  the  volume  of  a 
void  next  to  it.  When  the  volume  of  the  void  Is  larger  than  the  volume  of  the  molecule  the  viscosity 
is  low.  When  the  volume  of  the  void  is  smaller  than  the  volume  of  the  molecule  the  viscosity  is 
high.  Another  assumption  that  is  made  is  that  the  Doolittle  relation  applies  to  each  jumping  unit  in 
polymer  -  solvent  systems('*®).  Jumping  units  are  imagined  to  be  of  the  size  of  small  penetrant 
molecules  or  repeat  units  in  the  case  of  the  polymer.  Also  the  theory  relates  the  mobility  of  the 
polymer  and  the  solvent  jumping  units  to  the  free  volume  per  unit  volume  of  solution. 

The  Vrentas  -  Duda  theory  is  based  on  the  Cohen  -  Turnbull  resulK"*^)  showing  that  the 
mobility  of  a  molecule  in  a  pure  liquid  is  related  to  the  average  free  volume,  vf,  per  molecule 
through  the  following  equation: 

In  m  =  E  -  y  V*  /  Vf  (4)  where  y 

Is  a  numerical  factor  between  0.5  and  1,  also  known  as  overlap  factor,  that  is  introduced  because 
a  given  free  volume  may  be  available  to  more  than  one  molecule  and  v*  is  the  minimum  volume 
of  a  hole  that  a  molecule  must  find  in  its  vicinity  in  order  for  it  to  move.  According  to  Cohen  and 
Turnbull  molecules  are  resident  in  cages  bound  by  their  neighbors  for  the  majority  of  the  time. 
When  a  density  fluctuation  occurs  and  a  hole  of  enough  size  opens  within  the  molecular  cage,  the 
molecule  jumps  into  the  hole.  Diffusional  motion  is  considered  to  be  translation  of  a  molecule 
across  the  space  within  its  cage  and  bulk  diffusion  then  occurs  as  a  result  of  redistribution  of  the 
free  volume  within  the  liquid.  The  above  equation  was  derived  under  the  assumption  that 
redistribution  of  free  volume  (holes)  at  constant  volume  requires  no  change  in  energy.  Equation 
(4)  demonstrates  that  when  the  minimum  volume  of  a  hole  that  a  molecule  must  find  for  a  jump  to 
occur  is  larger  than  the  average  free  volume  per  molecule  then  the  mobility  of  the  penetrant  is 
small.  Conversely,  when  v*  is  smaller  than  vf  then  the  molecule  can  jump  to  a  neighbouring 
position  and  the  mobility  is  high. 

Vrentas  -  Duda(^3)  formulated  their  theory  based  on  the  following  assumptions  as  stated; 


(I)  Vf'  IS  the  free  volume  equal  to  the  space  not  occupied  by  the  actual  molecules. 

(II)  The  free  volume  available  to  each  jumping  unit:  solvent  molecule  or  polymer  segment,  is  equal 
to  Vf'  /  N  where  N  is  the  total  number  of  solvent  molecules  plus  polymer  segments  and  Vf'  is 
the  total  free  volume  of  the  solution. 

(III)  The  critical  hole  volume,  v*,  of  a  jumping  unit  is  equal  to  its  intrinsic  volume  at  zero  Kelvin. 

(IV)  The  partial  specific  volumes  of  the  solvent  and  polymer  are  independent  of  composition  i.e. 
there  is  no  volume  change  on  mixing  and  thermal  expansion  coefficients  are  given  by  their 
average  values  over  the  temperature  range  considered. 

(V)  As  the  temperature  is  increased,  the  increase  in  volume  comes  from  the  homogeneous 
expansion  of  the  material  due  to  the  increasing  amplitude  of  the  anharmonic  vibrations  with 
temperature  and  from  the  formation  of  holes  or  vacancies  which  are  distributed  discontinuously 
throughout  the  material  at  any  instant.  The  free  space  that  is  not  occupied  by  the  constituent 
molecules  and  is  not  associated  with  any  holes  or  vacancies  is  called  interstitial  free  volume.  The 
energy  for  redistribution  of  the  interstitial  free  volume  is  very  large  and  this  type  of  free  volume 
therefore  must  be  distributed  uniformly  among  the  molecules.  The  free  volume  associated  with 
the  discontinuous  distribution  of  holes  in  the  liquid  is  called  the  hole  free  volume.  This  free  volume 
can  be  redistributed  with  no  increase  in  energy. 


In  this  study  the  temperature  dependence  will  be  analyzed  using  the  WLF  equation  for 
diffusion(20),  which  is  also  commonly  placed  in  the  context  of  free  volume  theories.  The  well 
known  form  of  the  WLF  equation  is 


•ogio-f  = 

bO 


-Ci(T-TJ 

C2+T-T0 


(5) 


where  C  is  the  translational  friction  coefficient  of  a  penetrant  molecule  into  the  polymer  matrix,  Co 

is  the  value  of  C  at  Tq  and  ci  D  and  02°  are  the  WLF  coefficients  at  the  reference  temperature  Tq. 

The  diffusion  coefficient  ,D  .can  be  Introduced  into  this  equation  by  using  the  well  known 
expression : 


D  = 


kT 

C 


(6) 


where  k  is  the  Boltzmann's  constant.  The  WLF  equation  can  thus  be  applied  to  diffusion. 


The  local  segmental  motion  of  the  polymer  and  the  reorientational  motion  of  the 
penetrant  /  solvent  are  determined  using  Ti  NMR  measurements.  Previous  studies(2>21) 
have  shown  that  segmental  motion  and  diffusion  of  a  penetrant  in  PIB  may  be  correlated,  which  is 
not  generally  the  case  for  other  polymers.  This  was  attributed  to  the  very  efficient  molecular 
packing  of  the  polymer(‘').  This  high  density  of  PIB  relative  to  many  other  polymers  implies  that 
chain  reorientation  must  accompany  translational  motion  of  the  penetrant.  The  objective  of  this 
work  is  to  compare  the  time  scale  of  segmental  motion  of  the  polymer  with  the  self  -  diffusion 
coefficient  of  the  toluene  in  the  PIB  -  toluene  solutions  as  a  function  of  concentration  to  see  if  they 
are  indeed  correlated.  Segmental  motion  in  pure  PIB  has  been  examined  in  detail  in  a  recent 
report(2)  and  that  information  will  serve  as  a  guide  in  the  interpretation  of  the'  concentration 
dependent  results  to  be  presented  here. 

Experimental 

The  polyisobutylene  sample  used  for  the  study  was  obtained  from  Cellomer  Associates, 
Inc.  of  Webster,  NY  in  the  form  of  thick  sheets  of  high  molecular  weight  (Cat.  #40E,  Lot  #02,  mol. 
wt  1,000,000).  HPLC  Grade  toluene  was  obtained  from  RsherChemical  Fisher  Scientific.  The 
diffusion  measurements  were  performed  on  a  Bruker  MSL  -  300  spectrometer  equipped  with  a 
High  Resolution  Diffusion  Unit  and  gradient  probe,  with  a  proton  frequency  of  300.1  MHz.  For  the 
diffusion  measurements  the  phenyl  proton  peak  of  toluene  was  observed. 

The  well  known  Stejskal  -  Tanner  pulse  sequence(22)  %i  2  -  G(5)  -  Jt  -  G(8)  -  (echo)  was 
used  to  measure  the  diffusion  coefficient  of  toluene  in  PIB  solutions.  The  stimulated  echoC^^) 
pulse  sequence  n  /  2  -  G(5)  -  jc  /  2  -  x  -  jc  /  2  -  G(5)  -  (echo)  was  also  used  for  measuring  the 
diffusion  coefficient,  employing  longer  diffusion  time  observations.  The  strength  of  the  magnetic 
field  gradients  was  varied  from  14  gauss  /  cm  up  to  140  gauss  /  cm;  the  associated  length  5,  from 
500  psec  up  to  7  msec;  and  the  distance  between  the  two  gradient  pulses  A,  from  2  msec  up  to 
80  msec.  Those  values  of  A  correspond  to  diffusion  distances  from  approximately  1.1  pm  up  to 
9.1  pm. 


NMR  spin  -  lattice  relaxation  (T-| )  experiments  were  carried  out  on  Bruker  AC-200,  MSL- 
300  and  Varian  Unity  500  spectrometers,  which  correspond  to  50.33  MHz,  75.5  MHz  and  125.7 
MHz  "*30  Larmor  frequencies  respectively. 

Interpretation 

The  temperature  dependence  of  toluene  in  PIB  was  analyzed  with  the  following  form  of 
the  WLF  equation(20): 
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Dr.T  Cj  +  T-T„  (7) 

nmi  is.*™ 

To  '  self  diffusion  coefficient  at  Tn 

Equation  (7)  can  be  rewritten  as: 


a((|>J 


(T-TJ 


B.  ^  B.  -To)) 


The  pammeters  thal  were  floaled  ir,  order  to  fit  the  data  were  a(M  /  Bh  and  ffo  1  /  Ra  n  • 
analogous  to  the  B  in  equation  (1)  and  wili  be  defined  precisely  later  Tbe  f,  , 

experimental  self,  diffusion  coefficients  at  dillaii.nr  .  ■  ^  «  hts  along  with  the 

i"  Figure  , ,  and  the  reference  temperature  Tg  was  c^nrorgs  K  ^hTT'T 
ere  spawn  in  Table  , .  Mer  ofdaVBq  a.  «rBp  were  obtaL  X  “  oT^ 
02°  were  calculated  (from  ci  D  =  bh  /  2  3n-:i  fr,K  a  ^4  n  coefficients  ci  u  and 

shown  in  Table  1.  '  “  <^f))  end  they  are  also 

volume  theory(12)'*was*;^Tplted'  soncentratlon,  the  Fujita  free 

Doolittleds  •  18)  a,pa  ^  ='«''e,  the  theory  begins  with  the  result  of 

9-n  by  equation  -ight  is 

With  ^hiciti" r 

mo,ecu,arfrictioncoe«entandtheviscosityqZ^^^^^^^ 

ln(m)=A-Bd/f 

(8) 

.■r:r;  ;ir:, 


D  =  RTm 


(9) 


where  R  is  the  gas  constant.  T  is  the  absolute  temperature  and  D  is  the  self  -  diffusion  coefficient. 

In  general  the  free  volume  of  a  given  polymer  -  penetrant  solution  depends  on 
temperature  and  concentration  and  it  is  more  appropriate  to  denote  the  fractional  free  volume  as 
f(T,  (|)s)  where  again  (frs  is  the  volume  fraction  of  the  penetrant..  It  has  been  shown(25)  that  if  the 
increase  in  free  volume  by  the  addHion  of  the  diluent  is  proportional  to  the  volume  of  the  added 
penetrant,  then: 

vf  =  (vf  )2  +  yCO  v-i 


where  vf  is  the  average  free  volume  of  the  solution,  (vf  )2  is  the  average  free  volume  of  the  pure 
polymer.  yO)  is  a  proportionality  factor  that  may  be  compared  with  the  fractional  free  volume  of 
the  diluent,  v  is  the  total  volume  of  the  solution,  vi  is  the  volume  of  the  penetrant  and  V2  is  the 
volume  of  the  amorphous  polymer.  By  dividing  equation  (10)  with  equation  (11)  the  fractional  free 
volume  of  the  solution  f(T.  ^s)  's  found  to  be  a  linear  function  of  the  volume  fraction  of  the 

penetrant  <|)s: 

f(T,  (|>s)  =  fp(T.  0)  +  ( fs(T)  -  fp(T.  0) ) 

where  fp(T.  0)  is  the  fractional  free  volume  of  pure  polymer  and  fs(T)  is  a  parameter  that 
represents  the  fractional  free  volume  of  the  solvent  trapped  in  the  polymer  matrix.  fs(T)  -  fp(T.  0) 
represents  the  effectiveness  of  the  penetrant  for  increasing  the  free  volume  when  it  is  dispersed 
in  the  given  polymer.  Equation  (12)  can  be  rewritten  as: 

_Mn  +  (MI)_!EEl)(S  (13) 

“b:  B,  B,  B, 


If  equations  (12)  and  (8)  are  substituted  into  equation  (9)  considering  that  as  (()s  ->  0.  D  ->  Dq  we 
obtain: 


D  _  Bd(fs-fp)(l>s 
Do  il+ih-Ws 


(14) 


where  <|,s  is  the  penetrant  volume  fraction.  Dq  is  the  diffusion  coefficient  of  toluene  at  the  limrt  of 
zero  penetrant  concentration,  fp  is  the  average  fractional  free  volume  of  pure  polymer  and  fs  is 
the  fractional  free  volume  of  toluene  in  the  solution. 
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For  each  temperature,  the  diffusion  coefficient  of  toluene  as  a  function  of  concentration 
was  fitted  using  equation  (14).  The  parameters  varied  are  Dq,  (fp  /  Bj)  and  (fs  /  Bd).  The  reason 
for  using  (fp  /  Bd)  and  (fs  /  Bd)  as  fitting  parameters  and  not  fp  and  fs  is  to  avoid  making  any 
assumptions  about  the  parameter  Bd-  The  parameters  (fp  /  Bd)  and  (fs  /  Bd)  do  not  have  physical 
significance  whereas  fp,  fs  and  Bd  do.  The  parameters  produced  by  the  fit  are  shown  in  Table  2 
and  the  lines  corresponding  to  the  fit  in  Figures  2  and  3. 


From  the  diffusion  coefficient  of  toluene  at  the  limit  of  zero  penetrant  concentration  in 
Table  2  and  the  WLF  equation  (7),  WLF  parameters  can  be  obtained  similar  to  those  shown  in 
Table  1 .  Those  are  shown  in  the  Table  3 . 


The  self  -  diffusion  of  toluene  in  PIB  was  also  analysed  via  the  Vrentas  -  Duda  theory(^^' 
15,26-28  )  According  to  this  theory  the  self  diffusion  coefficient  of  a  penetrant  into  a  rubbery 
polymer  is  given  by: 


D,  =  Do  - 

RT  -  T„  +  T) + 0),  (!^)(K2,  -  T„  +  T) 


]  (15) 


where  Vj  is  the  specific  critical  hole  free  volume  of  component  i  required  for  a  jump  (i  can  be 
equal  to  1  or  2;  1  for  penetrant  or  solvent  and  2  for  polymer),  (oj  is  the  mass  fraction  of  component 
i,  Tgi  is  the  glass  transition  temperature  of  component  i.  Do  is  a  constant  pre-exponential  factor,  E 
is  the  energy  per  mole  that  a  molecule  needs  to  overcome  attractive  forces  which  constrain  it  to 
its  neighbors  and  y  is  an  overlap  factor  (between  1/2  and  1)  which  is  introduced  because  the 
same  free  volume  Is  available  to  more  than  one  molecule.  K-i  -j  and  K21  are  free-volume 
parameters  for  the  solvent,  while  K12  arid  K22  are  free  volume  parameters  for  the  polymer. 


The  free  volume  parameters  Ki  1 ,  K21 ,  Ki  2.  K22  are  defined  as  follows(1 5): 


K,,  =  V?(Tg,)[a,-(1-f°)a,J 
K 

a,-(1-fS,)a,, 


(16) 

(17) 


Ki2=V°(Tg2)[a2-(1-f°2)a,2] 


'H2 


^H2)®c2 


(18) 


(19) 


where  aj  is  the  thermal  expansion  coefficient  for  the  equiiibrium  iiquid  component  i,  aci  is  the 
thermal  expansion  coefficient  for  the  component  i  for  the  sum  of  the  specific  occupied  volume  and 
the  specific  interstitial  volume,  VjO  is  the  specific  volume  of  the  pure  equilibrium  component  i  and 
fHi°  is  the  fractional  hole  free  volume  of  the  component  i  at  Tgj.  Finally,  §  is  the  ratio  of  molar 
volume  of  a  solvent  jumping  unit  to  the  molar  volume  of  a  polymer  jumping  unit.  The  parameter  ^ 
can  be  a  measure  of  the  molar  volume  of  the  polymer  jumping  unit  when  the  entire  solvent 
molecule  performs  a  jump  during  the  process  of  the  self  -  diffusion. 

There  are  ninejndependent  parameters  in  equation  (15),  Dq,  E,  Kn/y,  K21  -  Tgi 
Kl2/y.  K22  *  Tg2,  V.J  ,  V2  that  need  to  be  determined.  However,  most  of  them  can  be 
determined  from  viscosity  data  or  sources  other  than  diffusion  data  and  in  favorable  cases  all  of 
them  could  conceivably  be  set  from  other  data.  In  order  to  fit  the  diffusion  data  of  toluene  in  PIB 
as  a  function  of  concentration,  the  parameters  shown  in  Table  4  were  used.  This  parameter  set 
was  chosen  from  a  number  of  literature  sources  and  is  discussed  in  some  detail  elsewhere[29]. 
Since  all  the  toluene  and  polymer  free  volume  parameters  mentioned  above  are  close  to  each 
other,  the  actual  choice  of  a  parameter  set  is  somewhat  arbitrary  and  doesn't  significantly  change 
the  quality  of  the  fit.  As  a  result  of  the  fit,  the  only  parameter  floated,  was  found  to  be  equal  to 
0.892  ±  0.08.  The  expenmental  self  -  diffusion  coefficients  of  toluene  in  PIB  along  with  the  fitting 

curves  (solid  lines)  are  shown  as  a  function  of  temperature  and  concentration  in  Figures  4(a)  and 
4(b). 


The  polymer  segmental  motion  was  characterized  using  the  j.j 
methylene  carbon  of  PIB  and  applying  the  model  based  on  the  Hall-Helfand  description[30]  of 
segmental  motion  combined  with  restricted  anisotropic  rotational  diffusion[2.31].  jhe  equations 

that  relate  the  spectral  density  function  corresponding  to  this  motional  model  and  the  spin  -  lattice 
relaxation  time  are  given  by: 


where 


:p  =  W„+2W„  +  W, 

'  1 


(22) 
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W„=S3y|y2S»^ 

The  methylene  C  -  H  distance  is  set  at  1 .09  A.  The  expression  for  the  spectral  density 
function  is  described  in  detail  elsewheret2>29].  According  to  this  model  two  correlation  times 
describe  the  segmental  motion  of  the  polymer.  One  of  them  is  tq  which  is  the  correlation  time  for 
single  conformational  transitions  and  the  other  is  the  correlation  time  for  cooperative  or 
correlated  transitions.  Each  of  these  two  correlation  times  was  given  an  Arrhenius  temperature 
dependence  with  activation  energy  EaO  and  Eai  respectively  and  corresponding  prefactors  Xoo. 
Also  £  is  the  amplitude  over  which  restricted  rotational  diffusion  occurs  and  Djr  is  the  rotational 
diffusion  constant.  The  form  of  the  temperature  dependence  of  £  and  Djr  is  shown  below: 


^  =  A  1^-5  deg 

(23) 

Djr  =  B  T  -  C  sec-1 

(24) 

The  results  of  fitting  the  T 1  field,  temperature  and  concentration  dependence  is  shown  for 
a  typical  case  in  Figure  5  and  the  results  are  shown  in  Table  5  while  the  rest  of  the  parameters 
have  the  values: 


(^00)0  =  2.1  1 0’l  ^  sec;  (xoo)  1  =  1 .0  1 0*^  ^  sec  (25) 

^  =  5.1 5  T<>-5  deg;  Djr  =  3.26  1 0^  T  -  3.5  1 0^  sec-1  (26) 


When  a  WLF  temperature  dependence  was  given  to  the  correlation  times,  there  are  six 
parameters  to  be  established:  (to)o.  (ti)o.  Djr,  c-j®  and  C2®  (c-j®  and  C2®  are  the  WLF 

coefficients).  The  reference  temperature  was  chosen  to  be  Tq  =  298  K  and  accordingly: 


To  _-2.303c?(T-To) 

('Co)o  c°+T-To 

T,  _-2.303c?(T-To) 

(xjo  c°+T-To 


(27) 


(28) 


c°  ^0 
'  2.303f((j)3) 

where  f(<l)s)  is  the  fractional  free  volume  of  the  solution  and  a  is  the  thermal  expansion 
coefficient  of  the  fractional  free  volume.  Bq  represents  the  minimum  void  size  necessary  for  a 
conformational  transition.  Comparable  fittings  to  those  obtained  using  the  Arrhenius  depedence 
are  obtained  and  the  fitting  parameters  are  given  in  Table  6.  Good  agreement  is  evident  for  the 
fractional  free  volume  of  the  solution  derived  from  equation  29  and  those  values  obtained  from 
penetrant  self  diffusion  analysis  shown  in  Table  1 . 

There  is  no  field  strength  dependence  in  the  T-|  data  for  the  ortho  carbon  of  toluene. 
Therefore,  the  extreme  narrowing  limit  equatfon  is  used  for  the  interpretation: 

=  fYHY§^'S(S  +  1)T^  (31) 

where  S  =  1  /  2,  t  is  the  correlation  time  for  molecular  tumbling  and  rcH  =  1 .08  A  is  the  distance 
between  the  proton  and  the  carbon  on  the  toluene  phenyl  ring. 

Again,  at  first  an  Arrhenius  temperature  dependence  was  given  to  the  correlation  time, 
where  E  is  the  activation  energy  for  rotation  and  -Coo  is  a  prefactor.  Two  parameters  were  floated  E 
and  Too.  During  the  fitting  procedure  it  was  found  that  Too  is  around  1.0x10'14  sec  and  changes 
little  as  a  function  of  concentration.  Thus  the  value  of  Too  was  set  equal  to  1 .0  10'^'^  sec  and  E 

was  the  only  parameter  that  was  floated  as  concentration  changed.  The  parameters  obtained  are 
shown  in  Table  7. 


(29) 

(30) 


Next  a  WLF  temperature  dependence  was  assigned  to  the  correlation  time 


„^t_-2.303c?(T-T„) 
to  c°  +  T-To 


(32) 


where  tq  is  the  value  of  the  correlation  time  t  at  the  reference  temperature  Tq  which  again  was 
chosen  to  be  298  K.  ci  ®  and  C2®  are  the  WLF  coefficients  at  the  reference  temperature: 


B 

2.303%) 


(33) 
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(34) 


where  f(<j»s)  is  the  fractional  free  volume  of  the  solution  and  a  is  the  thermal  expansion 
coefficient  of  the  fractional  free  volume  of  the  solution.  B  represents  the  minimal  void  size 
necessary  for  rotation.  For  the  fitting  process  of  the  Ti  data  C2®  was  set  equal  to  C2*^  the  WLF 
coefficient  obtained  from  the  temperature  dependence  of  the  self  -  diffusion  of  toluene  in  PIB  - 
toluene  solutions  (Table  1).  The  parameters  that  were  floated  in  order  to  fit  the  data  are  c-|®  and 
TQ.  the  values  obtained  are  shown  in  Table  8.  The  fittings  obtained  are  comparable  to  those 
obtained  for  the  Arrhenius  temperature  dependence. 


Discussion 


It  Is  clear  that  the  self-diffusion  of  toluene  in  PIB  is  well  described  by  free  volume 
theories.  Some  comparisons  between  the  different  free  volume  treatments  and  the  potential 
correlation  between  penetrant  diffusion  and  polymer  segmental  motion  can  now  be  explored  with 
the  data  in  hand. 


As  a  first  step  a  comparison  can  be  made  between  the  fractional  free  volume  determined 
from  the  temperature  dependence  of  the  translational  diffusion  of  the  toluene  with  that  determined 
from  the  concentration  dependence  of  the  same  quantity.  Such  a  comparison  is  made  in  Figure  6 
where  the  fractional  free  volume  of  the  solution,  f(T,  <|>s),  is  presented  as  a  function  of  the  volume 
fraction  of  toluene  in  PIB.  The  points  come  from  Table  1  and  they  were  obtained  from  the  WLF 
equation  (7)  where  the  temperature  dependence  of  the  toluene  self  -  diffusion  was  analyzed, 
except  the  filled  circle  that  represents  the  fractional  free  volume  of  pure  toluene  (at  293  K) 
obtained  from  viscosity  data(23).  The  solid  line  was  obtained  from  equation  (13)  using  the  values 
of  Table  2  for  fp(T)  and  fs(T)  (obtained  from  equation  (14)  where  the  concentration  dependence 
of  the  toluene  self  -  diffusion  was  analyzed).  The  error  bars  represent  an  uncertainty  of 
approximately  20%  for  f(T,  (t)s)  in  Table  1,  which  comes  from  quality  of  the  fitting  (namely  the 
fitting  of  the  temperature  dependence  of  the  self  -  diffusion  of  toluene  in  PIB).  As  can  be  seen 
from  Rgure  6,  the  values  of  the  fractional  free  volume  of  the  solution  that  are  obtained  from  the 
analysis  of  the  temperature  and  concentration  dependencies  of  the  self  -  diffusion  of  toluene  in 
PIB,  are,  in  fact,  close  to  each  other. 

A  free  volume  analysis  of  the  viscosity  of  the  toluene  -  PIB  solution  was  done  by 
Fujita(32).  |n  that  study  the  fractional  free  volume  of  the  solution,  f(293  K,  <|)s),  and  the  fractional 
free  volume  of  toluene  trapped  in  the  polymer  matrix,  fs(293  K),  are  reported  as  a  function  of 
concentration  at  293  K.  The  fractional  free  volumes  of  the  solutions  at  293  K  obtained  from  the 


self  -  diffusion  data  are  in  reasonable  agreement  wUh  the  values  obtained  from  the  viscosity 
datal  I.  This  establishes  the  utility  of  the  tree  volume  approach  in  relation  to  molecular  mobility 
of  a  penetrant  in  a  polymer.  Viscosity  is  macroscopic  measure,  while  the  diffusion  measurements 
presented  here  are  microscopic  measurements  looking  specifically  at  the  penetrant  molecule 
Two  drfferent  experiments  probing  two  diflerent  levels  of  dynamics  give  approximately  the  same 
fractional  free  volume  for  the  solutions,  when  the  appropriate  value  for  the  minimum  void  size 
necessary  for  diffusion  relative  to  viscosity  was  chosen. 

Another  result  that  has  been  confirmed,  is  the  decrease  of  the  fractional  free  volume  of 
toluene  when  the  latter  is  trapped  in  entangled  PIB.  This  is  seen  by  direct  comparison  between 

r  2)  and  the  data  for  neat  toluene  from  viscosity 

datal  .  J  The  same  conclusion  has  been  reached  from  the  viscosity  data(32).  with  a  completely 

different  approach.  Also  the  fractional  free  volume  parameters  for  pure  PIB  obtained  from  this 
work  are  similar  to  the  values  reported  in  the  literature  and  extracted  from  other  diffusion  or 
viscoelastic  measurements.  This  confirms  the  correctness  of  the  magnitude  of  the  free  volume  of 
the  particular  polymer.From  the  preceding  it  is  obvious  that  Bd.  the  parameter  that  represents  the 
minimum  void  size  necessary  for  a  penetrant  molecule  to  diffuse,  Is  important  if  results  obtained 
from  different  methods  that  observe  different  kind  of  motions  are  to  be  compared.  This  is  to  be 
expected  since  the  size  of  a  moving  unit  is  not  the  same  for  different  kind  of  motions.  The 
unfortunate  fact  is  that  it  is  not  possible  to  obtain  the  dimension  of  the  void  when  the  parameter 


The  stated  advantage  of  the  Ventras-Duda  theory  Is  that  the  parameters  can  be 
independently  determined  and  the  theory  has  potential  predictive  capability.  In  our  use  of 
Ventras-Duda  the  only  fitting  parameter  is  the  ratio  of  the  molar  volume  of  solvent  jumping  unit 
to  the  molar  volume  of  polymer  jumping  unit,  given  by: 

g  -  Vi°(0)  _  MX 

The  value  determined  for  4  of  0.892  ±  0.08.  corresponds  to  a  polymer  jumping  unit  of 
molecular  weight  94.3  as  derived  from  equation  (35)  when  the  toluene  moves  as  a  single  unit. 
This  value  Is  close  to  the  molecular  weight  of  toluene  (Mi  =  92.13).  It  is  thus  seen  that  the  size  of 
the  polymer  jumping  unit  Is  approximately  the  same  as  the  size  of  the  penetrant  jumping  unit.  In 
view  of  this  result.  It  would  be  appropriate  to  expect  that  both  polymer  and  penetrant  jumping 
units  would  require  approximately  the  same  minimum  hole  size  necessary  for  a  diffusion  jump  It 
has  been  shown(13)  previously  that  both  Fujita  and  Vrentas  -  Duda  free  volume  theories  can  be 


used  when  the  molecular  weights  of  the  two  jumping  units  are  the  same.  As  indicated  in  reference 
(33),  the  two  theories  are  identical  when  the  following  condition  is  true: 


V°  V? 


(36) 


and  when  Bd  (minimum  hole  size  required  for  a  given  molecule  to  permit  a  displacement, 
according  to  Fujita  theoryl'^2])  is  defined  as: 


(37) 


This  can  be  shown  by  rewriting  equation  (15)  as: 


4>i(- 


V°Y  V?Y 


) 


Di{0)  k(ll_k)] 

Y  Y  Y 


(38) 


where  Di  (0)  is  the  self  -  diffusion  coefficient  of  the  penetrant  at  the  limit  cot  =  0,  Vj®  is  the  specific 
volume  of  pure  component  i  at  the  temperature  of  interest,  Vj*  is  the  specific  critical  hole  free 
volume  of  component  i  required  for  a  jump,  (|)1  is  the  volume  fraction  of  the  solvent  and  fj  is  the 
fractional  hole  -  free  volume  of  pure  component  i  at  the  temperature  of  interest  and  the  fj's  are 
given  by: 


(K2i  +  T-T,,)K,, 

V? 


{K^  +  T-T^)K,, 


(39) 


Using  equations  (36)  and  (37),  equation  (38)  can  be  written  as: 

~^2 


(40) 


Equation  (40)  is  identical  to  the  equation  derived  by  Fujita (^2).  For  ^  =  0.892,  V2  =  1 .004 
cm^  /  gr,  V-|  *  =  0.91 7  cm^  /  gr,  V2®  =  1 .091  cm®  /  gr,  and  Vi  ®  =  1 .1 54  cm®  /  gr  we  find  that  the 
left  hand  side  of  equation  (36)  is  equal  to  0.821  and  the  right  hand  side  of  equation  (36)  is  equal 
to  0.794.  Therefore  we  see  that  since  the  two  ratios  are  close,  either  of  the  two  theories  can  be 
used  for  the  toluene  -  polyisobutylene  system.  Thus  it  is  no  surprise  that  both  approaches 
accurately  describe  the  temperature  and  concentration  dependence  of  the  self  -  diffusion 
constant. 


It  would  be  desirable  if  the  ideas  from  these  free  volume  concepts  were  able  to  shed 
some  light  on  the  rather  unique  behaviour  of  PIB  relative  to  other  polymers.  In  particular  the 
unusual  slow  diffusion  exhibited  by  PIB  with  respect  to  a  variety  of  penetrants  might  be  expected 
to  be  reflected  in  the  value  of  certain  free  volume  parameters.  This  is  seen  in  the  fact  that  the 
WLF  parameters:  ci  and  C2.  are  considerably  different  from  most  other  polymers,  which  can  be 
adequately  represented  by  the  so-called  universal  values.  The  parameter  a,  the  fractional  free 
volume  expansion  factor,  is  considerably  smaller  than  is  typical  for  other  polymers  and  does  give 
some  physical  insight  into  the  origin  of  the  lower  free  volume  at  temperatures  above  the  glass 
transition.  This  parameter  is  derived  directly  from  the  Fujita  analysis.  The  Vrentas-Duda 
parameters,  though  useful  in  their  potential  for  predictability,  do  not  appear  to  yield  any  physical 
insight  Into  uniqueness  of  PIB. 

If  we  now  include  consideration  of  the  relaxation  data,  three  different  kind  of  motions 
were  considered:  the  self  -  diffusion  of  toluene  in  high  molecular  weight  PIB,  the  molecular 
rotation  of  toluene  and  the  segmental  motion  of  the  polymer.  In  general  the  ci®  WLF  coefficient 
obtained  from  these  three  types  of  motion  should  be  different,  because  of  the  different  hole  sizes 
required  for  each  type  of  motion  to  occur,  i.e.,  B's  (in  general  c-jO  =  B/2.303f((l)s)  and  C2°  =  f(<})s)  / 

a).  For  the  case  of  diffusion  we  have  Bd,  for  the  toluene  rotation  B  and  for  the  segmental  motion 
of  the  polymer  Bq. 

In  Table  9,  (f((|)s)  /  Bq)  obtained  from  segmental  motion  and  (f((|»s)  /  Bd)  obtained  from 
translational  diffusion  are  compared.  From  Table  9,  it  can  be  seen  that  the  two  parameters  are 
close  to  each  other,  considering  their  uncertainties:  of  the  order  of  ±  0.015  at  low  concentration, 
up  to  ±  0.05  at  high  toluene  concentrations.  The  agreement  is  very  good  noting  that  the  values 
are  extracted  by  different  methods  from  two  completely  different  sets  of  data  that  describe 
different  components  of  the  solution.  The  toluene  self  -  diffusion  data  were  collected  by  observing 
directly  the  toluene  molecule  and  the  T i  data  were  collected  by  observing  the  methylene  carbon 
of  PIB  in  the  toluene  -  PIB  solutions.  From  Table  9,  since  the  fractional  free  volume  of  the 
solution,  f((t)s)  and  the  thermal  expansion  coefficient  of  the  free  volume,  a,  are  unique  and 
independent  of  the  type  of  motion  that  is  examined.  It  can  be  concluded  that  Bj  =  Bq.  Also  the 
value  of  obtained  from  the  Vrentas  -  Duda  theory  is  0.892  which  corresponds  to  a  polymer 
jumping  unit  with  molecular  weight  94.3.  This  Is  very  nearly  equal  to  the  molecular  weight  of 
toluene  which  Is  92.  The  fact  that  the  two  nearly  coincide  with  each  other  means  that  the  two 
jumping  units  will  need  the  same  minimal  free  volume  for  translation  (for  toluene)  and  a 
conformational  transition  (for  the  polymer)  implying  that  Bd  should  equal  Bq. 


Another  way  to  check  to  see  if  the  two  motions  are  correlated,  is  to  directly  compare  the 
correlation  time  for  translational  diffusion  xq  and  the  correlation  time  for  segmental  motion 
(correlation  time  for  single  conformational  transition).  The  first  can  be  obtained  from(^): 


where  D  is  the  self  -  diffusion  coefficient  of  toluene,  b  is  the  distance  of  closest  approach  which  is 
between  2.2  A  -  3.3  A.  Using  2.7  A  as  the  distance  of  closest  approach  and  the  self  -  diffusion  of 
toluene  in  PIB  -  toluene  solution,  the  correlation  time  for  translational  diffusion  is  calculated  and 
shown  in  Table  10  The  correlation  time  for  a  conformational  transition,  xq,  is  calculated  from  the 
Arrhenius  equation  and  also  shown  in  Table  10  It  can  be  seen  that  the  two  correlation  times  are  in 
agreement  at  the  0.045  toluene  volume  fraction  which  means  at  this  low  penetrant  concentration 
the  two  motions  (segmental  and  translational  diffusion)  are  correlated.  As  we  go  to  higher 
concentrations  the  difference  between  the  xq  correlation  time  and  the  xq  becomes  larger  up  to 
two  orders  of  magnitude  at  0.71 2  toluene  volume  fraction.  This  might  be  expected  since  as  we  go 
towards  higher  concentrations  a  toluene  molecule  has  fewer  polymer  molecules  in  its 
neighbourhood  and  it  is  not  necessary  for  a  polymer  segment  to  move  in  order  a  toluene 
molecule  to  translate. 

Thus  we  see  that  the  two  processes  (diffusion  and  conformational  transition)  require 
similar  minimum  free  volumes  in  order  to  occur  and  the  time  scale  of  the  two  motions  are 
correlated  at  low  concentratbns  of  penetrant. 

The  spin  -  lattice  relaxation  data  on  the  ortho  carbon  of  toluene  do  not  provide  as  much 
information  about  the  motion  of  toluene  as  the  corresponding  data  on  PIB.  However,  at  all 
concentrations  the  correlation  time  for  toluene  reorientation  is  two  or  more  orders  of  magnitude 
faster  than  segmental  motion  and  translational  diffusion.  Thus  the  reorientational  motion  of  the 
toluene  is  uncoupled  from  the  other  two  motions.  The  WLF  analysis  of  the  temperature 
dependence  also  leads  to  different  values  of  c-j  ®  and  thus  B  further  supporting  the  separation  of 
these  motions.  Toluene  reorientation  would  appear  to  be  a  much  more  rapid  motion  involving  a 
smaller  hole  size. 
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Table  1 

Parameters  obtained  from  the  temperature  dependence  fit  of  the  toluene  self  -  diffusion 

coefficient,  using  equation  (6), 


<l>s 

mfBd 

a(<|>s) !  Bd  (deg-"*) 

c-|B 

C2^(K) 

0.045 

0.0904 

5.8310-4 

4.803 

155.060 

0.104 

0.1050 

6.8410-4 

4.135 

153.509 

0.205 

0.1270 

7.2510-4 

3.419 

175.172 

0.327 

0.1550 

8.4410-4 

2.801 

183.649 

0.409 

0.1900 

9.6010-4 

2.285 

197.917 

0.712 

0.2700 

1.1010-3 

1.608 

245.455 

1.000 

0.3700 

1.7010-3 

1.174 

217.647 
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Table  2 


Parameters  obtained  from  the  concentration  dependence  fit 

coefficient,  using  equation  (14). 

of  the  toluene  self  -  diffusion 

T(K) 

fs/Bd 

fp/Bd 

Do  (cm2  /  sec) 

298 

0.273 

0.08735 

9.0810-9 

310 

0.300 

0.09780 

2.5510-9 

320 

0.333 

0.10700 

5.1510-8 

329 

0.364 

0.11630 

8.7010-8 

339 

0.381 

0.12500 

1.5510-7 

349 

0.398 

0.13200 

2.3010-7 

358 

0.423 

0.14100 

3.4010-7 

368 

0.449 

0.14450 

4.1010-7 

Table  3 

Parameters  obtained  from  the  temperature  dependence  fit  of  the  toluene  self  -  diffusion 
coefficient  at  zero  toluene  concentration  (Dq,  shown  in  Table  2),  using  equation  (7). 

(|>s  f{<j>s)/Bd  a((t)s)/Bd  (deg‘1)  ci^  02*^  (K) 


0.000 


0.0818 


5.3410-4 


5.308 


153.184 


Table  4 


Values  of  parameters  for  PIB  /  Toluene  systems. 


Vi*(cm3/gr) 

0.917 

V2*(cm3/gr) 

1.004 

Kii  /y(cm3/grK) 

1.45  10-3 

Ki2/y(cm3/grK) 

4.351  10-4 

K21  -Tgi  (K) 

-86.32 

K22*Tg2(K) 

-97.6 

0.892  ±  0.08 

Do  (cm2  /  sec) 

4.82  10-4 

E  (kcal  /  mole) 

0 

Mf 

92.13 

Vi®  (cm^/gr  mole) 

84.48 

Tgi  (K) 

117 

M2mono 

56.04 

0129 

9.33 

0229  (K) 

107.4 

Tg2(K) 

205 
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Table  5 

Activation  energies  for  conformational  transitions  as  a  function  of  concentration  of  toluene  in  PIB. 


Volume  Fraction  of  toluene  in 

PIB. 

EaO  (kJ  /  mole)  . 

Eai  (kJ/mole) 

0.045 

34.0 

34.5 

0.104 

33.4 

33.7 

0.205 

32.7 

32.0 

0.327 

31.9 

30.0 

0.409 

31.7 

28.4 

0.712 

30.4 

25.7 

Table  6 

Fitting  parameters  of  the  spin  -  lattice  relaxation  data  when  a  WLF  temperature  dependence  was 

given  to  the  correlation  times. 


Volume  Fraction 

of  toluene  in  PIB 

(to)o  (sec) 

(xi)o(sec) 

ClO 

C29  (K) 

0.045 

3.0010*8 

1.2010*8 

4.30 

155.06 

0.104 

1.8010*8 

8.40  10*9 

3.80 

153.50 

0.205 

8.85  10*9 

4.25  10*9 

3.40 

175.17 

0.327 

6.2310*9 

1.8810*9 

2.75 

183.65 

0.409 

3.4010*9 

9.70  10*19 

2.65 

197.92 

0.712 

3.0010*9 

4.00  10*19 

1.85 

245.45 
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Table  7 


Fitting  parameters  of  the  ortho  carbon  T-i  data  when  an  Arrhenius  temperature  dependence 

given  to  the  correlation  time. 


Volume  fraction  of  toluene  in 
PIBfe 

Too  (sec) 

E  (kJ  /  mole) 

0.045 

1.010-14 

22.73 

0.104 

1.010-14 

21.84 

0.205 

1.010-14 

20.76 

0.327 

1.010-14 

18.50 

0.409 

1 

o 

q 

16.77 

0.712 

1.010-14 

14.00 

was 


Table  8 

Fitting  parameters  of  the  ortho  carbon  T -i  data  when  a  WLF  temperature  dependence  was  given 

to  the  correlation  time. 


Volume  fraction  of 
toluene  in  PIB  <j)s 

TO  (sec)  at  To  298K 

ci  0  at  298  K 

C2^  (K)at£ 

0.045 

9.72  10-11 

2.45 

155.06 

0.104 

6.80  10-11 

2.35 

153.50 

0.205 

4.10  10-11 

2.25 

175.17 

0.327 

1.6710-11 

2.10 

183.65 

0.409 

8.40  10-12 

2.00 

197.92 

0.712 

2.90  10-12 

1.81 

245.45 
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Table  9 

Fractional  free  volume  comparison  of  the  PIB  -  toluene  solutions  obtained  from  the  diffusion  data 
analysis(''  7)  and  from  the  Ti  data  presented  in  this  work. 


Volume  fraction  of  toluene  in 

PIB  (|)s 

%s) !  Bd 

f(<|)s)  /  Bo 

From  Diffusion  data 

FromT-)  data 

0.045 

0.090 

0.101 

0.104 

0.105 

0.114 

0.205 

0.127 

0.127 

0.327 

0.155 

0.158 

0.409 

0.190 

0.164 

0.712 

0.270 

0.235 

26' 


Table  10 

Comparison  of  the  correlation  times  for  diffusion,  tq,  and  for  single  conformational  transitions  to 

for  0.045  toluene  volume  fraction. 


D  (cm^  /  sec)  Toluene 
Volume  Fraction  = 

T(K) 

TO  (sec) 

TD  (sec) 

0.045 

2.55  10-8 

298 

1.915  10-8 

2.859  10-8 

6.6610-8 

310 

1.126  10-8 

1.095  10-8 

1.2610-7 

320 

7.453  10-9 

5.786  10-9 

2.0410-7 

329 

5.25410-9 

3.574  10-9 

3.11  10-7 

339 

3.642  10-9 

2.344  10-9 

4.55  10-7 

348 

2.666  10-9 

1.602  10-9 

6.31  10-7 

357 

1.982  10-9 

1.155  10-9 

8.1410-7 

367 

1.451  10-9 

8.956  10-10 

Hgure  Captions 


Figure  1:  Self  -  diffusion  coefficient  of  toluene  in  PIB  solutions  as  a  function  of  temperature  at 
different  concentrations.  The  solid  lines  represent  fits  to  the  WLF  diffusion  equation 
presented  in  the  text.  The  zero  penetrant  (0%)  concentration  points  are  from  Table  2. 

Figure  2:  Self  -  diffusion  coefficient  of  toluene  in  PIB  solutions  as  a  function  of  concentration  at 
different  temperatures.  Solid  lines  represent  fits  to  Fujita  free  volume  theory 
presented  in  the  text. 

Figure  3:  Self  -  diffusion  coefficient  of  toluene  in  PIB  solutions  as  a  function  of  concentration  at 
different  temperatures.  Solid  lines  represent  fits  to  Fujita  free  volume  theory 
presented  in  the  text. 

Figure  4(a)  and  4(b):  Diffusion  coefficient  of  toluene  in  PIB  as  a  function  of  concentration  at 
different  temperatures.  Solid  lines  represent  diffusion  coefficients  calculated  from 
equation  (15)  of  Ventras-Duda. 

Figure  5:  T-|  of  the  methylene  carbon  of  PIB  in  toluene  -  PIB  solutions  as  a  function  of 
temperature  and  Larmor  frequency.  The  volume  fraction  of  toluene  in  PIB  is  0.045. 
An  Arrhenius  temperature  dependence  was  given  to  the  correlation  times.  The  solid 
line  corresponds  to  the  fit  for  the  Hall  -  Helfand  function  combined  with  anisotropic 
restricted  rotational  diffusion. 

Rgure  6:  Comparison  of  the  fractional  free  volume  of  the  toluene  -  PIB  solution,  obtained  from  the 
temperature  (points)  and  concentration  (line)  dependence  of  the  diffusion  coefficient. 
Solid  line  represents  values  obtained  from  equation  (13).  Riled  circle  represents  the 
fractional  free  volume  of  pure  toluene  at  293  K(^2)_ 


